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Editorial Notice 


Beginning with the January 1961 issues of the NBS Journal of Research, the covers of 
Journal sections have been redesigned so as to emphasize the subject content of each of the 
four sections. Thus the name of this section becomes RADIO PROPAGATION, Section D 
of the JOURNAL OF RESEARCH, National Bureau of Standards. However, it is recom- 
mended that the present method of citing the Journal remain the same except for the addition 
of the section title so that better continuity can be maintained with the referencing of prior 
issues; see example below. 

For the future the Bureau has adopted the following system as the preferred method of 
referencing. References are to be indicated in the text by the insertion in brackets of the 
author’s name and the year of publication, thus: [Trelease, 1951]. If the author’s name is 
part of the text, only the year is bracketed, thus: “. . . in a paper by Trelease [1951].” If 
there are two or more references citing different papers published in the same year by the same 
author, distinguish them by the letters a, b, c, after the year. If two or more references are 
given in the same brackets they should be given chronologically, thus: [Bean, 1953a and e; 
Harrison and King, 1955b]. At the end of the paper, list all references alphabetically by the 
authors’ names. Include in each entry the following: name of senior author, followed by his 
initials; names of junior authors, each preceded by his initials, title of paper (or book); title of 
publication or journal; volume number (and section title) ; issue number; inclusive page numbers; 
month and year of publication, thus: 


Bean, B. R., and B. A. Cahoon, The use of surface weather observations to predict 
the total atmospheric bending of radio waves at small elevation angles, Proc. 
IRE 45, No. 11, 1545-1546 (Nov. 1957). 

Brite, C. O., C. W. Tolbert, and A. W. Straiton, Radio wave absorption of several 
gases in the 100 to 118 kMe/s frequency range, J. Research NBS 65D, (Radio 
Prop.), No. 1, 17-22 (Jan. 1961). 

Slater, J. C., Microwave transmission, first ed., p. 24 (MeGraw-Hill Book Co., 
Ine., New York, N.Y., 1942). 


The provision of this complete information should be helpful to the reader in the following 
ways: the complete title provides a brief indication of the expected content of the reference; 
the beginning and ending pages indicate the length of the reference; and the volume number, 
issue number, month and year will assist the reader in locating the reference in libraries even in 
those cases where the cited reference contains a typographical or other errors. 

The above described preferred method of referencing has been adopted primarily to aid 
the reader by providing directly in the text information relative to the author and year of the 
reference. This method also is most convenient for the author in preparing his original manu- 
script since it provides a simple and unambiguous means of relating the references in the text 
to the bibliography at the end, and also facilitates the introduction of additional references in 
later drafts. 

In the case of papers in which all references are referred to only once and which also do 
not require a bibliography at the end, the literature references may be given as footnotes on 
the same page. 

Footnote references are to be designated by superscript numbers. In the case of short 
papers involving only a single literature reference, and in such cases as those described above, 
the literature references may be given as footnotes. 


Kennetu A. Norton, Guest Editor. 
November 18, 1960. 
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Additional Comments on 


IMPEDANCE MEASUREMENTS AND 
STANDARDS* 


G. A. Deschamps 


Investigations of magnified and squared SWR responses and the use of 
auxiliary tuners with directional couplers both led to a modified reflectometer 
technique permitting accuracies of better than 0.1 percent in SWR measure- 
ments [Beatty, 1958}. 

The use of auxiliary tuners with a directional coupler also permits the 
construction of a highly accurate phase shifter. Combination of this with the 
above modified reflectometer vields an impedance meter similar to a bridge 
and both magnitude and phase of the reflection coefficient may be independ- 
ently determined to high accuracy [Engen, 1959]. 

Possible errors in the use of a sliding short-circuit as a phase shift standard 
have been discussed and a method for avoiding them described. Accuracies 
of 0.3 degree have been obtained [Magid, 1958]. 

Techniques for the measurement of the equivalent network parameters of 
discontinuities coupling two or more modes in multimode waveguides have 
been described [Felsen et al., 1959; Lewis, 1959]. 

Measurement of standing waves in a specially constructed line was used 
to evaluate high dielectric constants at UHF [Williams and Foster, 1957]. 

Rapid measurement of impedance or transmission properties over a wide 
range of frequencies continues to be of great interest. A system for measuring 
phase and attenuation through a component [Linker and Grimm, 1958] from 
8.7 to 9.6 kMe/s and an automatic impedance plotter [Watts and Alford, 1957] 
covering the frequency ranges 180 to 900 Me/s and 50 to 250 Me/s were 
developed. 

Automatic continuous measurement of phase at microwave can be done 
by shifting the frequency to audiofrequency [Mittra, 1957]. A broadband 
direct reading measuring device has been constructed where the LO tracts the 
RF in order to produce a fixed intermediary frequency [Dropkin, 1958]. 


*This material belongs after the fourth paragraph of the published report (p. 598, J. 
Research NBS, 64D, 6,) and must be used with the references on that page. 
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Incoherent Scattering by Free Electrons as a Technique 
For Studying the Ionosphere and Exosphere: Some 
Observations and Theoretical Considerations’ 


Kenneth L. Bowles 


(September 1, 1960) 


Incoherent scattering by the free electrons of the ionosphere has been suggested as a 
technique for measuring the electron density profile both below and above the F region 


maximum. 


This paper reports observations which confirm the existence of the incoherent 


seatter and show that its intensity is essentially the predicted value. The observed Doppler 


broadening is considerably smaller than originally predicted. 


In the second part of the 


paper, an explanation for the reduced Doppler broadening is offered. The scatter is explained 
as arising from statistical fluctuations of electron density, the distribution of which is controlled 


by the positive ions. 


Introduction 


The ability of a random distribution of free elec- 
trons to exhibit weak scattering of radio waves has 
been cited by Gordon [1] as a possible means of 
measuring electron densities and temperatures con- 
tinuously to heights of several thousand kilometers 
above the earth’s surface. Gordon’s prediction is 
based upon the assumption that the free electrons 
may be considered to scatter mutually independ- 
ently. His development is essentially a combina- 
tion of the Thompson scattering property of the free 
electron with Lord Rayleigh’s theory of the scatter- 
ing of light by the particles of a gas [2]. 

Use of such incoherent scatter to study the ion- 
ized regions of the outer atmosphere should have the 
great virtues of theoretical simplicity, and the 
ability to observe a complete vertical cross section 
of electron density and other parameters continu- 
ously for long periods of time with ground-based 
equipment. With these considerations as motiva- 
tion, and having available a large high-power pulse 
transmitter, the National Bureau of Standards under- 
took experiments to test the existence of the pre- 
dicted scatter and to verify its utility. Results of 
the first preliminary observation were reported by 
the writer in a short communication [3]. Since 
then, more sensitive observations have been made 
that are described in this paper, along with a theo- 
retical interpretation. Much the same material 
was originally circulated among a limited number 
of workers in this field as an unpublished report of 
the National Bureau of Standards [16]. 


! Contribution from Central Radio Propagation Laboratory, National Bureau 
of Standards, Boulder, Colo. 


Part I of the paper is devoted to the experimental 
observations. The characteristics and special fea- 
tures of the equipment are outlined. Observations 
in which the incoherent scatter may be identified 
are illustrated using original oscillograph records. 
The incoherent scattered power is found to be with- 
in a factor of two of the value predicted using 
Gordon’s approach. Representative electron den- 
sity-versus-height profiles are given for two of the 
oscillograms, and shown to agree reasonably well 
with profiles for the lower part of the F region ob- 
tained from conventional ionograms. The possi- 
bility that any one of several other well-known 
propagation modes might have resulted in ambig- 
uous observations is examined and rejected. 

In spite of the apparent verification of Gordon’s 
prediction for echo power, a distinct departure is 
found in terms of the observed echo spectrum. The 
spectral broadening of the scattered echoes is found 
to be considerably less than the Doppler broadening 
predicted using Gordon’s approach. While a spec- 
tral width of approximately 100 ke/s might have been 
expected for F' region echoes, the observation is that 
the transmitted bandwidth of about 9 ke/s must be 
broadened only slightly by the scatter process. 

In Part IT we propose that the scatter should be 
regarded as arising from random fluctuations of 
electron density, rather than from electrons scattering 
individually. In this way we develop a simple ex- 
planation for the slight echo broadening, based upon 
Coulomb interactions between ions and electrons. 

In the past the density fluctuation approach has 
been used successfully to predict the intensity of 
scatter from aggregates of gas particles even in cases 
where the distributions of particle density were not 
entirely random. However, little attention appears 
to have been paid to the question of spectral broaden- 
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ing using density fluctuations. It is generally as- 
sumed that when plasma oscillations can be neg- 
lected, as in the case of high frequency “‘reflections”’ 
from the ionosphere, the spectrum of scatter from 
aggregates of free particles must be directly pro- 
portional to the distribution of line-of-sight veloc- 
ities of the several particles. We offer a simple 
illustration to show that this proportionality may 
not necessarily hold true when there is a degree of 
mutual correlation among the positions of the parti- 
cles. We argue that such correlation is automati- 
cally introduced if there are many particles within a 
zone whose depth corresponds to approximately 180 
degrees of phase or less. This zone defines a pre- 
ferred “‘scale’’ among the density fluctuations, and 
the rate at which irregularities of this scale change 
from one state of fluctuation to another provides 
a means of estimating the spectral breadth of the 
scatter. 

Obviously the picture of density fluctuations 
breaks down when the individual scattering particles 
are so sparsely distributed that only one can be 
found in a volume large relative to the characteristic 
scale. Such a case is found with radar ‘chaff,’ 
where there is every indication that the chaff par- 
ticles contribute to the scatter independently. Evi- 
dently some criterion must be established to deter- 
mine when it is better to consider the scatter as 
arising from density fluctuations, and when from 
independent particles. Two possible “independence 
criteria’ are examined—one based on mean spacing 
between neighboring particles, the other based on 
mean-free-path as proposed by Gordon. The pres- 
ent experiment offers an opportunity to distinguish 
between these criteria. 

It is shown that, apart from a proportionality 
factor close to unity, the spectral width of radiation 
scattered from a randomly distributed gas is the 
same, using either density fluctuations or independ- 
ent particles. From the experimental results, we 
draw the conclusion that the distribution of free 
electrons in the ionosphere is not entirely random 
(even within a volume small enough that the mean 
density gradient can be regarded as negligible). 

We examine the effect of Coulomb interactions 
between ions and electrons on the fluctuations of 
electron density. For scales of irregularities large 
compared with the Debye shielding distance, Xp, 
we find that the fluctuations of ion density and those 
of the electrons are strongly correlated. The mag- 
nitude of the fluctuations is close to that of a neutral 
gas whose density is equal to the electron density. 
Because the ions have so much more momentum, 
the rate of change from one state of fluctuation to 
another is controlled almost exclusively by the ions. 
We therefore estimate, for example, that a spectral 
broadening of the order of 500 c/s would be correct 
for echoes from the F region at our 41-Mc/s operat- 
ing frequency. Since the results of recent experi- 
ments indicate that this estimate is correct, we infer 
that the density fluctuation approach is correct in 
the F region. This suggests that the correct in- 


dependence criterion may be based upon the mean 
spacing between neighboring particles. 

The predicted effects of the earth’s magnetic field 
on the observations are discussed briefly. It is 
pointed out that these effects may permit one to 
identify the ionic constituents and to measure the 
magnetic field strength in the earth’s outer atmos- 
phere. 


Part I. Experimental Observations 
1. Equipment 


Except for the dimensions involved, the equip- 
ment is in most respects a conventional radar 
directed vertically. Operating parameters are sum- 
marized in table 1. 


TaBLE 1. Parameters of radar equipment used 

Operating 40.92 Me/s 

Peak pulse power- - ------ (4 to 6) X 10° w or 1X 10° w (see 
text 

Piise 3 to 150X10-® sec available 
120X10-* sec for all data 
shown in this paper 

Average 4X 104 w maximum 


Pulse repetition frequency. Flexible within above limitations. 
25 to 40 p.p.s. used for most 
work shown in this paper (see 
text) 


Receiver noise (power) 9, 13, or 35 ke/s 
bandwidth 
Receiver noise factor. 2 


Up to approximate factor 100 
(20 db) (see text) 

116 by 140 m (1.6 104m?) 

Linear north-south 


Signal/noise improvement 
by averaging 

Antenna cross section 


Antenna gain, efficiency, See text 
ete. 
NBS Long Branch, Illinois trans- 


mitting station—W WI 
Latitude 40.2° N. 
Longitude 90.0° W. 


The transmitter is a newly built unit of a hitherto 
untried design. Not all malfunctions in its opera- 
tion had been removed at the time the experiments 
here reported were carried out. Absolute measure- 
ments of the incoherent scatter intensity made 
during February 1959 (described in detail later in 
the paper) make it clear that the transmitter was 
generating roughly 1-Mw peak power at the operat- 
ing frequency, rather than 4 to 6 Mw, during the 
work reported in reference 3. The difficulties re- 
sponsible for this malfunctioning have been cor- 
rected, but a new problem appeared before the 
February tests. Due to damage from a severe frost, 
the transmitter could only be operated using its 
driver amplifier with an output power of 1.0 Mw. 
This power was carefully calibrated using a probe 
technique for comparison with the known output of 
a lower power amplifier. Approximately 1 db of 
attenuation was experienced in the temporary con- 
nection of the driver amplifier to the antenna, and 
this is reported below as part of the total loss of 
antenna efficiency. 

The receiver is a communications unit of conven- 
tional design, modified to permit pulse operation. 
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A preamplifier and converter provide low noise fea- 
tures not available in the commercial receiver. 
Bandwidths suitable for pulse work were obtained 
by loading the IF transformers. The switching 
arrangement of the receiver then permitted the 
9-kc/s band to have the characteristic of 6 cascaded 
L-C circuits. The 35-ke/s band had a flat top, 
with skirts roughly the same as those of the 9-kc/s 
band. The bandwidths quoted were measured by 
comparing and equating the output of the receiver 
for calibrated inputs from a noise diode source, then 
from a ew standard signal generator, taking the noise 
statistics and detector operation into account. The 
bandwidths, between frequencies at which the re- 
sponse is down 6 db from the maximum, differ only 
slightly from the noise bandwidths. 

Averaging of the receiver output was used to 
improve the quality of the scatter measurements. 
The averaging device. or “integrator,” is built 
around a multiple-deck mercury-jet commutator. 
Each of 240 independent capacitors is connected to 
the receiver output through a large resistor only 
during an assigned time interval following the trans- 
mitter pulse. Thus the voltage appearing across 
each capacitor represents the RC running average 
of receiver output corresponding to a specific range 
interval—usually 15 to 25 km wide. The vertical 
d-c amplifier of an oscilloscope having high input 
impedance is also connected to the capacitors in 
sequence. The presentation on the oscilloscope is 
therefore a histogram representing the running aver- 
age of many conventional A-’scope sweeps. The 
limitations of the device are such that a maximum 
improvement in effective signal-to-noise ratio of 
about 20 db (a factor of 100) can be obtained. 
Several different time constants were used in taking 


the photographs presented in this report, in each 
case the signal-to-noise ratio improvement being 
somewhat less than 20 db. 

The antenna is a simple broadside array of 1,024 
half-wave dipoles located 0.16 wavelength above 
ground. The dipoles are situated in approximately 
a square array and fed in phase through open wire 
transmission lines. By attaching the _ half-wave 
elements to the transmission line at one wavelength 
intervals, transposition of the lines is avoided. A 
total of 32 elements are connected effectively in 
parallel and fed at the center of each such arrange- 
ment. Those elements located at alternate halt 
wavelength intervals are fed by the adjacent trans- 
mission line. There are a total of 32 such transmis- 
sion lines spaced at half wavelength intervals. The 
neighboring combinations are then connected by 
pairing successively through equal length feeders 
until the whole array is connected to one pair of 
terminals. The reflection coefficient of the ground 
has been increased and stabilized by placing copper 
wires on the ground, parallel to the dipoles. The 
T-R and A-T-R switches are built in open wire line, 
using open tungsten spark gaps in the conventional 
configuration. 

By a fortunate coincidence, a very suitable celestial 
radio source is available as a check on antenna 
performance. The station is located at latitude 
40.2° while the radio source known as Cygnus-A is 
at 40.5° [4]. For practical purposes Cygnus-A can 
therefore be considered as passing through the center 
of the main beam of the antenna. A number of 
total power records of the Cygnus source as it 
transits the beam were taken as a means of measuring 
the east-west antenna pattern and calibrating the 
gain. One such record is shown in figure 1. 


RECEIVER NOISE 


FicurE 1. 


RELATIVE POWER ABOVE! NTERNAL. 


Example of the kind of total-power record obtained using the large dipole-array during transit of the radio 


source Cygnus-A. 
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Although scintillations were often quite severe, 
several passes were obtained when scintillations were 
less than 10 percent of the total power received. 
Hence, a good estimate of the maximum power could 
be made. As may be seen from figure 1, the antenna 
exhibits the behavior expected in the main lobe and 
lowest order side lobes for a uniformly illuminated 
aperture. The measured beam width between half- 
power points is 3.1°, and between first nulls is 7.2°, 
agreeing with the values expected for an aperture 
16 wavelengths wide [5]. It may be noted that there 
appears to be little or no broadening of the trace of 
Cygnus due to the broad angular component known 
as Cygnus-X. Using the best available measure- 
ments of the flux of power from Cygnus-A, those of 
H. W. Wells [6], it is estimated that 1.75<10-" w 
per c/s should be available at the antenna terminals 
at the peak of the Cygnus trace, assuming perfect 
efficiency on one polarization. It was found that an 
equivalent noise diode current of 110 ma through 50 
ohms gives a comparable deflection. Since this 
current generates an input power of approximately 
4.4X10-" w per c/s, the antenna is estimated to be 
25 percent efficient or to have a degradation in gain 
of about 6 db. The accuracy of the measurements 
is certainly not better than 1 db. The loss in effi- 
ciency may be explained by the following three 
factors: 1 db. is lost in coaxial lines and matching 
networks connecting the T—R system to the receiver 
input terminal (or transmitter output terminal). 
Perhaps half of the remainder is represented by 
resistance losses in the poorly screened ground, and 
in the transmission lines, matching networks, etc. 
The remaining efficiency is presumably lost through 
spurious radiation from portions of the open-wire 
transmission lines on which matching was not 
attempted, portions of wires and guys only partially 
removed from the circuit by quarter-wavelength 
shorted stubs, etc. This radiation must appear as 
energy in sidelobes randomly distributed over the 
hemisphere illuminated by the antenna. 


2. Scatter Observations Identified 


The kind of observations obtainable using the 
conventional “A’scope” display is illustrated in 
figure 2._ Cosmic noise is represented by the constant 
“grass” level, peaking at about one-third of the 
vertical scale. The apparent rise in the grass level, 
peaking broadly at about 350 km, is due to the in- 
coherent scatter echo. The sharp peak at about 
80 km appears to arise from the same kind of scatter 
responsible for at least part of oblique path VHF 
scatter propagation [7]. Further observations of 
this D region scatter have been made and it is 
hoped will be discussed in a separate paper. 

The kind of record obtainable using an A’scope 
presentation of the “integrator” output is_ illus- 
trated in figure 3, a direct photograph of the oscil- 
loscope-screen. ‘The cosmic noise level in this case 
is just under the first horizontal line above the lower 
border of the photograph. In this record, and others 
on the same presentation shown below, the broad 
features of the electron-density profile of the F 


region of the ionosphere may be clearly seen. Note 
that while the F region electron density must be 
close to one-tenth the F region maximum, the 
echo intensity varies only slightly between the 
heights of 100 and 350 km. This results from the 
decay of system power sensitivity with the inverse 
square of the height. 

Before positively identifying the observed echoes as 
incoherent scatter, one must reject the possibility 
that they might arise from sidelobe response to some 
other known mode of propagation. Overlapping 
meteor echoes observed via the sidelobes are an obvi- 
ous possibility. Figure 4 illustrates the appearance 
of a typical meteor echo. The sharp spike in the 
midst of the broader but weaker rise in the grass level 
is a meteor echo which lasted a fraction of a second 
compared with the 4-sec exposure time. Figure 5 
shows the integrator presentation for a time close to 
0600 hours local time. As would be expected for 
the F region, the main features of the broad trace are 
considerably reduced in amplitude compared with the 
midday record of figure 3. Between 100- and 200- 
km range, however, there appears to be stronger 
echo intensity with a much rougher distribution, due 
to the greatly increased meteor activity typical of 
times near 0600. Adjacent integrator intervals show 
intensities which vary randomly in an uncorrelated 
fashion. Rejection of meteor echoes as a serious 
contributing factor therefore seems justified, if allow- 
ance is made for their occasional appearance on the 
records, particularly at ranges less than 700 km. 

Long-range F layer-propagated ground-backscatter 
occasionally produces strong interfering echoes in the 
system. Examples of such echoes areshown in figures 
6a and 6b, where slightly more than one complete 
radar cycle is displayed. A second transmitted pulse 
occurs in 6a and 6b at about 5,000-km indicated 
range. Figure 6a is a record obtained in the custom- 
ary manner using the large broadside antenna both 
for transmitting and receiving. Figure 6b was ob- 
tained by substituting, for reception, a simple half- 
wave dipole antenna located one-half wavelength 
above ground. Since the illumination of the scatter- 
ing medium was the same for both examples, it is 
clear that the response of the large antenna at low 
radiation angles approximates the response of the 
dipole within a few decibels. Figures 6c and 6d 
were taken using expanded range and intensity scales 
for making the same comparison. Figure 6c, taken 
with the large antenna used for both transmitting and 
receiving, shows the usual F' region profile. Figure 
6d, taken with the large antenna used for transmit- 
ting, and the dipole for receiving, shows no sign of the 
F region profile. By utilizing this comparison, and 
by taking continual care to adjust the spacing be- 
tween successive transmitter pulses, it has been 
possible to detect and eliminate any contribution to 
the interesting part of the records due to long-range 
backscatter. A similar comparison of the dipole 
with the large antenna has permitted the rejection of 
the possibility of important echo contributions due 
to magnetic field alined irregularities, such as from 


the aurora. 
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Figure 2. A’scope photographs obtained using 4-sec exposure. 
Near midday, October 22, 1958. 
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Ficure 3. A’scope photograph using the “integrator” display. 
1135 1.t., February 27, 1959. 
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Fiaure 4. A’scope photograph showing transient trace due 


to a meteor echo, in addition to incoherent scatter. 
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Figure 5. Integrator A’scope photographs taken at about 


0600 L.t., February 28, 1959. 
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Fievure 6. Long-range backscatter records obtained using: (a) Large array, (b) single dipole—for re- 
ception. Short range, high gain records using: (c) Large array, single dipole—for reception. . 
I 
While it is possible that other interfering modes of 800 
propagation might appear, a convincing final argu- 
ment for the identification of the echoes as incoherent 
scatter is provided by a comparison with known F woo 
region characteristics. Assuming the height of the | 
F region maximum electron density to remain roughly 600 
constant, the maximum scatter intensity should be 
proportional to the square of the F2 region critical 
frequency observed on a conventional vertical in- 500 
cidence sweep frequency sounder. Unfortunately, E 
such a sounder had not been installed at the Long = | N 
Branch, Illinois station at the time of this experi- 3% 400 | T 
ment. The scatter records were compared with w | 
critical frequencies observed on the sweep frequency | ae 
sounder at Fort Belvoir, which is located at about the aes 2 | 
same latitude. Critical frequencies were estimated al 
from the scatter records by normalizing the system 200 sit | 
sensitivity such that at one particular hour the Yy 
maximum indicated electron density corresponded ff | 
with the Fort Belvoir critical frequency for that hour. 100 
The Fort Belvoir critical frequencies were then. found | 
to lie within one-half megacycle of the values calcu- | : 
lated from the scatter records, on the several days in 0 Pt 
February 1959 when comparisons were made. Fig- RELATIVE SCALE 
ure 7 is a replot of figure 3 to account for the height : 
variation of sensitivity, and should be proportional 4 
to the electron density versus height profile. The | Fiaure 7. Electron density profile obtained from A'scope : § Fic 
dotted curves superimposed on this profile are com- display of figure 3. 
puted from the Fort Monmouth and Boulder iono- a D 
grams for the same day and local time using the NBS |  gouider tonograms normalized to the same Monmouth and Bou 
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“‘true height” technique [8]. A second pair of obser- 
vations, made in the same way, are shown in figures 
8 and 9, respectively. This pair, made shortly after 
sundown, shows the rapid decay of ionization below 
F max compared with figures 3 and 7. Although 
neither the agreement of the scatter profile with the 
ionograms, nor the ionograms with each other, is 
exact, the results exhibit a convincing similarity. <A 
more precise comparison has since been made utiliz- 
ing a conventional sounder at the Long Branch 
station [26]. 
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Integrator A’scope photograph taken at 1940 1.t., 
February 27, 1959. 
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Figure 9. Electron density profile obtained from A’scope 
display of figure 8. 


Dotted curves are “true height’’ profiles computed from Fort Monmouth and 
Boulder ionograms, normalized to the same foF 2. 
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3. Comparison of Experimental Observations With 
Theoretical Predictions 


Comparison of the observed scatter power with the 
values predicted theoretically using Gordon’s method 
[1] is possible, using the system parameters as out- 
lined above. Magneto-ionic effects in the iono- 
sphere, which produce Faraday rotation of polariza- 
tion [9], degrade the system sensitivity by3db. This 
may be seen as follows: Somewhat over 50 complete 
rotations of polarization would be expected for a 
41-Mc/s wave passing twice through the entire day- 
time ionosphere. Thus more than one-half revolu- 
tion can be expected over any given interval of 18-km 
height corresponding to the transmitted pulse length. 
The net effect is similar to a randomization of polari- 
zation so that on reception the linearly polarized 
antenna can extract only half of the energy in the 
incident scattered wave. 

Based on all these considerations, the peak scatter 
power at the receiver terminals corresponding to an 
F2 region critical frequency of 12.8 Me/s (210° 
electrons/em*) should be 3.5X10-" w. The height of 
the F2 region maximum electron density is here 
assumed to be 330 km, in good agreement with the 
observed height. The observed scatter power at the 
receiver terminals, corresponding with the F2 region 
maximum, was 4.3107'® w when the critical fre- 
quency was about 13 Mc/s. This must be considered 
as excellent agreement between experiment and 
theory, inasmuch as the overall accuracy of the 
computation cannot be better than plus or minus 3 
db. For receiver bandwidths of both 9 and 35 ke/s 
the observed scatter power was the same. 

In spite of the agreement of the observed echo 
intensity with the prediction, one major difference is 
apparent from the equality of the echo strength for 
9- and 35-ke/s bandwidths. Doppler broadening of 
the echo due to independent scattering from free 
electrons should be about 100 ke/s at F region heights. 
The F region temperature is here assumed to be 
approximately 1500 degrees Kelvin [10]. Since the 
9-kce/s bandwidth apparently admits as much echo 
power as the 35-ke/s bandwidth, the echo must be no 
broader than about 9 ke/s. The spectrum of the 
received energy must be at least 9 ke/s wide since that 
is the spectral width of the transmitted energy. One 
must therefore conclude that the spectral broadening 
due to the scattering process must be considerably 
less than 9 ke/s, otherwise the received spectrum 
would be considerably wider. To further test this 
conclusion, observations were made with the receiver 
slightly detuned from the transmitted frequency. 
One such sequence is shown in figure 10. Figure 10a 
was obtained in the same manner as the records 
shown above with the receiver correctly tuned. 
Figures 10b and 10c were obtained with the receiver 
tuned 15 ke/s above and 15 ke/s below the transmitter 
frequency, respectively. The latter two records 
show no sign of an echo from the F region, leading 
to the conclusion that the bandwidth of the echo 
power must be considerably less than half the re- 
ceiver pass band of 9 ke/s. Since the average ve- 
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Fiacure 10. Integrator A’scope photographs taken: (a) With 
receiver tuned to transmitted frequency, (b) receiver tuned 15 
ke/s high, (c) receiver tuned 15 ke/s low. 


locity of the electrons is expected to be proportional 
to the square root of their kinetic temperature, the 
1500 degree estimate for the F region would have to 
be high by a factor of more than 100 to account for 
these observations on the simple theory. 

Other measures of ionospheric temperature appear 
to be consistent with the 1500 degree estimate well 
within one order of magnitude. The second part 


of this paper is a discussion of the observed echo 
bandwidth. 

Recently several theoretical papers have appeared 
which predict that the intensity of the scatter from 
the ionosphere should be half that predicted by 
Gordon’s method [11, 12, 13, 14, 15]. The reduced 
scatter power predicted by these theories is within 
the experimental accuracy of the present experiment 
and ne disagreement can be inferred. Another 
experiment has also been performed by Pineo and 
others [17] at 440 Me/s, and results very similar to 
ours have been obtained. Their measurement of the 
spectrum of the echoes was much more precise than 
possible with our equipment, and is discussed further 
in Part IT, 3. 


Part II. Modified Theory 
1. Independent Particles Versus Density Fluctuations 


There is a considerable body of literature on the 
subject of scattering of electromagnetic radiation by 
small particles. Following the early work of Lord 
Rayleigh [2], attention has been given to the criterion 
by which individual particles can be considered to 
scatter independently or, using Rayleigh’s term, in 
arbitrary phase. The word “independence” is used 
in this paper in the statistical sense and implies a 
lack of statistical correlation. 

For particle densities greater than those for which 
such independence is correct, calculation of the 
scatter properties has been found to depend upon an 
understanding and description of the statistical 
fluctuations of particle density. The reader is re- 
ferred to Rocard [18] for a comprehensive summary 
and bibliography of the early developments in this 
field. A better statistical description as well as 
a bibliography of more recent work are given by 
De Boer [19]. Using the De Boer approach, the 
writer suggests the following arguments as a basis 
for a modified theory of the incoherent scatter. 

The density fluctuation theory begins with a sta- 
tistical description of the spatial distribution of the 
scattering particles. When neutral particles are 
separated by many particle diameters so that inter- 

article forces are small, a random distribution of 
ocations is established. (This is strictly true, of 
course, only when the gas has no strong gradient of 
density in a distance comparable with the scale size 
under consideration.) The statistical fluctuations 
of particle density create a region containing weak 
irregularities of refractive index which are capable 
of scattering radio waves. When the particle distribu- 
tion can be considered random over all scale sizes 
comparable with the wavelength of the incident radia- 
tion, De Boer [19] finds the scattered intensity to be 
equal to that of the simple independent scattering case. 
De Boer implies that the scattered intensity would 
not depend upon the scale size chosen for the irregu- 
larities. On the other hand, we shall see that a 
specification of scaie size is important to an under- 
standing of the spectral broadening of the scattered 
radiation. 
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On the other hand, when the spatial distribution 
of the electrons is such that the phases of their re- 
spective contributions cannot be considered mutually 
independent, the simple Doppler profile of the previ- 
ous paragraph is no longer correct. This pros ce 
is supported by the fact that the observed frequency 
broadening of the incoherent scatter echo is much 
smaller than that implied by the above Doppler 
profile when any reasonable kinetic temperature is 
assumed for the ionosphere. (As in the previous 
section, we have assumed that the temperature in 
the / region is of the order of 1500 degrees Kelvin.) 
Assuming that the correct explanation of the dis- 
crepancy is to be found in terms of a collective 
description of the spatial distribution of the electrons 
—the fluctuation in density of the electrons due to 
the Brownian motion—an estimate of Doppler 
broadening which is consistent with the observations 
is obtained below. 

A simple illustration will serve to show how a group 
of particles, having a Maxwellian distribution of 
thermal velocities, could under the proper conditions 
produce a resultant having little or no Doppler shift 
by comparison with the broadening which the parti- 
cles would produce, were they independent. Figure 
11a is an idealized sketch of a radar system intended 
for study of the ionosphere. The antenna beam is 
approximated by a uniform response within the 
sector, and zero response outside this range. 

We now consider the relation between the concept 
of density fluctuations as scattering centers and 
the concept of individual electrons as scattering 
centers. In all cases of interest, the incident elec- 
tromagnetic field is modified so slightly by scattering 
that as a first approximation the individual electrons 
must give rise only to single scattering of the incident 
energy. That is, each electron is considered to be 
coupled electromagnetically only to the radar system 
and not to other free electrons. Hence, a detailed 
summation of the contributions of all the free elec- 
trons individually would lead to the correct descrip- 
tion of the scattered field. Such a summation would 
have to account for the phase of the contribution 
from each electron individually—an enormous com- 
putational problem. By taking advantage of ways 
in which large groups of electrons—and their indi- 
vidual contributions—exhibit collective behavior, 
it is possible to simplify the problem considerably. If 
the phase of the contribution from each free electron 
is uniformly distributed throughout a range of at 
least several times 27 radians and is statistically 
independent of the phase of the individual con- 
tribution from each other free electron, the approach 
used by Rayleigh [2] and Gordon [1] is correct. In 
this case the scatter due to each particle appears 
in the resultant signal at the Doppler-shifted fre- 
quency corresponding to its line of sight velocity, 
S—fo=2v/r. The spectrum of the scattered energy 
(arising from a monochromatic source) would have 
the functional form of the probability density func- 
tion describing the number of electrons found at 
each line of sight velocity. Ordinarily this would 
reduce to the usual Gaussian, or Maxwellian [12], 


form 


mip 


S(f)=S(0) exp] — 


where S(f) is the power spectrum function of the 
scattered energy in terms of the spectral frequency, f. 
The transmitted frequency is fo, 0 is the mean speed 
of flight of the electrons, \ is the free space wave- 
length of the radio energy, 7’ is the kinetic tempera- 
ture, k is Boltzmann’s constant, and m, is the mass 
of the electron. 

As shown by the phase vector, or phasor, diagrams 
in figure 11(b), the phase P(A) of an echo arising from 
from a scatterer at a distance A from T is 360 degrees 
advanced relative to the phase P(B) of the echo 
from another scatterer lying a distance B=A-+d/2 
from T. In three dimensions A and B are the radii 
of the spherical surfaces of constant phase with cen- 
ters at T. 


= §(0) exp] — 
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Figure 11. Schematic sketch of radar system. 


Phasors are referred to the transmitted phase—T. 


We now consider our group of particles to be initially 
located at points within the slab AA—CC which has 
thickness d < )/2. The phasor resultant signal 
due to this group of particles is shown in figure 12 
along with the individual components, for represent- 
ative initial particle locations. 

We allow enough time to pass that the particles 
have diffused into a slab of somewhat greater thick- 
ness than AA-CC. Figure 13 is a phasor diagram 
illustrating how the individual components and the 
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Figure 12. Phasor diagram illustrating the initial distribution 
of contributions due to particles within the slab AA-CC. 


~ CC 


Figure 13. Phasor diagram illustrating how the distribution 
of electron contributions of figure 12 may have changed after 
a short time has elapsed. 


resultant phasor might have changed during that 
time. We see that the phases of the individual com- 
ponents have rotated considerably in the new pic- 
ture, but the phase and amplitude of the resultant 
have changed only slightly. Of course, the degree 
to which the resultant changes is a function of the 
distribution of rates of change of the components. 
In the simple illustration above, the reduced rate of 
change of phase of the resultant can be viewed as an 
average of the rates of change of the components. 
The reason for discussing rates of rotation of phasor 
components is, of course, that the shift of frequency 
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of any component relative to the transmitted fre- 
quency is directly proportional to the rate of phase 
rotation of that component. In the simple illustra- 
tion the phases of the components would rapidly 
become randomized as diffusion set in. Until this 
happened, the probability density function describing 
the equivalent spectrum of the resultant signal 
would be markedly reduced in width relative to the 
spectrum of the individual components. This simple 
illustration is analogous to the scatter by free elec- 
trons in an underdense meteor trail. The echo from 
a short length of such a trail is not Doppler broadened 
by thermal motions of the electrons. 

We now ignore, temporarily, the continuous motion 
of the electrons and consider the statistics of the 
signal scattered from a portion of the ionosphere in 
which the electrons are suddenly rendered motionless. 
We consider for two limiting cases the distribution 
of signal strength arising from a slab AA-—CC, as in 
figure 11 above, on a number of successive such 
instants, or ‘‘realizations.” In both cases we con- 
sider that the particles are randomly distributed 
within the slab, all positions having equal probability. 
When, as above, the thickness, d, of the slab is taken 
to be small relative to d/2, or more specifically when 
d</4z, the resultant amplitude, 2, is almost equal 
to the arithmetic sum ot the contributions from the 
individual particles. That is, the resultant echo 
power, P~n?r,, where r, is the scattered echo power 
arising from one electron and n is the total number 
of electrons within the slab. The phase of the 
resultant is somewhere between the phase cor- 
responding to AA and that corresponding to CC, 
depending upon the distribution of particle density 
within the slab. The distribution of the resultant 
signal amplitude over a large number of realizations 
is simply proportional to the distribution of n. By 
virtue of their close spacing in phase, the particles 
in this case contribute to the resultant collectively, 
i.e., coherently. 

The other limiting case is for d>)/2. This is 
the two-dimensional random walk case considered 
by Rayleigh [2], in which 2,,.=nr, and the resultant 
for any given realization is Rayleigh distributed. 


2 R? 
P(R) dR=2,R exp 


where P(R)dR denotes the probability that P lies 
within the range 2 to R+dR. In this case the phase 
of the resultant is uniformly distributed throughout 
2x radians, and without some detailed knowledge of 
the positions of the individual components the phase 
cannot be specified within that range. Judging from 
the result of the argument for a thin slab above, we 
can satisfactorily describe the properties of the signal 
from a thick slab by dividing it into components 
having d¢)/4zm i.e., knowledge of the position of 
every particle is not necessary. This is the import 
of the arguments given by Rocard [18] and De Boer 
[19] which were referred to above. 

Evidently the transition from the collective be- 
havior of the thin slab, to the incoherent behavior 
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which describes the thick slab, is by no means 
sudden. As d increases from }/4m to larger values 
there is a transition region in which the coherent 
aspects of the thin slab are gradually lost. It should 
suffice to say that there should be some maximum 
d for a thin slab to exhibit essentially coherent be- 
havior, and some minimum d for a thick slab to 
exhibit essentially incoherent behavior. We make 
the approximation that for some intermediate value 
of d, the electrons within each slab contribute 
coherently to the scatter, while the contribution from 
adjacent slabs is uncorrelated statistically in phase 
or in amplitude. This intermediate value of d will 
be termed the scale of the irregularities of electron 
density and designated by the letter L. We note, 
for the case of a radar, that L is of order \/4. 

A much more precise handling of these concepts is 
embodied in the theoretical work done in recent years 
with specific application to the problem of radio 
wave scattering from irregularities arising in ion- 
ospheric or tropospheric turbulence [20, 21, 22]. 
Here the intensity of the scatter is found to be 
proportional to the integral over the scattering 
volume of the spatial autocorrelation function of 
refractive index, taking phase into account. For a 
more precise account of this development the reader 
is referred to the papers cited. What is important 
to our present problem is that this integral is formally 
equal to the spectrum of the irregularities evaluated 
at the wavenumber k for the case of direct backscat- 
tering as observed by radar. The reciprocal of this 
wavenumber, 1/k, is identified with the scale L 
multiplied by a constant close to unity. 

As noted by Balser [13], the spectrum of scatter 
from a large volume containing independent scat- 
tering centers is simply the sum of the spectra 
associated with the individual centers. In other 
words, in our approximation the spectrum of in- 
coherent scatter arising from a large volume of the 
ionosphere should have the form of the spectrum 
arising from any given slab of scale L (provided the 
mean characteristics of the volume are uniform 
throughout—temperature, density etc.). Since we 
have approximated the contents of such a slab as 
scattering coherently, the instantaneous field strength 
arising due to the slab is represented as proportional 
to the total electron content found at any instant 
within the confines of the slab. 


2. Role of the Ions in Determining the Criterion for 
Independence 


We now face the problem of deciding under what 
conditions the scattering is best considered to arise 
from independent particles, and under what condi- 
tions the representation using irregularities of particle 
density gives a more accurate result. We shall see 
that the ionosphere, being a partially ionized gas, 
provides a unique opportunity to answer this 
question. 

Gordon [1] made the assumption that the charac- 
teristic scale Z is meaningful only when the mean 
free path, 1, of the scattering particles is small com- 
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pared with L. He uses the converse of this assump- 
tion as the criterion for independent scattering from 
individual particles—i.e., when L<l. We shall 
refer to this as the ‘7 criterion.”” An alternate and 
very different assumption is that the criterion could 
be set by the relationship of the scale ZL to the mean 
separation of neighboring particles. This says that 
when L>n7™'*, one must use density fluctuations 
of scale L to calculate the scattering properties. 
Here n~* is the approximate mean separation be- 
tween neighboring scattering particles and n is the 
number of particles per cubic centimeter. We shall 
refer to this as the ‘“n~'“ criterion.”’ We note that, 
in the F region and above, / is of order 1 km or 
greater, while a typical F region electron density of 
10° per cm* corresponds to n~!“=0.01 em. The 
scale L, for the experiments reported in the first part 
of this paper, is of order 2 m. Thus we have 
n'8<L<l, ie., both criteria are satisfied. How- 
ever, since the independent particle approach sug- 
gests that a Doppler broadening much greater 
than actually observed should be found, we are led 
to inspect the consequences of the density fluctuation 
approach for the ionosphere. 

We now offer a crude argument which demon- 
strates that the distribution of the free electrons 
cannot be entirely random. We consider the case 
of scattering from a randomly distributed gas of 
neutral particles having approximately the same 
total particle density as that of the F' region of the 
ionosphere. Long-range forces among the particles 
are considered to be negligible so that only during 
collisions do the particles affect the motions of each 
other. Referring to figure 11 we consider the con- 
tents of a scattering slab AA-CC the thickness of 
which is equal to L. Since we have approximated 
that all particles located within this slab scatter 
coherently, the amplitude of the scattered contribu- 
tion from the slab is proportional to its total content 
of particles. Beginning at a time ¢ we note that at 
a time t4)-+7, the content of the slab has almost com- 
pletely changed due to the line-of-sight motions of 
the particles. Those particles which had been within 
the slab at t) have almost all moved to new locations 
and have been replaced within AA-CC by other par- 
ticles. The time 7 is therefore given approximately 
by the time taken for the average particle to move 
a line-of-sight distance equal to L, or 7=L/i,, where 
dv, is the root mean square line-of-sight velocity of 
the particles. Another way of saying this is that the 
content, or state of fluctuation, in the slab AA—CC 
loses correlation with its state at a previous time, if 
at least + seconds have elapsed since that previous 
time. Using the well-known relation between the 
autocorrelation function and the power spectrum [23], 
we therefore estimate that the bandwidth of the 
scattering from the slab of thickness Z is approxi- 
mately Af=1/r=2,/L. This corresponds to a band- 
width of Af=4/AX0.963,=0.960,/L from eq (1), not- 
ing that ,=1.0850 (see reference [12]), for the case 
of independent scattering from the particles. We 
see that the broadening expected for the two ap- 
proaches is the same. Turning to the actual case of 
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incoherent scatter from the F region of the iono- 
sphere, we believe that the velocities of the free 
electrons are equivalent to the supposed kinetic 
temperature of about 1,500 degrees Kelvin. Since 
the observed frequency broadening of the scatter is 
much smaller than expected for a fully random dis- 
tribution of electrons at this temperature, we look 
for an explanation in which some mutual correlation 
among the electron positions and velocities is found. 

We now explain the observed echo spectrum 
through the use of the density fluctuation approach. 
Insight regarding the distribution of the free elec- 
trons can be gained from a paper published by Pines 
and Bohm [24]. These authors considered the case 
of an electron gas embedded in a uniform smear of 
positive charge. Since the mobilities of electrons 
are so much greater than the mobilities of ions, this 
case is a good approximation to the situation in the 
ionosphere. Pines and Bohm found that the mean 
square intensity of fluctuations of electron density 
of scale size L is reduced, compared with that for the 
case of random motions, by a factor (27)?d?,/L?, for 
all {Z>Xp. In this relation Ap is the so-called 
Debye shielding distance, Ap=(kT/4ane?)” and T 
is the temperature in degrees Kelvin, k is Boltz- 
mann’s constant, n is the electron density and e is 
the charge on the electron. Only for the condition 
L<hp do the free electrons behave as though 
randomly distributed. Since Z is of the order of 
one or more meters for the present and other pro- 
posed radar experiments, and Xp is less than a few 
centimeters at all heights of interest, the independent 
random fluctuations in electron density of scale L 
must be viewed as negligible. The intensity of 
scatter arising from such fluctuations must also be 
negligible. 

Similar considerations hold for the case of a gas 
of positive ions embedded in a uniform smear of 
negative charge. In both cases the mutual repulsion 
of charges of like sign causes the reduction of fluctua- 
tions. However, Spitzer [25] has shown that in a 
real plasma, the mobility of the electrons is so great 
as to cause their equilibrium distribution to balance 
the fluctuations of positive charge density nearly 
completely. Once the fluctuations of positive charge 
density are balanced by equal and opposite fluctua- 
tions of negative charge density, the ions may be 
distributed nearly randomly, as they would be dis- 
tributed if they were neutral particles. The fluc- 
tuations of the free electrons from their equilibrium 
distribution are strongly reduced, in accordance with 
the argument of Pines and Bohm, inasmuch as the 
positive ion distribution departs only slightly from 
the approximation of a uniform smear of charge. 
Thus one is left with a picture of scattering from 
fluctuations of free electron density of scale LZ such 
that the magnitude of the scatter is the same as for 
a randomly distributed neutral gas, and its rate of 
change is appropriate to particles having the mean 
mass of the ions. The rate of fading must therefore 
be computed as if the positive ions, rather than the 
electrons, produced the scattered echo. This should 
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not be taken to imply that the ions make a direct 
contribution to the scatter, since the free electrons 
have a considerably greater radar cross section. 

For example, if we assume that the ionized con- 
stituent of the F layer is atomic oxygen, we estimate 
the Doppler broadening 


A c/s 


for an operating frequency of 41 Me/s and a temper- 
ature of 1,500 degrees Kelvin. 

The result of the arguments given in this section 
is that the characteristics of incoherent scattering at 
VHF, for n"°<L, if also }»><L, may be esti- 
mated by substituting for the ionosphere a fictitious 
gas of randomly distributed neutral particles with 
special characteristics. These particles would have 
the radar cross section of the free electron, and the 
mass (and consequent velocity distribution) of the 
ions. For n-“*L<Xp the Doppler broadening 
should approach that of independent electron scatter. 
For both cases, the scattering is best described in 
terms of fluctuations of particle density. Only when 
nL do we expect to find that the simple 
independent particle description of the scattering is 
correct. 


3. The Echo Spectrum 


In an informal report [16] and in the present 
paper, approximate calculations for the fading spec- 
trum are presented. Two more rigorous theories 
have been developed by other workers [11, 12] 
which, although using different approaches, predict 
the same spectral distribution. This theoretical 
spectrum has been experimentally confirmed by 
Pineo et al., at 440 Mc/s [17]. Two other theoret- 
ical papers dealing primarily with the intensity of 
the scatter have also appeared [13, 14]. 

At 41 Me/s, the bandwidth of the spectrum pre- 
dicted by these authors is about 1,000 c/s for 
echoes from the F region. Recent statistical studies, 
of the fading of echoes observed on the equipment 
described in this paper, show that the bandwidth 
of 41 Me/s echoes from the F region is certainly 
greater than about 600 c/s. This observation 
agrees with the results of Pineo et al. [17], and with 
the above mentioned theories. 

The author’s approximate calculation of spectral 
width differs from the precise theories by more than 
a factor of two. The difficulty apparently resulted 
from the assumption that echo contributions from 
adjacent volumes of scale Z can be considered to be 
statistically independent. It is true that the fluc- 
tuations in density of such adjacent volumes are 
independently distributed. However the electrons, 
found in a given volume at one instant, move into 
adjacent of antl while the phase of the individual 


contribution of each electron various continuously. 
This continuous variation of phase implies a sort of 

correlation among many irregularities of 
The assumption basic to the approximate 
None of the 


term 
seale L 
calculation was therefore not correct. 
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theoretical treatments mentioned above suffer from 
this difficulty. In spite of this difficulty with the 
author’s treatment, the heuristic arguments of the 
foregoing sections offer an approximately correct 
explanation of the physical principles involved. 
The rate of rotation of the resultant phasor repre- 
senting the contribution of the volume V must be 
related to the average of the rates of rotation of its 
component phasors. Since the individual compo- 
nents must have rates of rotation distributed in 
both advancing and receding senses relative to the 
reference, the average must be less than the mean 
absolute rate of phase rotation of the components. 


4. Magnetic Field Effects 


In the foregoing sections it has been assumed that 
the effects of the earth’s magnetic field can be ig- 
nored. In regions where the mean free path is 
large compared with the wavelength this can only 
be true for propagation parallel with the lines of 
force of the earth’s field, and for frequencies large 
compared with the gyrofrequency. The latter con- 
dition will always be satisfied in the earth’s iono- 
sphere and exosphere at practical operating frequen- 
cies. The quantity @ used in section 2 is related to 
the Maxwellian distribution of particle velocities 
and to the components of that distribution in the 
line of sight. If the magnetic field does not seriously 
alter the Maxwellian components of velocity parallel 
to the field, then the results of the preceding sec- 
tions will be approximately correct for a vertical 
incidence radar in temperate or arctic zones. 

If the propagation is nearly perpendicular to the 
lines of force of the earth’s field, magnetic control 
of the particles can be serious. The ions now 
execute paths with a circular component of motion 
the plane of which includes the line of sight. If Z 
is large compared with the radius of this circular 
motion, the particles are restricted to paths nearly 
parallel to the surfaces of constant phase. In this 
case the fading spectrum should be calculated as in 
the case ef short mean free path [11] using an equiv- 
alent coefficient of diffusion, D,,, at right angles to 
the magnetic field. When ZL is small compared with 
the gyroradius, there should be little effect on the 
main spectrum of fading. 

An interesting case arises when Z is of the same 
order of magnitude as the average gyroradius of the 
ions. All ions of the same mass must have the same 
gyrofrequency within a volume of scale L. The 
rate of migration of the center of rotation of any 
given particle out of the volume, Z thick and bounded 
by surfaces of constant phase, is given approximately 


by 


] 


where RF is the range to the volume in meters and the 
expression in brackets is the length of the first Fresnel 
Zone. The fading due to this component of motion 
would be only a few cycles per second at an operating 
frequency of 50 Me/s. The gyrofrequency of the 
largest ion expected, 03, is of the order of 10 ¢e/s at 
F region heights over the magnetic equator. The 
gyrofrequency for protons, likely to be encountered 
in the exosphere, is of the order of 250 c/s. Thus 
each ion may be expected to execute a minimum of 
several turns during the period taken for its center 
of rotation to move a small fraction of 360° phase 
from the radar station. We may therefore expect 
that the scatter echo will exhibit a tendency toward 
periodicity having a frequency equal to the ion 
gyrofrequency. This periodicity should be easily 
measurable using autocorrelation analysis of se- 
quences of pulse echoes. 


From the approximate frequency of the periodic 
component one may infer the mass of gyrating ions. 
If several ionic constituents are present, several 
periodic components should appear, the relative 
power in each being proportional to the number of 
ions of the appropriate constituent. In this way 
the scatter radar can be used as a sort of mass 
spectrometer to determine the variation of ionic 
constituents with height. Once the ionic constituents 
have been identified, precise determination of the 
gyrofrequency should provide a measure of the 
profile of the magnetic field intensity as a function of 
height. At the same time, the thermal broadening 
of the echo spectrum will be only slightly reduced 
by the magnetic field when ZL is of the order of 
magnitude of the gyroradius. Knowing the mass 
of the ionic constituents, it should be possible to 
make a precise determination of temperature. A 
more refined theory to take the magnetic effects into 
account is under development. 
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100 to 117 kMc/s Frequency Range’ 
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This paper describes a technique for the measurement of the absorption characteristics 
of millimeter radio waves over a wide pressure range and presents the results of measure- 
ments made with this technique for four gases in the frequency band from 100,000 to 117,000 
megacycles per second. The four gases investigated in near pure state were carbon monox- 


ide, nitrous oxide, sulfur dioxide, and nitrogen dioxide. 


follows: 
Predicted or Measured Meas- 
Gas previously New measure- absorption | ured line 
measured /[ment of frequency coefficient (a)| width 
frequency 
Me/s Me/s em-! Me/s per 
2 mm Hg 
Carbon monoxide-._- 8. 7510-5 4. 25 
Nitrous 1. 2510-4 4, 25 
Sulfur 104,046. 
105, 333. 94 | 9X<10-6 17.0 
Nitrogen dioxide_---- 105, 158. 4 No line detected 


1. Introduction 


The Electrical Engineering Research Laboratory 
of the University of Texas has made a number of 
measurements throughout the low-centimeter and 
millimeter wavelength regions of the absorption of 
radio waves due to atmospheric gases [1].2_ By re- 
peating these tests, the attenuation was plotted as a 
function of water vapor content. In almost every 
instance, it was possible to interpret the losses as 
being due to absorption by water vapor and oxygen. 

In one set of measurements, however, variable 
losses were noted which were not associated with 
changes in water vapor concentration. These losses 
occurred in the frequency range between 100,000 and 
118,000 Me/s and were attributed to some of the 
rarer gases in the atmosphere [2]. 

It is the purpose of this present paper to investi- 
gate the absorption properties of several gases with 
measured or predicted absorption lines in the fre- 
quency region in question. In order to isolate the 
effect of each line, the gases were studied in the nearly 
pure state. The absorption characteristics obtained 
by measuring the gases individually cannot immedi- 
ately be used to predict their absorption when di- 
luted with air. Further studies are in progress 
which, we bope, will assist in making this transition 
from individual gas loss to loss due to mixtures. 

1 Contribution from Electrical Engineering Research Laboratory, The Uni- 


versity of Texas, Austin, Tex. Prepared under Office of Naval Research Con- 


tract Nonr 375(01), N R 371032. 
2 Figures in brackets indicate the literature references at the end of this paper. 


The results of these tests are as 


2. Instrumentation and Measuring 
Technique 


The millimeter wave absorption was measured 
by the use of a circular tube 500 ft in length and 

in. in diameter. The cell may be evacuated and 
then filled with the gas under test. For the measure- 
ment of the larger absorption coefficients, a dia- 
phragm was inserted between sections of the pipe 
and only part of the length was filled with the test 
gas. The remainder of the tube was filled with 
nonabs»rbing gas. 

A Teflon diaphragm at either end permitted the 
injection and extraction of the electromagnetic 
energy. Horn antennas spaced back of the ends 
of the pipe provided proper illumination across 
the diaphragms. The electric field across the end 
of the cell was in a vertical direction with the 
intensity varying horizontally from maximum at 
the center line to zero at the edges. This field 
approximates that of the TE,, mode. Distortion 
from higher-order modes was very small since the 


attenuation in the guide for lossless gases was only 


db. 


The signal source consisted of a Bell Telephone 
Laboratories M-1805 klystron operating over a 
50,000- to 58,500-Mc/s range, a modified L-5506 
‘“in-guide” crystal harmonic generator to provide 
the second harmonic for the tests, and the associated 
power supplies and necessary waveguide components. 
The maximum signal available was of the order of 
10-7 w. 

The receiver also utilizes an M-1805 klystron, 
an “in-guide” wafer crystal for harmonic generation 
and mixing, and the necessary amplifiers and power 
supplies. The receiver was frequency modulated 
to reduce the d-c amplifier requirements and to 
compensate for klystron frequency instability. 

The usual measurement procedure was to evacuate 
the cell and determine the attenuation at the fre- 
quencies of interest as the test gas pressure was 
increased. A continuous frequency scanning tech- 
nique was used at low pressures in the vicinity of 
the maximum absorption. 


| 
| | 


3. Carbon Monoxide 


As is the case of each of the various gases studied, 
the frequencies of the absorption lines of carbon 
monoxide had been predicted from a study of the 
infrared molecular energy levels. The various 
resonant frequencies have been tabulated by Ghosh 
and Edwards [3]. Gordy and his associates have 
measured with great accuracy the frequency of the 
carbon monoxide line in the spectrum of interest 
in this paper. A comparison of these values is as 
follows: 


115,270.2 Me/s 
115,270.56 Me/s. 


Predicted: 
Measured (Gordy): 


Equipment was not available at this laboratory 
for absolute frequency measurements to this preci- 
sion. The above absorption line at the upper end 
of the band and one for nitrous oxide at the lower 
end of the band, which was also measured with 
accuracy by Gordy, were used as calibration points 
on the wavemeter. By assuming linear variation 
between these points, a measure was obtained of 
the intermediate frequencies. 

The absorption versus frequency measurements 
in the vicinity of the carbon monoxide line are 
shown in figure 1. Pressures ranging from 4.5 to 
27 mm of mercury were used in obtaining these data. 
The lines drawn on figure 1 were calculated by the 
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Figure 1. Carbon monoxide absorption spectra. 


Absorption coefficient as used in the figures is in nepers per centimeter. 
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use of the Van Vleck-Weisskopf equation [5] with 
the maximum attenuation in the equation being 
adjusted to fit the maximum measured value at 
the pressure in question. A line breadth propor- 
tional to pressure with a value of 4.25 Mc/s per mm 
Hg was obtained from the half power points at 
very low pressure where continuous sweep of the 
frequency was used. This line breadth constant 
provided a fairly good fit of the measured data 
when the maximum absorption was adjusted to 
fit the measured maximum at each pressure. This 
variation of maximum absorption with pressure is 
not implied in the Van Vleck-Weisskopf equation, 
but it may possibly be explained by saturation 
effects [6]. This decrease of absorption at the 
lower pressures has been attributed to the power 
saturation of the molecules and the possibility that 
lines could not be completely resolved due to a 
small amount of frequency modulation present in 
the measuring signal. Power variations over a 
range of several decibels did not, however, seem 
to alter materially the shape of the maximum 
absorption-pressure curve. 

To further investigate the pressure dependence 
of the absorption, this peak attenuation was 
measured over a pressure range from a few milli- 
meters to nearly 1 atm. These data are shown in 
figure 2. The absorption is seen to increase rapidly 
at very low pressures and slowly at higher pressures. 
The higher pressure increase is thought to be due 
to the “skirts” of higher frequency lines. The 
absorption coefficient, a, is chosen as the value at 
which the curvature changes abruptly since the 
gradual increase with pressure above the knee of 
the curve is attributed to contributions from other 
lines. The saturation effect is assumed to be 
confined to the very low pressures. 


1.5 x104 1 1 1 130.29 
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Ficure 2. Absorption versus pressure for carbon monozide. 


Measured values. 
4. Nitrous Oxide 


Gordy and his associates [4] had measured with 
precision the frequency of the nitrous oxide line 
predicted for the spectrum covered by this report. 
A comparison of the measured and predicted values 
is as follows: 


Predicted: 
Measured (Gordy): 


100,493.0 Me/s 
190,491.76 Me/s 


| 
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As indicated earlier, the frequency of this line was 
used as a reference in the calibration of the wave- 
meter for the frequency range in question. 

The nitrous oxide absorption spectra for a number 
of pressures are shown in figure 3. The curves in 
the figure were drawn using the Van Vleck-Weisskopf 
equation with a line breadth proportional to pressure 
for a proportionality constant of 4.25 Me/s per mm 
Hg. This line breadth constant was also determined 
from the half-power width at very low pressures. 
As in the case of carbon monoxide, it was necessary 
to adjust the peak absorption for the particular 
pressure used. 

The variation of the peak absorption as a function 
of pressure for nitrous oxide is shown in figure 4. 
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Fiaure 3. Nitrous oxide absorption spectra. 
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Ficure 4. Absorption versus pressure for nitrous oxide. 
Measured values, 
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The low-pressure change may be explained in part 
by saturation effects and the high-pressure change 
by the effect of the skirts of other lines. The absorp- 
tion coefficient applicable to the line is chosen as the 
value at the point of most rapid change in slope as 
in the case of carbon monoxide. 


5. Sulfur Dioxide 


Ghosh and Edwards [3] list four predicted sulfur 
dioxide lines in the vicinity of the spectrum covered 
by this report as follows: 


100,911.75 
104,028.28 
105,333.94 
116,790.07 


ep 40 additional lines between 44,000 and 200,000 
Me/s. 

The sulfur dioxide spectra were examined in detail 
over the range from 102,000 to 106,000 Me/s and 
the results are shown in figure 5. The effects of the 
first three lines listed above are readily seen in this 
figure. In addition, the complex spectra of sulfur 
dioxide was evident in the flattening of the absorption 
spectra over the region investigated. 
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Figure 5. Absorption coefficient of sulfur dioxide versus 
frequency at different pressures. 


The line in the vicinity of 104,000 Mc/s was 
examined at low pressures by frequency scanning 
and was found to consist of two lines with a separa- 
tion of 4.5 Me/s centered approximately at 104,042 
and 104,046 Mc/s. Each line had a width of 13.7 
Mc/s per mm Hg and the absorption of lower fre- 
quency component was approximately 1.4 times 
more intense than the higher frequency one. 

The maximum absorption of the combined line 
as it blended together at increased pressures was 
measured and is shown in figure 6. The lower line 
in this figure is the absorption which would be associ- 
ated with the double line in the vicinity of 104,000 
Me/s. The frequency for which the calculations 
were made falls between the two lines and therefore 
goes essentially to zero at very low pressures. The 
difference in the measured and calculated values is 
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Ficure 6. Maximum absorption coefficient for two observed 


lines of sulfur dioxide versus pressure. 


apparently due to the attenuation caused by the 
skirts of the many other lines. The absorption of 
the 105,189 line was so low that it could not be 
resolved at the lower pressures. Its effect is evident, 
however, at higher pressures (fig. 5). 


6. Nitrogen Dioxide 


An absorption line had been predicted at a fre- 
quency of 105,158.4 Me/s for nitrogen dioxide. A 
continuous search was made from 104,200 to 105,600 
Mc/s at very low pressures and spot measurements 
were made at 104,720, 104,920, and 105,120 Me/s. 
In none of the tests was there any indication of 
appreciable attenuation. It was felt that if a line of 
nitrogen dioxide exists with significant strength of 
absorption, its effect would have been noted. 


7. Conclusions 


The measurements described in this paper have 
confirmed the presence of absorption lines for carbon 
monoxide, nitrous oxide, and sulfur dioxide in the 
frequency range from 100,000 to 117,000 Me/s. A 
predicted line for nitrogen dioxide could not be 
detected. The data can be represented by the Van 
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Vleck-Weisskopf [5] line shape function using a line 
breadth proportional to pressure and a maximum 
loss at a given pressure adjusted to fit the measured 
maximum. One of the predicted lines for sulfur 
dioxide was found to be a doublet. The numerical 
results are summarized in the following table: 


Predicted or | New measure- Measured | Measured 
Gas previously ment of absorption line 
measured frequency coefficient width 
frequency (a) 
Me/s Me/s Me/s per 
2 mm Hg 
Carbon monoxide 8. 7510-5 4.25 
13.75 
Sulfur 104,046. | | 
105, 333.94 | 105,189_____....-- | 9 X10-6 17.0 
Nitrogen dioxide___-. 105, 158. 4 No line detected 
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(May 31, 1960) 


Formulas are developed for the diffracted field around a perfectly conducting cylinder 
with a dielectric sleeve of arbitrary thickness. These formulas represent the field due to a 
unit electric line source parallel to the cylinder (either inside or outside the dielectric sleeve) 
as a spectrum of radial eigenfunctions. It is shown that in each case the field in the region 
containing the source can be expressed as the sum of two terms, the first of which is a spectral 
representation of the field when the outer as well as the inner surface of the dielectric is a 
perfectly conducting boundary. The second term of the sum, which alone involves the 
properties of the dielectric, is an integral that converges rapidly at high frequencies. Using 
these general results, perturbation calculations are carried out for three limiting examples 
of plane wave scattering: (1) Thin dielectric; (2) low-refractive index; (3) very small sur- 
Pe curvature. In the latter approximation the correspondence with optical results is 
shown. 


1. Introduction 


The scattering of plane electromagnetic waves by cylindrical obstacles has received detailed 
theoretical treatment only in the two special cases for which the scalar wave function, ¥, 
satisfies either the Dirichlet or the Neumann boundary condition at the cylindrical surface. 
In electromagnetic terms these problems correspond respectively to the scattering by a per- 
fectly conducting cylinder of a wave whose electric or magnetic vector is parallel to the cylinder 
axis. For a cylindrical obstacle whose surface is not perfectly conducting a more complicated 
boundary condition of the type ¥+ad0,¥=0 is required (sometimes called an impedance 
boundary condition). In such cases the theory has not been so extensively elaborated.* 

This paper is concerned with some of the theoretical features of the scattering problem 
associated with a cylindrical obstacle of this more general kind. In particular we shall examine 
the scattered wave produced by a dielectric cylinder of circular cross section having a perfectly 
conducting coaxial core. One aspect of this problem has been investigated experimentally by 
Tang [1] * for an incident wave with electric vector polarized parallel to the cylinder, and in 
the same paper the Fourier expansion of the wave function has been used for numerical calcu- 
lations. It seems, however, that no general discussion of the characteristics of the wave 
solution has been given; it is this task which will now occupy our attention. 


2. Derivation of the Green’s Function 


The cylinder which we shall consider (fig. 1) is composed of a perfectly conducting core 
of radius a surrounded by a dielectric layer of inner radius a and outer radius 6. In the sub- 
sequent analysis the relative dielectric constant of the layer will be taken to be a real scalar 
greater than one, ¢-=e€/@>1. That is to say, the refractive index of the dielectric is a real 
number greater than that which characterizes the surrounding medium (air or vacuum). For 
a lossy or anisotropic dielectric the theory is unaffected except by the substitution throughout 
of a complex or tensor dielectric constant for «. Of course, the detailed calculation of field 
quantities will be much more complicated in such cases. 

If the composite cylinder described above is excited by a harmonic line current parallel 
to its axis (electric line source, exp[—iwt]), the task of finding the resultant electric field at 


1 Contribution from the Division of Engineering, Brown University, Providence, R.I. 
2 The research reported in this paper has been supported in part by the Office of Naval Research and the David W. Taylor Model Basin 
under contract Nonr 562(24) and in part by the Electronics Research Directorate of the Air Force Cambridge Research Center under contract 


AF 19-(604) 4561. 
3 Certain aspects of such problems have been considered by J. R. Wait, Electromagnetic radiation from cylindrical structures, Chapter 


16 (Pergamon Press, London, 1959). 
4 Figures in brackets indicate the literature references at the end of this paper. 
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FicurE 1. The composite cyl- 
inder and coordinate system. 


every point outside the perfect conductor resolves itself into the determination of the Green’s 
functions which satisfy appropriate partial differential equations and boundary conditions. 
For the moment we shall limit our considerations to the case in which the line source is outside 
the dielectric. Denoting coplanar position vectors relative to an origin on the cylinder axis 
by r’ and r for source and observation points, respectively, we require solutions in the zy-plane 
for the scalar wave equations 


a<r<b,r’>b (2.1b) 


subject to the following boundary conditions at r=a and r=b: 


G,(a,r’)=0; (2.2a) 
G)(b,r’) (b,r’); (2.2b) 
0,Go(b,r’) =0,G,(b,r’). (2.2c) 


Here the wave number, k=w/v, is related to the dielectric constant and the refractive index, 
n, of the dielectric according to 


k=kyvVe,=nkp. (2.3) 


The mathematical problem defined by eqs (2.1) and (2.2) may be solved in a variety of 
ways. We shall employ a technique developed and discussed by Titchmarsh [2] and Marcu- 
vitz [3] for separable equations. The procedure consists in constructing so-called characteristic 
functions which satisfy the related system of ordinary differential equations 


E ¢ (r,r’ +) (2.4a) 
E |p =0, (2.4b) 
g (0,0’ 4) =—6 (06—6’). (2.40) 


These equations are found from (2.1) by separating the polar variables (r, #) with \ taking the 
role of separation constant. In addition to (2.4) the functions f, and f, are required to satisfy 
boundary conditions corresponding to (2.2), while the angular function g is made unique by 
certain convergence conditions to be specified. Once the characteristic functions have been 
explicitly determined from (2.4) and the boundary conditions, the Green’s functions required 
by (2.1) are obtained by an integration in the complex \-plane. For it can be verified that if 
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C is a circular contour of infinite radius in this plane, then 


1 , , 
ay (r,r’ ,A)=ré(r—r’), (2.5a) 
1 
arg (0,0’r) =6 _(2.5b) 


With the help of (2.4) and (2.5) direct substitution in (2.1) will show that the required Green’s 
function is 


G(r’) =5— | 9 (0,8" (2.6) 


where T is a contour which encloses the singularities of one of the functions f or g, the sense of 
the contour being given either by (2.5a) or (2.5b) as the case may be. 

In carrying out the procedure outlined above, we shall begin by considering the charac- 
teristic function which satisfies (2.4c). It has the form 


9 (6,0" A) = (2.7) 


where v?=) and @s, 6< are respectively the larger and the smaller of the two variables 6 and 6’. 
The positive exponential solution has been chosen to ensure the regularity of g for all values of 
d on the principal sheet of the two-sheeted \ surface. On this sheet Im »>0, and since 


0< arg v= arg <n, (2.8) 


the region of regularity corresponds to the upper half of the v-plane. The constant A in (2.7) 
is determined by the fact that dg/d@ has unit discontinuity at 6=6’ for a delta function source. 
We thus find from (2.4c) that the characteristic function is 


+. ‘ 
90,0" el ol O< arg (2.9) 


It is easy to verify that (2.9) satisfies (2.5b). Since the only singularity of g(@, 6’, \) is a 
branch point at \=0, the clockwise contour C may be deformed to an integral around the branch 
cut which is conveniently taken along the positive real axis of the \-plane. If the deformed 
contour is now mapped onto the real axis of the v-plane by the transformation \=»’, the 
resulting integral is seen to be a representation of the 6-function. Thus 


The characteristic Green’s functions which satisfy (2.4a) and (2.4b) are constructed in a 
similar fashion but from linear combinations of Bessel functions of order y= X . At large dis- 
tances from the source, r>r’, they must represent outgoing waves. Furthermore, the interior 
function must vanish at the conducting surface r=a, and the exterior function must behave 
like a “standing wave” in the region 7r’>r>6. In terms of Hankel functions of the first and 
second kinds the most general expressions which exhibit these characteristics are 


= Bol Hy” (kor<) (kob) (kor) (kor>); (2.10) 
AQyr’ (ka) H? (kr) —H? (ka) HY (kr) H (hors); aS <b, >b. (2.11) 


Two of the three integration constants appearing in (2.10) and (2.11) are determined by the 
pair of linear algebraic equations which express the continuity of the radial function and its 


first derivative atr=6. The remaining constant is fixed by the jump in ») at 
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Making use of the Wronskian relation 


(BH (8) (8) (2.12) 
U 
and introducing the*notation 
a=ka, B=khkb, Bo=kyb; (2.13) 


°C, (ker) =H,” (a) H,” (kr) (kr), 


Ci (kr) =H” (ker) —H,? Hy? (kr); (2.14) 
and 


we find that the characteristic function for an exterior source at 7’ is 


H;” 


fo(r,r’ [ er me r>b; (2.16) 
, CO, k ky 


A direct calculation shows that. these functions satisfy (2.5a). 
It remains to determine G(r,r’) from (2.6). First we note that the only singularities of 
f, and f; in the d-plane are simple poles at the set of complex values, \,, for which 


=0. (2.18) 


That these functions have no branch point is easy to verify. With the relations, 
and 


we learn that 
C_,=C, and X_,=X,. 


It follows that f) and f; are both single-valued over the entire \-plane. In accordance with 
(2.6), we now form the integral 


i rs : 


where the contour T is chosen to enclose the singularities of f in a positive sense. Since the 
product function has a branch point at A=0, I must be indented around the branch cut, which 
we shall take along the real axis as shown in figure 2. In the v-plane the resulting contour 
maps toy. By virtue of (2.8) and the asymptotic behavior of the Hankel functions, the contri- 
bution to (2.19) from the semicircular segment of y may be made arbitrarily small if the radius 
of the arc is made large enough. In this way we see that both a discrete and a continuous 
representation of G@ are possible, 


5 It would be more explicit to write C,(a, kr), but since subsequent operations on C, affect it only as a function of kr, it is perhaps clearer to 
use the abbreviated notation. 
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Ficure 2. The complex v and X planes. 


d- plane 


3. The Discrete Spectrum 


The discrete representation of the Green’s function is obtained by evaluating the residues 
of f(r,r’,A) from (2.16) and (2.17). A more symmetrical form results if one makes use of the 
relation 


1 


which, at v=v,;, becomes simply 


4 1 
(2) 
HP Bo) (3.2) 


where the subscript indicates the 7 root of (2.18). With this change the residue sums corre- 
sponding to (2.19) form discrete spectral representations of the total field due to an exterior line 
source aur’ >b. They are 


eer) HP (he) 


It is now clear from considerations of symmetry that if the line source is inside the dielectric, 
a<r’ <b, then 


It is noteworthy that the representations (3.3a) and (3.4b), with source and observation points 
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in the same region, are unaffected by an interchange of the variables r and r’. If r and 7’ are 
in different regions, however, their interchange turns (3.3b) into (3.4a) and vice versa. 

The radial eigenfunction expansions derived above have certain characteristic features 
which should be mentioned. For fixed r and r’ each term of the series represents (in the 
terminology of Franz) a “creeping wave’ propagating in the 6-direction with exponentially 
damped amplitude. Since @ covers the interval (—, ©), the total field at any point is the 
sum of the contributions at 6+2mz for all integral values of m. At high frequencies, however, 
the imaginary part of v is a large positive number so that all terms other than the one for which 
m=0 are exponentially small. Asymptotically, (3.3) and (3.4) represent the total field in the 
range 0<|6—@’|<2r. It is to be expected, moreover, that the set of v; are ordered with in- 
creasingly large imaginary parts so that if |@—6’| 40, only the first term of each residue expan- 
sion will contribute significantly to the field at high enough frequencies. However, the problem 
of finding even the first root of (2.18) would seem to be a formidable one. 


4. The Continuous Spectrum 


An alternative representation of the Green’s function as an integral over a continuous 
spectrum is obtained from (2.19), using (2.16) and (2.17). Since g(6, 6’, \) has a branch point 
singularity at A=0, the contour integral becomes 


where the indicated path is around a branch cut along the real axis of the d-plane. With the 
help of (3.1) and the change of variable, \=v’, we find from (4.1) that 


dv Hy (koro) Hy (kors) (4.2) 


Similarly the field inside the dielectric is 


(1) , 


In this form the Green’s function of (4.2) has some interesting features. The first term 
satisfies all the conditions associated with the field of an electric line source in the vicinity of 
a perfectly conducting cylinder of radius b when the medium is characterized by the propaga- 
tion number kp. One way to see this is to close the integration contour in the upper half of the 
v-plane. The residue sum that results is 


where 
(By) =0. (4.5) 


From the Wronskian, (2.12), one finds that 


HP 
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so that (4.4) can be rearranged in the form 


H® (kor) H® (kor’) civ 


G(r, H® (Bo) © (4.7) 
H® (By) H®’(B,) (4.8) 


These are the results for a conducting cylinder which have been derived in other ways by T. 
Wu [4] and by the author [5]. 

We see then that the field of a composite cylinder has as one part the creeping waves of 
a conducting cylinder. Only the second integral of (4.2) depends upon the parameters of the 
dielectric (thickness and refraction index) and this dependence resides in the factor X!”. 
Indeed as a6, X{” grows without limit and the integral vanishes, leaving the Green’s function 
for a conducting cylinder of radius 6. On the other hand, as k—>k, the limiting value of the 
second integral of (4.2) is 


H® H® 


and on substituting we find that the Green’s function reduces to that for a conducting cylinder 
of radius a. These limiting forms will be used later for perturbation calculations. 

When the line source is inside the dielectric sleeve, the derivation of the Green’s function 
proceeds in a similar way. We begin with characteristic functions for the radial equations 
having the form 


ho(rr’ A) (kr<) (hers); aXr’<b, (4.9) 
hi(ryr’ )=AiC, (kr) (krs) + BH® (krs)]; a<r,r’ <b. (4.10) 


Both functions satisfy the boundary condition at rc=a, and ho represents an outgoing cylin- 
drical wave as rs. The constants can be evaluated as before, and after a certain amount 
of manipulation, we find that 


N= C,(8)H® (Bo) (4.11) 

A)= 4 (8) (krs) Hy? (kr>) C,(B)C(B) (4.12) 


The exterior function ho(r,r’,) of (4.11) corresponds exactly to f, of (2.17) after an interchange 
of the variables r and r’, so the related Green’s function corresponds to (4.3). When both 
source and observation points are inside the dielectric, the convolution of h,(r,r’,A) with (2.9) 
around the branch cut of g(0,0’,A) leads to the following integral representation of the Green’s 
function: 
Add 4.13 


The first term of this Green’s function vanishes at re=a and r,=6, and one can easily verify 
that it corresponds exactly to the field of an electric line source inside a dielectric filled coaxial 
transmission line with perfectly conducting walls. The second term represents the modification 
of that field due to radiation into the exterior medium. 
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The residue sum associated with (4.11) has already been given in (3.4). To complete the 
discussion we shall write down the symmetrical representations for a coaxial line, analogous to 
(4.7) and (4.8). They follow from the convolution of the first term of (4.12) with g(@,6’,\) 
around the simple poles located at the points of the \-plane for which 


C,, (8)=C, (6) =0. 


It is found that 
Cs 


14 
C, (ker) C, (kr’) 41 
where we have made use of the Wronskian relations, 
(kr) (ker) (kr) C, (kr) HS (@). (4.16) 


5. Plane-Wave Excitation 


In order to render more manageable the integrals appearing in (4.2) it is convenient at 
this stage to turn one’s attention to special cases. Since the general class of problems asso- 
ciated with the scattering of a plane wave is of considerable interest, we shall first consider the 
question of reducing the field to that of an incident plane wave of unit amplitude. This may 
be done by considering the first term of (4.2), which we shall call G,, 


The Hankel functions, H,(p), grow as exp(» In =") when y+, so the individual terms of this 


integral diverge and cannot be treated separately. However the divergences may be removed 
by first putting (5.1) in the form, 


and then making the substitution, 
(p)=2J,() — Hy” (p). 


With these changes it is found that 


(5.3) 
Each of these integrals is absolutely convergent. 

To establish a correspondence between (5.3) and the required plane wave solution, we 
follow Wu’s procedure and sum the creeping waves for the first term. Then by interchanging 
the order of summation and integration (or by invoking the Poisson summation formula) it 
can be shown that 


re 


& 
3 
| 
| 


where the function on the right is just the free space Green’s function representing a unit line 
source. If we now let rs and make use of the asymptotic form, 


2 


both sides of (5.4) can be normalized to an incident plane wave of unit amplitude by multi- 
plying with the factor —42(2)~"?(kor>)exp —t(kors— 7/4). Making these changes in (5.3) and 
setting 6’=0 gives the total field, Y,(r), when a plane wave from the right is scattered by a con- 
ducting cylinder of radius 6. This field is 


0 0 (Bo) 


The procedure that led to (5.4) will also show that 


When (4.2) has been normalized in accordance with (5.5), the total field outside a dielectric 
clad cylinder illuminated by a plane wave from the right becomes 


6.7) 


The last term is the modification of the scattered field introduced by the dielectric. It can also 
be put into the form, 


d H (kor) Sine 
yo |? cos v6, (5.8) 


which converges for all values of 6. 


6. The Surface Distribution 


In this section the distribution of y on the dielectric surface will be examined in a number of 
special cases. This choice of field point does not affect the methods of approximation that will 
be used, but does simplify the general expression for y to a certain extent. For on the circle 
r=b the first term of (5.7) vanishes, so that the integral to be considered is simply 


2 (0-5 


In terms of the index of refraction, n=k/ko, and the quantity, 6=$—a, the conditions under 
which approximations to (6.1) will be investigated are (a) 6/81; (b) |n—1\/<1; and 
(ec) 


6.1. Thin Dielectric: 6/8<1 


Since all the dielectric effects are contained in the function, 


Cy (a, B) H,”"(Bo) 
Cy(a, HY (By)’ 


(6.2) 


we require the expansion of 1/X, when a~b. It is clearly advantageous to treat X, as a function 
of the variable a@ (i. e., the inner radius @) since that parameter appears only in C,. Moreover, 
H» (Bo) has no zeros for real v, so the first term of (6.2) predominates for small 6. We may 
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therefore put 1/X,~F,, where 


Since lim F,(a)=0, the Taylor expansion of (6.3) has the form 


F,(6—8) (8) +5 (6.4) 


With (6.3) and (6.4) it is a straightforward matter to develop the required series explicitly. 
We need certain relations of the Wronskian type [6], and up to terms of order 6‘ the result of 
the calculation turns out to be® 


2B 3 B 4 3 Vv Vv 
It is clear that the first group of terms represents the function —In[1—(6/8)], so that 
] b 186 
245% +0 (5) (6.6) 


On substituting (6.6) in (6.1) a series of integrals is obtained, the first of which is 


2° 
(no) 


The calculation of 0,y,(b) from (5.5) shows that except for the constant, In (6/a), the last result 
is just the current distribution on the surface of a perfectly conducting cylinder of radius 6. 
Therefore, to a first approximation for thin dielectric sleeves of any refractive index and 
radius, the electric field distribution on the dielectric surface is proportional to the magnetic 
field distribution on a conducting cylinder of the same size when illuminated by a plane wave. 
The first term involving the refractive index is 0(6/8)* and is proportional to n?. 


6.2. Low-Refractive Index: |n—-1|<1 


In order to expand the surface distribution, 


dy 


in the neighborhood of ko, let us consider the function, 

G(k) (Bo) (6.7) 
Its Taylor expansion may be written in the form 

=G, (ko) + (n—1) (ko) + O(n—1)?, (6.8) 
where (n-1)ky=k—kp. Using (6.2) and the Wronskian (2.12), we find at once that 


C,(a, H (a) 
=— 4 ( H? (Bo) A® (ae (bo) (6.9) 
6 The successive derivatives of Fy(a), evaluated at a=, are 
Fy" (@)= (1+6?—»2) 
B 
i B 
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Thus the first term in the expansion of (6.1) is just ¥,(b), the total field at r=b when the plane 
wave is scattered by a conducting cylinder of radius a in the medium kp. 
The second term is found by evaluating 


(Ko) = (Bo) G5 (Ko) (6.10) 


where 0, signifies the differential operator 0/ok. Treating C, as a function of two dependent 
variables a, 8, and replacing C) by OsC,, we have 


OF (a, B) 


a,.X=a, E 


It follows, with the help of (6.9), that 


2 (1) 


Since 0,=a0a+5dg,, the final bracket of (6.11) becomes, after a little rearrangement, 
I=kab(C,0a%C,— — Bb(OsC,)? + bC,(OsC,+ BO3C,). (6.12) 


If we now make use of Bessel’s differential equation for C,, the last parenthesis of (6.12) 
is seen to be 


+ ?—B)C.. (6.13) 


Moreover, the bracketed expression in the first term of (6.12) can be represented by the 
Wronskian determinant 

C, OC, | 


(6.14) 


In the usual way [7] one can show that 


which implies that BW is independent of 8. To find the value of the constant we let B=a. 
Then C,=0, and 

lim 0. 

Boa 


In the limit the derivatives are ordinary Wronskians like (2.12), and it follows that 


16 
(6.15) 


The result of substituting (6.13) and (6.15) into (6.12) is 


161 


With this result, eq (6.8), (6.9), and (6.11) can be combined into an explicit expression for the 
first term of the expansion of G,(k) near k=ky. Substitution in (6.1) yields the following 
perturbed distribution at r=) when n~1: 


+3) | a6) 
29 


& 


If one is interested in the perturbation introduced by the dielectric at points outside the 
sleeve, it is only necessary to put H{” (kor) for H{ (89) wherever that function appears explicitly 
in (6.16). Whether the field point is on the surface or outside the sleeve, however, it is readily 
verifiable that all the integrals converge and that the bracket vanishes when b=a, leaving 
only the field of the conducting cylinder. In general, it seems that none of the correction 
terms in (6.16) can be neglected. 


6.3. Very-High Frequencies: k,a>1 


The last method of approximation to be dealt with is useful when both inner and outer 
cylinders are large compared to the wavelength in either medium. By combining the longer 
wavelength with the shorter radius these conditions are satisfied by the single relation, a>1. 
When this requirement is met, no further restrictions need be imposed, a priori, on either the 
thickness of the dielectric sleeve or on its refractive index. Since we shall consider only the 
contribution to (6.1) from that part of the integral for which |»|/+-3|»|!/°<ap, it is possible to 
use the Debye approximation for all Hankel functions. To this order of accuracy we may 
expect that the integration will yield the optical rays, but no diffraction effects. 

For a typical Hankel function the appropriate asymptotic form is 


H® sin mo) exp iBof (no), (6.17) 
where mp is the real angle 
no= cos” *(»/Bo), (6.18) 
and 
(m0) =sin COS 19. (6.19) 


The corresponding approximation for the derivative with respect to the argument, §p, is 
(Bo) sin noH;” (Bo). (6.20) 


With (6.17) and (6.20) it follows that 
(Bo) 


Bo ~ sin No- (6.21) 
In a similar way, by defining the angles 
n=cos~'(v/8); §=cos™'(v/a), (6.22) 
and the function 
F(v)=Bf (n)—af(§), (6.23) 
one can show that 
B sin 7 (6.24) 
If (6.21) and (6.24) are now combined according to (6.2), we find that 
1 1 ip 1 — e2iF (v) 
where p and 7 signify the ratios, 


Now for all 63, sin sin m,s0 thatO0>p>1. Consequently, 0<7<1 and |r exp 2iF(v)|<1. 
The denominator of (6.25) can therefore be represented by its Taylor expansion, and it follows 


that 
1 1 a 
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7) 


This result already makes evident the ray nature of the approximation. Successive 
terms of the series correspond to the rays which are, (1) reflected by the sleeve, (2) reflected 
once at the conducting cylinder, (3) reflected once at the conducting cylinder and once inter- 
nally at the dielectric surface, etc. The phase function of each term corresponds to the number 
of path segments a ray traverses in the dielectric (with a minus sign added for each reflection 
at the conducting cylinder), and its amplitude function includes the appropriate transmission 
and reflection coefficients for each contact with the dielectric surface. Since these quantities 
are all functions of the integration variable, successive terms can be evaluated in detail only 
by more and more intricate stationary phase approximations. There is, however, one limiting 
case for which it is relatively easy to establish the correspondence between the terms of (6.27) 
and optical rays. When the radii a and b become indefinitely large, while their difference, 
d=b—a, remains constant, the solution for the cylindrical problem at each point b= (6,6) of the 
dielectric surface approaches that for a plane wave scattered by a conducting plane sheet 
with a dielectric coating of thickness d, the dielectric surface being tangent to the cylinder at 
the point in question. 

To see how this comes about let us substitute (6.27) and (6.17) in (6.1). The first integral 
of the resulting asymptotic series becomes 


Q\'2 «(4 p ‘ 
»b)=(=) dy(By Sin ne) i+, (v), (6.28) 


where the limits on the integral define the interval within which the Debye approximations 
are valid and where 


=—Bof (n) vr, 0<86. (6.29) 


For large 8) the principal contributions to this integral come from the neighborhood of the 
stationary points of g“(v) in (—A, A). Since the asymptotic integrand decays exponentially 
outside this interval, one may extend its limits to (— ©, ©) without appreciable error. Then 
if a stationary point of g"(v) is located at v=», the Taylor expansion of (6.28) in the vicinity 
of this point gives as its first approximation, 


(1) -1/2 1 i- (1) 
~( Gin drexpi[ a9) | (6.30) 
where A=v—», and the subscripts indicate functions that are to be evaluated at v=. 

The value of »v required in (6.30) is that which reduces 0,g°(v1)=0 to an identity. It 
can be shown that this occurs for the angle ():=37—8, corresponding to 


sin 0 (6.31) 


According to the condition specified prior to (6.17), the angle 6 must therefore be such that 
sin 6-+38-2 (sin (6.32) 


It is clear that this inequality cannot under any circumstances be satisfied at the shadow 
boundary where 6=72/2. With (6.31) and (6.26) we determine that 
d2g (»,) 
By cos 
and 
2p1 2B cos 6 
1+p: cos 0+8 cos 0, 


where RF is the reflection coefficient for a plane dielectric surface [8] and the angles of incidence 


1+R, (6.33) 
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and refraction are related by 
sin @ _ 


=n. 
sin 6, 


With these results, we have 


(b)~3 Je By cos 6 cos 6 34) 


The value of the infinite integral is given in standard tables, and (6.34) can be reduced to 


(b) ~ (1+ 0088, (6.35) 


which corresponds in amplitude and phase to the sum of the incident and reflected rays at the 
dielectric interface for angles 0<t (fig. 3). Note that further restrictions on 6, such as (6.32), 


are no longer required. 


Ficure 3. Depicting ray geometry. 


The second term in the expansion of (6.1) is 


(b) --(?) e4 dv (By Sin o)~"/? (1+ p)? exp (»), (6.36) 


where 
(vy) =2F (v) +g” (). (6.37) 


Again, after finding the point v, at which 0,g(v.) vanishes, we can approximate (6.36) with 


where \=v—y, and the subscripts indicate functions that are to be evaluated at v=y». 
According to (6.37) the required stationary point is the solution of 


0,[2F'(v) (v)]=0. (6.39) 
With the help of (6.23) and (6.29) it may be shown that this equation is equivalent to 


1 
—2(n—) + (6.40) 


In order to discuss the solution of (6.40) in the plane limit we shall assume that (7—¢) is small so 
that sin (y—f) is approximately equal to its argument. Then 


(6.41) 


sin 7 B 
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where 6=kd and terms of order (6/8)? have been neglected. Accordingly, if (y—¢) is to be a 
small quantity it is sufficient that 6/8=d/b< <1, provided 70 or what amounts to the same 
thing, v?~6?. This requirement is, of course, the same as that already imposed by the original 
Debye approximations. With this understanding, (6.40) can be put into the more convenient 
form 


cos [ =sin 6) (6.42) 


Expanding the left-hand side of this relation up to terms of order 6/8, we obtain 


where as before 


Bo Cos 


sin 8; pi= 


Since in this limit ». is a perturbation on », all functions of » in (6.38) have Taylor expansions 
in the neighborhood of »,. Retaining terms in 6/8 if they affect the phase of y® but not other- 
wise, we find, for example, that 


2 

=— By cos 6+ 2(B cos (5) (6.44) 

To this degree of approximation (6.38) becomes 
and we finally obtain (6.45) 
y?(b) ~ (1+ R)(1—R) exp [—7 cos 6+ 276 cos 6,1, (6.46) 

where 
= 2B cos 6, = 

By Cos cos 9, (6.47) 


The quantity, 7,,, is the transmission coefficient of a ray incident on the dielectric surface from 
the interior at an angle 6, with the normal. If the transmission coefficient from the external 
region is designated by 7=1-+-R, then the sum of (6.35) and (6.46) has the form 


These two rays are depicted in figure 3. It is an elementary exercise in geometry to show that 
their amplitudes and phases are exactly those predicted by ray optics, provided the dielectric 
surface is nearly plane over the arc segment illuminated by the bundle included between rays 
land 2. Additional rays for this quasi-plane example may be written down by inspection from 
the further terms of (6.27). 
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JOURNAL OF RESEARCH of the National Bureau of Standards—D,. Radio Propagation 
Vol. 65D, No. 1, January-February 1961 


An Atlas of Oblique-Incidence Ionograms (A Digest)’ 


Vaughn Agy, Kenneth Davies, and Roger Salaman 


(June 24, 1960) 


A brief description is given of an atlas (published separately) of oblique-incidence 
ionograms obtained on two paths, 1,150 and 2,370 kilometers, 


The purpose of this note is to describe an atlas of 
oblique-incidence ionograms recently compiled by 
the authors [1].?. The atlas is intended to serve a 
twofold purpose: First, to provide a reasonably 
comprehensive introduction to this type of iono- 
gram for those workers who are not now familiar 
with them, and second, to present records which 
illustrate the various propagation characteristics of 
the specific paths used by the National Bureau of 
Standards. An analytical study of the technique 
and the experimental results was given by Agy and 
Davies [2]. To date, NBS sweep-frequency ex- 
periments have been carried out over two approxi- 
mately east-west paths [2]: Washington, D.C.—St. 
Louis, Mo. (about 1,150 km) [3, 4], and Washington, 
D.C.—Boulder, Colo. (about 2,370 km) [5]. A ver- 
tical incidence ionosonde was operated at the path 
midpoint in each case. 

The atlas is divided into four parts, namely: 
1. Introduction; 2. Washington—St. Louis routine 
records; 3. Washington—Boulder routine records; 
and 4. Washington—Boulder experimental records. 

The introduction contains a discussion of the 
equipment together with charts necessary for the 
identification of certain modes of propagation and a 
brief description of the records. 

The next two sections illustrate diurnal and sea- 
sonal variations, layer formation and disappearance, 


1 Contribution from Central Radio Propagation Laboratory, National Bureau 
of Standards, Boulder, Colo. 
2 Figures in brackets indicate the literature references at the end of this paper. 
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spread echo, MUF extension and several other phe- 
nomena of interest. The magnetic K figures are in- 
cluded to provide a basis for judging the effect of 
magnetic disturbance. The accompanying figure is 
a sample page from the section dealing with the 
Washington—Boulder path. 

The final section is devoted to various short term 
experiments over the Washington—Boulder path. 
Expanded time-delay and frequency scales were used 
as well as shorter pulse durations than in the previous 
cases. Equipment effects such as pulse width and 
receiver gain are considered. Detailed sequences of 
oblique ionograms during magnetic storms are shown 
together with the midpoint vertical-incidence iono- 
grams. A sequence of A-scan records is presented 
which shows how the field strength varied near the 
MUF during one sweep. 
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JOURNAL OF RESEARCH of the National Bureau of Standards—D, Radio Propagation 
Vol. 65D, No. 1, January-February 1961 


A New Approach to the Mode Theory of 
VLF Propagation’ 


James R. Wait 


(July 14, 1960) 


An attempt is made in this paper to present a concise derivation of the mode theory 


of VLF propagation. 


Taking note of the fact that the important modes for long-distance 


propagation are near grazing, suitable approximate forms of the wave functions are intro- 


duced at the outset, rather than at the end, of the analysis. 
for the influence of earth curvature in a relatively concise manner. 
Finally, numerical results for the attenuation 


earth’s magnetic field is also discussed. 


It is thus possible to account 
The influence of the 


and the phase velocity of the dominant mode are presented. 


1. Introduction 


Despite the extensive work [1—5]? on the wave- 
guide mode theory of VLF propagation, certain 
aspects of the subject are not yet resolved. The 
difficulty appears in finding suitable approximations 
to the spherical wave functions in the rigorous solu- 
tion. It is the purpose of this paper to discuss 
this problem in some detail. For a background of 
the subject the reader is referred to a recent com- 
prehensive paper on the general subject [6]. Here 
the solutions will be obtained from a somewhat 
different point of view. 


2. The Groundwave Field 


We will start by assuming that the source is a 
vertical electric dipole located on the surface of a 
smooth spherical earth of radius a, conductivity o and 
dielectric constant «. Spherical coordinates (r,6,¢) 
are chosen with the dipole located at r=a and 6=0. 
For harmonic time dependence the radial electric 
field component is written in the form 


sin 6)'/? (1) 
apart from a constant factor.’ In the case of an 
airless earth, in which the ionosphere is neglected, it 
is well known that V) may be written in the form 
[7,8] 


Vor 


e7 ite w,(ts—y) 
W,(ts) » (2) 


where 
ka 1/3 1/3 ka 1/3 
6, k(r—a), ig=( ) A, 


1 Contribution from Central Radio Propagation Laboratory, National Bureau 
of Standards, Boulder, Colo. 

2 Figures in brackets indicate the literature references on page 46. 

3 Equation (1) is normalized such that Vo approaches 2 for ka and A-0. 


and 


lew 
o+tew o+tew 


The coefficients ¢, are solutions of the equation 
wi (t)—qu (t)=0, (3) 


where w,(t) is an Airy integral and the prime in- 
dicates a derivative with respect to ¢. In terms of 
Hankel functions of order one third, 


w,(t)=exp (—2nt/3) (—mt/3)"? [(2/3) (—t)*”]. 
(4) 


The above formula for Vo is usually called the 
residue series representation for the groundwave 
field. It is often used as the basis for field strength 
calculations on a spherical smooth earth on the 
assumption that the ionosphere may be neglected 
or separately accounted for. It is valid when 
ka>>1 and r—a<<a. It is also required that the 
values of |¢,| for the important modes are not large 
compared with unity. Such approximations as these 
are certainly valid for VLF waves. 

For the purposes of the subsequent analysis Vo is 
written as a contour integral in the manner 


The contour encloses the poles at t=t, in a clockwise 
sense. The equivalence of eq (2) and (5) is easily 
verified on noting that 


wy’ (ts) —qwy (ts) = (ts— Wy (ts) 
since 
w;' (t)=twi(t), 
and 
wi (ts) =qui(ts). 
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3. The Sky Waves 


We'will now enclose the earth by a concentric 
reflecting shell located at r=a+h as indicated in 
figure 1. The electrical properties of this layer will 
not be specified at the moment, but it is assumed 
that an upgoing wave will be converted to a down- 
going wave. Thus after one reflection 


e~ ap, (t—y) 
is converted to 
A(t)e~‘"w,(t—y), 
where 


w,(t) =exp(2mi/3) (6) 


and A(¢) is unknown function of t. The boundary 
condition at r=a-+A is written in the form 


dy qiw(t—y) a (7) 
where 
> 1/3 
kh 
and ia=(F A; The quantity A; involves the 


properties of the layer beyond r=a+h; for the 
moment it is not given explicitly. F ormally A,=Z, Mo 
where Z is the radial surface impedance at r=a+h. 
Most generally A; (or Z) is a function of ¢ but usually 
it may be taken as a constant [6]. On identifying 
w(t—y) with w,(t—y)+A(t)w,(t—y), it readily fol- 
lows from eq (7) that 


wi (t—Yo) + qiwi(t—Yo) (8) 


A(i)=— W;(t—Yo) + qiwe(t—Yo) 


The downcoming wave characterized by the 
function w.(t—y) is now reflected at the ground and 
thus it generates a new upgoing wave of the form 

A() Bw (t—ye"™ 


- The boundary condition at r=a may be written 
dw(t—y) 
+ quwit—y) —0. (9) 


Then on identifying w(t—y) with the sum of the 
downcoming and the upgoing wave it follows that 


The process may be —_ any number of times. 
The resultant field can thus be written as a sum 
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(REFLECTing LAYER) 


(souRCE) 


(70 (GROUND) 
Figure 1, A sketch of the model employed in this paper. 


The surface of the ground and the lower edge of the reflecting layer are-con- 
centric spherical surfaces. 


tkad 
(11) 
where 
V= 2 V;, (12) 
and 
13 
for 7 even, while 
/2 ittey,(t—y) m 
14 


for j7 odd. Formally these are geometrical progres- 
sions so they can be summed. On interchanging the 
order of integration and summation, this leads to the 
integral representation 


e~ [wi (t—y) + A(t) w2(t— 


[w; (t)— qua (t) |[1—A(t) 
_ nf e~ (t—y) + A(t)w.(t—y)] 
w(t) — qui (t) + A(t) (t) — 
(15 


Now the contour is to enclose the complex poles which 
occur at ¢t=¢, where t, is a solution of 
1— A(t) B(t)=0. (16) 


The residue series representation for the total field 
is thus given by 


re: 


(20 (7,8, gy 

I 
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I 
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V= 


(17) 


As in the residue series formula for the groundwave 
field, eq (17) is valid for ka>>1 and h< <a. 

The interpretation of the preceding results is now 
discussed and certain simplifications made to facili- 
tate computation. If (—t)>>1, it follows that 


w,(t) — qwe(t) 
—1(—t)!?—q 


48/2 
/ 


(18) 


1/3 
On making the substitution (—)=(*) C and 


my. 
noting that iq~(#) A, it follows that 


(19) 


If Cis identified as the cosine of an angle of incidence, 
(C—A)/(C+A) can be recognized as a Fresnel 
reflection coefficient. Similarly, if (yp—t)>>1 


where C’ may be identified 


as the cosine of the angle of incidence at the iono- 
sphere. It is also noted that 


2h 1/2 


The factor (C’—A,;)/(C’+A,) is, of course, a Fresnel 
reflection coefficient referred to the bottom of the 
layer at r=a+h. 

The modal eq (16) which determines the coeffi- 
cients ¢, may thus be written in the approximate form 


(21) 


R,R, exp (22) 
where 
and 
t= 


ka 
(oo 


then fand [~2khC which cor- 
C+A; 


a to the modal equation for the flat-earth 
case [5]. 
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4. Evaluation of the Residues 


For numerical work it is desirable to express the 


function 
ra) 
E |, 


explicitly in terms of Airy functions. 
First, it is noted that 


w,(t)— qw,(t) —21(t,—¢) 
dé (Wilts) —Quilt) 


(23) 


where use has been made of the relations 

(1) w;’ (t)=twi(t) and 

(2) w,(t)w,(t) —w;(t)w.(t)=—2i which are valid for 
any value of ¢. Similarly, 


wi (t—Yo) + (t—Y) 
Ot Wo(t—Yo) + qiWe(t—Yo)_] 


—Yo— 


Thus 
re) 
E 
—2i(t,—@q’) 
~ [wi(tn) — — (tp) — qw2(t,) | 
4 Gi) [Wo(tn) — (25) 


(tn) (tn) |[We(tn—Yo) + JiW2(tn—Yo)}? 
Therefore, the complete residue series may be written 
V= (wx) 


[We(tn—Yo) + QiW2(tn—Yo) }? 

(26 


Numerical evaluation of the series is straightforward 
if tables of the Airy functions of complex argument 
are used. In the situation where 


(— tn) 1 


the Airy functions may be replaced by the first term 
of their asymptotic expansions. Under this approxi- 
mation and after considerable algebraic reduction, it 
is possible to write, for y=0, 


/2 


aR, R, 
oie, 


n=0,1,2.. 
aot 


| 
“43 
a 
4 
3 
& 
% " 


where 


C’—A, 


and 


(27) 


If the height of the source dipole is 23=7)—a and the 
observer is z=r—a, the F,, in the above expression 
should be replaced by 


F, (z)fn (Zo) 


where 


0 
0 
=R,-” exp 01422 


exp [ (C2422 (28) 


If 22/a<<C,? it is permissible to use 
exp (tkC,z)+R,'” exp (—ikC,z) 


and if R,~1 corresponding to a perfectly conducting 
ground, 


In(2)& cos kC,z. 


The value of C, to employ in the residue series is 
related, of course, to t, by 


Under the condition that (—t,)?>>1, the root- 
determining equation becomes 


where V is defined by 
Noting 
1/2 


~2kh if 
~2k(2ah)”? if 2h/a>>C3, 
it is seen that the second term within the square 


bracket on the right-hand side is small if kA is reason- 
ably large and if the losses are small. 
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5. Discussion of Various Approximations 


As a check on the algebra it may be noted that 
under the flat-earth approximation (i.e., 2h/a<C,”), 
the residue series formula may be written 


(2x0 vite 
(32) 
where 


p=a6, and For practical 
purposes, and 6,21 for n=1, 2, 3. 
Equation (32) agrees with an earlier analysis (5, 6]. 

For a first approximation it is probably permissible 
to neglect the curvature altogether. To improve on 
this it is important to retain curvature dependent 
terms in a consistent manner. For example, at VLF, 
the restriction C2>>2h/a is not as stringent as 


Thus eq (27) for V is probably not any better (if as 
good) as eq (32). To effect a real improvement one 
should retain the Airy functions of argument t, 
without approximation, but allow the Airy functions 
of argument t,—yo to be replaced by the first terms 
of their asymptotic expansions. This is justified 
since Y is reasonably large compared to unity in the 


VLF range. With this approach the root-deter- 
mining equation may be written 
(33) 
where 
w; (t)— (t) 
3 JIL 


If the ground is assumed to be perfectly conducting 
R,=1 and q=0. Then, on writing 


the mode equation reduces to 
arg Fy, (35) 
where now 


For long-distance propagation in the VLF range, 
only the mode corresponding to n=1 is of importance. 


(36) 
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Furthermore, since C is small, for a first approxi- 
mation 

wit)? 


also, 8 is essentially a constant and is approximately 
equal to 1/A;. Thus 


This formula should be a fairly good approximation 
if |C?| is somewhat less than 2h/a. For example, if 
B=0, corresponding to R;=—1, and if h/a=10-? 
the above may be written 


Cex (37) 


»0.041248 
C 0.01334. 

For h/A in the region from 2.5 to 5, |C?| is never 
greater than 5X 107°. 

Using perturbation methods, higher order approxi- 
mations for the solution of eq (35) may be obtained. 
— of such computations will not be described 

ere. 


6. Certain Extensions of the Theory 


The preceding analysis has referred specifically to 
a vertical electric dipole. If the source is a vertical 
magnetic dipole the results derived above are still 
applicable if the reflection coefficients R, and FR; are 
replaced by their counterparts 2? and &} for hori- 
zontal polarization. The only other change is that q¢ 
is now defined by 


Since [q] >1, under all practical conditions, 


(t) 


which is to be used instead of F,. The modes ex- 
cited by a magnetic dipole are of a TE (transverse 
electric) type whereas for an electric dipole they are 
of a TM (transverse magnetic) type. 


7. Application to Specific Ionospheric 
Models 


In the case of a homogeneous isotropic ionosphere 
with a sharp lower boundary, the Fresnel reflection 
coefficient for small C’ may be written 


C’—A 
fm 


(39) 


4l 


1 
where par: In terms of the plasma frequency wo 
t 


and the collision frequency »v 
B& (1—tw,/w) (—iw,/w)—/? 


where w;=w,"/y under the condition that wv. 
Unfortunately, the presence of the earth’s mag- 

netic fields renders the ionosphere anisotropic. The 

reflection coefficient is thus a matrix of the form [1] 


Ry 

where the coefficients Ry, AR), and ,R, indi- 
cate the complex ratio of a specified electric field in 
the wave after reflection to a specified electric field 
in the wave before reflection. The first subscript 
denotes whether the electric field in the incident 
wave is parallel (;) or perpendicular (,) to the plane 
of incidence and the second subscript refers in the 
same way to the reflected wave. It is understood, 
of course, that [FP] is a function of C’, the cosine of 
the angle of incidence at the ionosphere. Following 
an earlier analysis [6], the generalization of eq (33), 
to include the influence of anisotropy, may be 
expressed in the form 


[1—RyR 
x [1 — RR vie] 


— Ri Ry RF vic— (40) 
where R,, Rt, F,, Ft, and V have been defined. 
Strictly speaking, equation (40) is only valid when 
there is perfect azimuthal symmetry about the source 
dipole. 

It may be noted that if the conversion coefficients 
,2, and ,R, vanish, the modal equation splits into 
two parts which are 


1—)RyR, =0 (41) 


and 


1—, Ry (42) 
These two uncoupled equations give the solutions for 
the TM and TE modes, respectively. In the more 
general case the modes are coupled. However, for 
low-order modes and particularly for the dominant 
mode (n=1), the coupling is extremely small and a 
good first-order approximation is obtained on the 
assumption that coupling may be neglected. Fur- 
thermore, the TE mode is only of academic or 
secondary interest since it is not excited to any ex- 
tent by a vertical antenna. Thus the constants of 
the first-order mode may be obtained from eq (37) 
with an appropriate value of 8. 
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For propagation in the magnetic meridian or for 
propagation in polar regions the QL (quasi-longitu- 
dinal) approximation [9] to describe conditions within 
the ionosphere is valid. The refractive index yu, in 
this case, is given by [1] 


=1—1(w,/w) exp (43) 
where 
w, longitudinal component of gyrofrequency 
tan ¢,=—= 
v collision frequency 


and Thus for highly oblique 
incidence (i.e., C is small) it may be shown that [6] 

Re (44) 
where 


p= 
CoC. 


where yo and u, denote the two values of x», and 
Co=(1—1/uo)# and C,=(1—1/u,”)3. 


The QL approximation mentioned above is only 
valid when [10] 


) 


where w, and wz are the longitudinal and transverse 
components of the (angular) gyrofrequency. As 
mentioned above, this condition is satisfied for 
propagation in the magnetic meridian and/or for a 
steeply dipping field. Clearly, this conditien is 
violated when the transverse component of the 
earth’s magnetic field is large such as for propagation 
around the magnetic equator. The case of a purely 
horizontal and transverse field has been considered 
by Barber and Crombie [11]. Using their results, 
which are applicable to a sharply bounded ionosphere, 
it may be readily shown that [6] 


=0 (45) 


where now 


and 


For east-to-west propagation, y is positive, while 
for west-to-east propagation, y is negative. In the 


present situation it is convenient to define an angle 
42 


such that 
y=-+L|tan ¢7| for east-to-west propagation 
= —L|tan ¢7| for west-to-east propagation. 


For propagation at any angle with respect to the 
earth’s magnetic field and for an arbitrary angle of 
dip J, it is necessary to seek more general solutions. 
First, one must solve a quartic equation for the 
refractive indices (ordinary and extraordinary com- 
ponents). Then on matching tangential field com- 
ponents at the air-ionosphere boundary, an explicit 
expression for the reflection coefficient may be ob- 
tained [12]. Using this quite general approach 
Johler and Walters [13] have obtained extensive 
numerical results for the coefficients Ry, 
and ,F) for various (real) angles of incidence and 
range of values for the electron density, N, collision 
frequency v, dip-angle J and the azimuthal direction 
of the wave ¢, (measured clockwise from the north). 
Formally, the expressions for these reflection and 
conversion coefficients may be substituted in eq (13) 
which, in principle, can be solved for the sabe in 
the general case. For VLF and at highly oblique 
angles of incidence (e.g., 82°), the quantity |. Xx 
.F))| is of the order of 10~* or 10~* and thus coupling 
between the TE and TM modes is negligible. At 
least this is true for the range of values of the iono- 
spheric parameters considered here. Thus the un- 
coupled modal equation may again be used for an 
adequate first approximation after having determined 
the appropriate value of 8. 


8. Presentation of Results 


Using the above approaches numerical results 
have been obtained for the complex values of C for 
the dominant mode.* From this quantity, it is a 
simple matter to compute the corresponding value of 
S which describes the propagation in the horizontal 
direction. For example, the attenuation per unit 
distance is proportional to the quantity—H Im S 
and the phase velocity relative to ¢ is 1/Re S. In 
the present paper, these are the quantities of prime 
interest since they are adequate to describe the 
propagation characteristics of VLF radio waves at 
large distances from the source (i.e., >2,000 km). 
At shorter distances higher order modes should be 
considered. Their influence has been discussed 
elsewhere [5]. 

In figure 2 the attenuation coefficient, for the 
dominant mode (n=1), is plotted as a function of 
H(=h/d) for a curved earth where h/a=0.01. The 
upper- and right-hand scales in this figure and in 
many which follow, are a frequency scale and at- 
tenuation scale in decibels per 1,000 km of path 
length applicable to a height of 70 km. The quasi- 
longitudinal approximation is employed and ¢, is 
either 60° or 0°; the latter, of course, corresponds to 
no magnetic field. The constant B on the curves is 


4 For sake of brevity, the subscript 1, to indicate the first mode, is dropped in 
what follows. 
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defined by are not expected to be accurate to more than a few 
we parts in 
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Ficure 2. Attenuation coefficient,—H Im 8, as a function of 
H (=h/A) for the dominant mode. 


The upper and the right-hand scales are for h=70 km. _ 


The four B values on the curves (0.02, 0.05, 0.1, and 
0.2) correspond to w, values of 13, 5.4, 2.7, and 
1.3X10°, respectively. For a given value of B, it is 
seen that the attentuation depends on ¢,. For large 1.000 
values of B (i.e., poorly conducting ionosphere) the 
influence of the terrestrial magnetic field is to lower 
the attenuation whereas the converse is true at 
smaller values of B. Since B is in the range 0.05 to 
0.01 typically for daytime conditions, the influence 
of the magnetic field is small under the quasi- 
longitudinal approximation. 

. = The relative phase velocity of the dominant mode 
is shown in figures 3a and 3b for ¢,=0 and 60°, 0.998 
respectively. The (upper) frequency scale is for a 
height of 70 km. From these it appears that the 
influence of the magnetic field is quite small. It is 
rather interesting to note that for a frequency of 0.997 
14 ke/s the phase velocity is very near, c, the velocity 
of light.’ It should be remembered, of course, that 
the phase velocity is referred to the surface of the 


0.999 


ground in this particular analysis. N | | | | I | 
Since one is interested in relatively small deviations acai 3 4 5 6 

of 1/Re (S) from unity, considerable care is needed in 

approximating the wavefunctions involved. The H 


results shown in figures 3a and 3b for a curved earth | Ficures 3a andb. The relative phase velocity, 1/Re S as‘a 


unction of H for the dominant mode. 
5 Recent experimental data appear to confirm this point (M. L. Tibbals, f f fi 
private communication). The upper-frequency scale is for h=70 km. 
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The calculation of the mode characteristics for a 
magnetic field of arbitrary direction is, of course, 
quite involved. However, approximate results can 
be obtained if coupling between TE and TM type 
modes is neglected [6]. Actually, in the case of a 
purely transverse field, the coupling is zero anyway. 

In the following the influence of an arbitary mag- 
netic field is demonstrated by calculating the ratio, 
P, of the attenuation with a magnetic field to the 
attenuation without the magnetic field. Various 
values of the dip-angle J, and the azimuthal 
direction of propagation, ¢,, are chosen. The ratio, 
Q, of the phase-velocity deviations with and without 
a magnetic field is also considered. Specifically, these 
ratios are 


___ Attenuation with magnetic field Im 8 
~ Attenuation without magnetic field [Im 
(48) 


and 


Q= Phase velocity deviation with magnetic field 
~ Phase velocity deviation without magnetic field 
(49) 


1 
| 
Re S 2a 
where J ~ 


| 
E =0 


To obtain numerical values of P and Q, computed 
data from J. R. Johler [14] were requested for the 
reflection coefficient |), for a real angle of incidence 
of 82° and at frequencies of 10, 12, 14, 16, 20, and 
22 ke/s. The appropriate value of 8 was then 
calculated from 


(50) 
where C’=cos 82°. The electron density values for 
the sharply bounded model were 10° and 310° 
electrons/cm*, respectively. The collision frequency 
was taken to be 2X10’ and the strength of the 
earth’s magnetic field was 0.5 gauss. The angle 
¢ which defined the direction of propagation 
measured clockwise from north assumed the values 
0°, 60°, 120°, 180°, 240°, and 300°. As a check on 
the method, values of 6 were calculated directly 
ay eq (40) which is applicable to ¢,=90° and 180° 
or J=0. 

The results for the attenuation ratios are shown 
in figures 4a and 4b for a dip-angle J of 0°, in figures 
5a and 5b for an J of 45°, and in figures 6a and 6b for 
an I of 84.3°. 

These results would indicate that nonreciprocity 
is most evident in propagation around the geomag- 
netic equator. North-to-south or south-to-north 
propagation appears to be almost reciprocal. It is 
also evident that for this range of electron densities, 
the nonreciprocity for attenuation is more marked 
for the smaller densities. It is also apparent from 
these results that the nonreciprocal effects are 
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Ficures 4a andb. The ratio of the attenuation rates with 
and without a (transverse) terrestrial magnetic field, 
¢c is the direction of propagation with respect to geomagnetic north. 2 
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Figures 5a and b. The ratio of the attenuation rates with and 
without a terrestrial magnetic field for a dip angle of 45°. 


1.2 
F 1,0 
0.9 
0.8 
0.7 
0.6 
P 
P 14 N= (10)? 
1=0° f= 22kc 
Vv = 
7 \ 
hi 
1.0 
0.9 \ 
os 
(b) : 
0.7 
a 
F210 k 
10 10 kc 
f=22kc 
0.8 (a) - 3 ( 
07 
P 1.3 
=22kc 
45° SY 2 
lo 
09 
(b) 
0.7 i 
| 
Qa 
F 


diminished as the magnetic field approaches the 
vertical. Qualitatively the curves are very similar 
and there is a strong suggestion that nonreciprocity 
is proportional to the component of the earth’s 
magnetic field which is horizontal and transverse 
to the direction of propagation. 
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Figures 6a and b. The ratio of the aitenuation rates with and 
without a terrestrial magnetic field for a dip angle of 84.3°. 
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Figures 7a and b. The ratio of the phase-velocity deviations 
with and without a (transverse) terrestrial magnetic field. 
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Figures 8a and b. The ratio of the phase-velocity deviations 
with and without a terrestrial magnetic field for a dip angle 


of 45°. 
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Ficures 9a and b. The ratio of the phase-velocity deviations 
with and without a terrestrial magnetic field for a dip angle 
of 84.3°, 


The corresponding curves for the ratio Q of the 
phase velocity deviations are shown in figures 7a and 
7b for 7=0, in figures 8a and 8b for 7=45°, and in 
figures 9a and 9b for J=84.3°. It is rather interesting 
to see that these are not symmetrical about the 
north-south geomagnetic axis. It is also evident 
that the nonreciprocity in the phase velocity 
becomes greater at the lower frequencies. It is 
noticed that Q is never more than 4-percent greater 
than unity. This corresponds to a change of phase 
velocity of only a few parts in 10°. 
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9. Concluding Remarks 


In this paper no direct reference has been made to 
experimental data. In another paper [15] the 
experimental data on attenuation rates at VLF in 
daytime [16] were compared with calculations using 
formulas from the present paper. The agreement 
between theory and experiment was quite good and 
the nature of the nonreciprocity between east-west 
and west-east propagation was satisfactorily ac- 
counted for. Unfortunately, the experimental data 
on phase velocity is quite sparse. It would be 
interesting to compare the computed results given 
in this paper with absolute measurements of phase 
velocities at various frequencies over a given path. 


The author thanks Kenneth Spies for carrying out 
the calculations and Mrs. Eileen Brackett for her 
assistance in preparing the manuscript. 
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10.2. List of Symbols 


(r,0,¢) =spherical coordinates of observer, 
a=radius of the earth, 
h=height of reflecting layer, 
k=2zx/wavelength, 
@ is central angle, 
x= (ka/2)!4 @, 
y= (2/ka)'8 (r—a), 
(ka/2)1* A, 


a=( teow teow 
o+tew 
w,(t) is Airy function of first type defined by eq (4), 
w2(t) is Airy function of second type defined by eq 


6), 

tis argument of Airy function, 

Ai=Z/n, 

no 1207, 

Z is surface ee of ionospheric reflecting 
layer at r=a-+h, 

iqi= (ka/2)'4 Aj, 

yo= (2/ka)' kh, 

w(t) is some linear combination of w;(t) and w(t), 
A(t) is defined by eq (8); it is related to the iono- 
spheric reflection coefficient, 
Bid is defined by eq (18) ; it is related to the ground- 
reflection coefficient, 

C= (2/ka)14(—1t)12; C may be geometrically in- 
terpreted as the cosine of the angle of in- 
cidence at the ground, 

C’ = (C?+ 2h/a)!2; C’ may be geometrically inter- 
preted as the cosine of the angle of inci- 
dence at the ionosphere, 

R, and FR; are reflection coefficients at the ground and the 
ionosphere, respectively, 


V is a dimensionless quantity proportional to the 
vertical electric field, 
V; is defined by eq (13); it may be interpreted 
loosely as the j’th order skywave, 
t, is the n’th root of eq (16); it determines the 
attenuation and phase of the n’th mode, 
p=aé; it is the great circle distance between 
source and observer, 
5n is defined by eq (32); it is a coefficient in the 
mode sum, 
F, is defined by eq (34); it is a wave correction to 
the ground-reflection coefficient, 
@r= 
wo) is the (angular) plasma frequency, 
v is the collision frequency, 
B is defined by eq (39); it characterizes the iono- 
reflection coefficient Ri, 


me | this is a matrix which describes 
reflection from an anisotropic layer, 
uw is defined by eq (43) ; it is a double-valued index 
of refraction for the ionosphere subject to 
the QL approximation, 
wo and yu, are the two components of yp, 
wz is the longitudinal component of the gyrofre- 
quency, 
wr is the transverse component of the gyrofre- 
quency, 


Dol 


wrh 
H=h/d 


[RJ= 


=tan- 
P and Q are defined by eqs (48) and (49), 
I=dip angle of the earth’s magnetic field, 
¢, is the azimuthal direction of propagation. of the 
wave; ¢.=0 corresponds to propagation 
from ‘south-to-north, ¢a=90° corresponds 
to propagation from west-to-east, etc. 
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East-West Effect on VLF Mode Transmission 
Across the Earth’s Magnetic Field ' 


D. Dobrott and A. Ishimaru ” 


(July 5, 1960) 


The effect of a constant transverse magnetic field is analyzed for the propagation of VLF 
electromagnetic waves about the earth. These waves are considered to be propagating by 
modes in a parallel plate waveguide. The lower boundary of the guide is considered to be a 
perfect conductor, while the upper boundary is assumed to be a sharply defined, semi-infinite, 
homogeneous plasma with a constant magnetic field applied. The source of these waves is 
an infinite number of short, horizontal, cophased dipoles, uniformly spaced parallel to the 
constant magnetic field vector. Admittance relations are derived for the upper boundary of 
the guide by considering the waves to be of grazing incidence. It is found that these admit- 
tance parameters depend on the direction of propagation. These admittance parameters are 
plotted versus frequency for various heights. Expressions for eigenvalues and the eigen- 
functions are found as functions of the admittance. By employing a two dimensional 
Green’s function, the amplitudes of the various modes due to the dipole source are evaluated. 
It is noted that the propagation constants differ depending on the direction of propagation, 
thus offering an explanation of the east-west effect of VLF transmission. A numerical 
example is calculated and field strength versus distance values are found to correspond to 
some experimental results. 


The solution of the radiation from a dipole in this 
parallel plate waveguide with anisotropic boundary 
conditions is, in general, difficult. However, one 


1. Introduction 


The mode}theory of VLF radio transmission has 


evolved from the study of electromagnetic wave 
propagation between the earth and the ionosphere. 
Wait [1] * and others [2, 3] have solved the VLF 
electromagnetic problem using mode theory and as- 
suming a sharply bounded ionosphere, while neglect- 
ing the earth’s magnetic field. Friedman [4, 5] has 
made a modal analysis which differ in method from 
that of Wait, but arrives at similar conclusions. 
Friedman’s analysis utilized an upper boundary of 
scalar impedance Z. 

All of these solutions have the same propagation 
characteristics regardless of the direction of propaga- 
tion. However, Crombie [6] has noted evidence that 
VLF propagation differs depending on whether the 
propagation is from east to west, or from west to 
east. He also indicated that this so-called ‘“east- 
west effect”? may be partly due to the earth’s mag- 
netic field. Recently, Wait [7, 8] and Crombie 
{[6, 9] have shown analytically that the rate of at- 
tenuation for east-west propagation behaves as had 
been indicated by the latter author. 

This paper attempts to show analytically the east- 
west effect on the basis of mode theory and iono- 
sphere-freespace boundary admittance. Because of 
the presence of the earth’s magnetic field, the admit- 
tance of the sharply bounded ionosphere is not a 
scalar, but a tensor. Further, because of these mag- 
netic effects, the tensor depends upon the direction 
of propagation. 


1 Contribution from Boeing Scientific Research Laboratory, Seattle, Wash. 

? Present address of Professor A. Ishimaru: Department of Electrical Engi- 
neering, University of Washington, Seattle, Wash. 

§ Figures in brackets indicate the literature references at the end of this paper. 
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may solve an analogous two-dimensional problem. 
This simplifies the arithmetic while maintaining the 
qualitative structure of the three-dimensional case. 

The two-dimensional problem corresponds to the 
source of electromagnetic waves being an infinite 
number of horizontal, cophased, short dipoles, uni- 
formly spaced parallel to the direction to the earth’s 
magnetic field. This situation is displayed in figure 1. 

In the sections that follow, the radiation. from 
these sources are examined in detail. A numerical 
example is shown also. 
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Fiaure 1. Cophased sources for TM waves. 
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Further, the nature of these sources produces only 
transverse-magnetic waves. Considering these facts 
and eq 2.1, Maxwell’s equations in the plasma be- 


2. Dielectric Tensor of Magnetically In- 
fluenced Plasma 


The permittivity of an ionized gas may be deter- | ©°™®: RS 
mined by examining the ion motion under the influ- ee 
ence of electromagnetic waves. Certain restrictions OH, OF, _ ou, H EB 
must be placed on this gas if the evaluation of the Or Og J Molin 2 
equation of motion is to remain simple. 3 
For the purpose of this analysis, the gas will be oH. fe 
considered to be a homogeneous plasma. That is, = —jweo(e2,—jnE,), & 
the gas is uniform throughout and the number of oz = 
positive charges is equal to the number of negative e f 
charges. The plasma exists in a state where colli- . ee Fa 
sions are relatively few. Since the positive and neu- nE,+ (3.1) 
tral particles have more mass than the electrons, é 
their mobilities are correspondingly smaller and it 
shall be considered that their velocity is zero. Fur- Solving eqs 3.1 in terms of the longitudinal electric © 
ther, the number of electrons per cubic meter N, | field to obtain the following equation & 
avior of a volume element of gas is the same as that 2 (&—n’) a 
of an “average” electron obeying the laws of classical (3.2) I 


co 


dynamics. These are the assumptions of the ‘‘Mag- 
neto-Ionic Theory” [10]. 

In the special case considered in figure 1, where a 
magnetic field is considered perpendicular to the di- 
rection of propagation, the ionosphere’s relative di- 
electric tensor is given by, [11] 


The solution to this equation is of the form 


E,=K,< (3.3) 


€2 0 —jn where p and y are, in general, complex and the plus . 
tat 0 or the minus sign designates a wave traveling in the | ® 
“ (2.1) negative z-direction or positive z-direction respec- 
0 ‘ tively. 
Mm . Substituting eq 3.3. into 3.2 one obtains the charac- 4 
. teristic equation of 3.2 3 
(2.2) 
(2.3) 
Jol? + (gw+r)?] fi 
In the free-space region, Maxwell’s equations may — sg 
aint be solved for the transverse-magnetic case in terms 3 
= a “° (2.4) | of the longitudinal electric field yielding the following 
In expressions 2.2, 2.3 and 2.4 w, is the plasma fre- > Oe 
quency, v is the collision frequency, w, is the cyclotron | atsote E,=0. (3.5) 
frequency, and w is the frequency of the electro- Ou 02 


magnetic wave propagating in the plasma region. 

sing 2.1, one is able to solve for the electromag- 
netic fields in the plasma region. The solution of the 
fields is required to examine the conditions at the 
plasma-freespace interface. 


The solution for the homogeneous wave equation 
depends upon the boundary conditions, and since the 
tangential electric field must be zero at a perfect 
conductor, the solution is of the form, 


3. Admittance Tensor of Plasma at 
Interface 


E,=K;, sin que*™ (3.6) 


In order to solve for the electromagnetic fields for 
the problem described by section 1, one must first 
consider the restrictions placed upon the problem 
by the particular geometry. Since the sources are 
considered to be uniformly distributed along the y- 
direction axis, there is no y-variation for the fields. 
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where g and y are complex. 
Equation 3.6 has a corresponding characteristic 
equation given by, 


(3.7) 


| 
| 

| 

T 
of 

| 

= 


& 5 


From eqs 3.1 and 3.3 the fields in the plasma region 
are related by the following relations 
for the +-2-direction, 


ey—jpn 
1 —jnY¥—e.p 
(3.8) 
for the — z-direction 
—exy—jpn 
Jweo 
(3.8) 
(cont’d) 


From the two tangential fields E, and H,, an admit- 
tance relationship may be defined as 


H, 
YY 


(3.9) 


From the two expressions of eqs 3.8 and relation 3.9, 
one may obtain two expressions for the admittance 
of the plasma boundary. 


(3 0) 
where Y,-» and Y,_, are the relative admittance 
parameters for east-west and west-east propagation, 
respectively. 

If grazing incidence is considered, upon the plasma- 
freespace interface, it is found that yk, where the 
sign has already been determined by the form of eq 
3.6. This corresponds to the far-field situation. If 
this condition were not imposed, the admittance 
would be a function of the free-space eigenvalue gq. 
Thus one obtains the following two expressions for 
the admittance of the plasma boundary, 


e-w 
w-e ep nk (3 ) 


The following restriction must be placed upon p, 


Re (p)>0, Im (p)>0. (3.12) 


_ With these restrictions in mind, the final expres- 
sions for the admittances in terms of the constituents 
of equation 2.1 are 


€2 


(3.13) 
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A numerical solution of these admittance parame- 
ters for height versus frequency is shown in figures 2 
and 3. The values of collision frequency and elec- 
tron density have been obtained from papers by M. 
Nicolet [12] and R. E. Houston [13], respectively. 
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Figure 2. Real part of admittance (Y,2). 
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Fiaure 3. Imaginary part of admittance (Y12). 
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4. Eigenvalues and Propagation Constant 
In order to solve for a general expression for the 
eigenvalues in terms of the admittance matrix for all 
waves, one must also solve Maxwell’s equations for 
the transverse-electric waves. This solution, in terms 

of the longitudinal magnetic field H,, is given as 
H,=Y )K3e*™ cos qz. (4.1) 


Recalling eq 3.6 and Maxwell’s equations for TM 
and TE one can obtain, 


Y sin gx (4.2) 
and 
cos qa. (4.3) 


Relating these four fields at z=a by the admittance 


matrix Y, 
H, Yu TE, 
Yn YooJ LE, 
Substituting the solutions 3.6, 4.1, 4.2, and 4.3 into 
4.4 and combining terms 


(4.4) 


(x 12 SiN cos a) sin 
jk. 
(Y sin ga) (cos sin K; 
(4.5) 
For 4.5 to be true, the determinant 
sin cos a) sin a) 
=0. (4.6) 


(Y sin ga) (cos sin a) 


Evaluating the determinant 4.6, the following general 
expression for the eigenvalues is obtained in terms of 
admittance parameters, 


(4.7) 


A, tan ga+cot 


where A,= Yu Yi2 Ya. 
For the special parallel plate case under considera- 


tion, the only admittance parameter needed is Yj. 
Thus eq 4.7 becomes 


cot ga= (4.8) 


Rewriting relationship 4.8, g=qn,n=1,2,3... 


2 
( a 1 (4.9) 


If the term q,/k Yj. is large as is the case for the 
“electric” wall, the following approximation for the 
eigenvalue g, may be obtained. 


_ jka 
For the ‘‘zeroth mode” case, that is n=0, 
qk 1/2 
? (4.11) 


and if n is greater than zero, g, may be approximated 
nm in the right-hand side of 4.10 yielding, 


(4.12) 
The propagation constant 
(4.13) 


may be obtained from the results of 4.11 and 4.12. 
Restrictions on 4.13 are that 


Re >0, Im (7,)>0; 


when 4.13 is evaluated for n=0, the attenuation 
constant, Re(y,)=a, is very small, and for n>2, 
the phase constant, Im (yn)=8,, is small represent- 
ing evanescent modes. 

One can see that if Y;.—©, as in the case of a 
perfect conductor, expression 4.12 reduce to the 
familiar result, nz/a. Expressions 4.11 and 4.12 
are used, along with the data of Houston and Nicolet, 
to plot mode decibei loss versus distance on figure 4 
for n=0, 1,2. The numerical values for the admit- 
tance used were 


(4.14) 


Y g-w=22.5—37.5), (4.15) 


(4.16) 
This corresponds to a radian frequency of w=2 104 
and a height of 90 km. 

A synopsis of experimental results by Crombie [6] 
indicates that the ratio of electric field strengths 


W-—E 
(4.17) 
which is certainly true for the numerical case cal- 
culated and graphed on figure 4. 
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5. Mode Solution—Source Present 


With the source of current moment JL eq 3.5 
becomes, 


(2+ (5.1) 


where 6 (x—b) and 6 (2) are Dirac ‘‘delta-functions.”’ 
This problem is analogous to the Green’s function 
problem if one relates the longitudinal electric field 
to the Green’s function as L,=—jouILG, the cor- 
responding Green’s function problem may be solved, 


2 2 


The boundary conditions for the Green’s function 
are the same as for the electric field F, 


(G’+Y@)|,-.=0. (5.3) 

It has been shown in section 4 that the solution 
to the homogeneous wave equation for the longi- 
tudinal electric field and hence the Green’s function 
is of the form 


G= > sin (5.4) 
0 


Thus if 5.2 is true we must be able to expand the 
delta-function in terms of eigenvectors, [14, 15] 


B,, sin (5.5) 
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Multiplying through by sin q,,x and integrating over 
the range of x, an expression for B,, 


B 4q, SiN 


"~(2qnd—sin 29,4) 


Making the substitution 


4qn 
"a= q,a—sin 


where 7, is often called the normalization factor. 
Integrating 5.2 over the source with respect to z 
yields 

_ SiN gud 
(5.8) 


and 


G=> 24 


gj 
Yn(2Qnd—sin 2q,0) sin sin 


(5.9) 


Therefore the expression for the electric field is given 
by 


= qopo2lLgn 
n(2q,4—sin 


sin g,5 sin gp. 
(5.10) 


The other fields may be evaluated for Maxwell’s 
equations for the TM case. 
The expression for the vertical component of the 
electric field with a ‘‘normalized” current is given as 
2 sin g,b 
> 


2gn4—sin 29,0 


cos gna. 


If the transmitter is in an airplane at a height of 
a/4 and the receiver is on the ground, 5.11 becomes, 


For the first three modes 5.12 is of the form 
where 
ga) 


2q,a—sin 2q,4 


The magnitude of the field strength at the receiver 
is given by the expression, 


|E,|= {|Co|e- 202+ |C, 2022 
cos [(Bo— 2+ (C.—¢,)] 
cos Be) 2+ (C,—C2)] 


cos [(8,— Be) pus 
(5.15) 


B : 
| 
3 
: 
| | 
| 
| 
| 
d 


where |C,,| and C, are the magnitude and argument 
of C,, respectively. 

Equation 5.15 is plotted on figure 5. These 
results depend not only on attenuation, but also on 
field strength magnitude which in turn depends on 
the direction of propagation and the height of the 
transmitting antenna. This corresponds qualita- 
tively with the far field experimental results men- 
tioned by Crombie [6]. 
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Figure 5. Field strength versus distance. 


6. Conclusion 


The effect on east-west, west-east propagation of 
electromagnetic transverse magnetic waves across 
the earth’s magnetic field has been analyzed as a 
boundary value problem. Certain idealizations were 
made to solve this problem. Admittance relation- 
ships were derived from field relations existing in a 
homogeneous anisotropic plasma, and it was found 
that the admittance parameters depended upon the 
direction of propagation. These parameters were 
plotted versus frequency for various heights. Ex- 
pressions for the eigenvalues were derived. The 
attenuation over distance was plotted for a special 
case. It was found that the ratio of the electric 
field for the individual modes of west-east/east-west 
propagation was greater than one. This corresponds 
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to some experimental data. A two dimensional 
Green’s function was derived and a numerical calcu- 
lation of the corresponding electric field strength was 
plotted versus distance. 
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Magneto-Ionic Propagation Phenomena in Low- and Very- 
Low-Radiofrequency Waves Reflected by the Ionosphere’ 
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LF-VLF ionosphere reflection coefficients which illustrate the dependence of the ampli- 
tude and phase of the reflected wave upon the direction of propagation relative to the direc- 
tion of the earth’s magnetic field are presented. The calculations are based on a plane, 
sharply bounded, model ionosphere with plane wave excitation, but employ full use of the 
magneto-ionic formulas for complex directions of propagation in the ionosphere such that the 
influence of the earth’s magnetic field in the different directions of propagation is demon- 
strated. A special table of values applicable to VLF is presented. 


1. Introduction 


The classical magneto-ionic theory implies a de- 
pendence of the electromagnetic waves on the direc- 
tion of propagation relative to the direction of the 
earth’s magnetic field. 

The application of the classical magneto-ionic 
theory to the evaluation of low- and very-low-radio- 
frequency reflection coefficients for a sharply bounded 
model ionosphere by using an assumed electron 
equation of motion, together with Maxwell’s equa- 
tions was investigated in a previous paper [1]? 
The importance of the orientation of the earth’s 
magnetic field vector, H,,, was discussed in consider- 
able detail, and the effect of various orientations of 
the vector (relative to the direction of propagation) 
on the precise value of the reflection coefficient was 
noted. This paper further pursues the subject of 
“directional effects’? at low- and very-low-radio- 
frequencies. The forces acting on the electron in the 
ionosphere excited by an electromagnetic field are 
described by an electron equation of motion, Johler 
et al. [1] (1960). A particular term of this equation 
poe(VXH,) called the “Lorentz force,’ Lorentz [2] 
(1906) relates the action of the earth’s magnetic 
induction, 4/7, or an electron (charge e) traveling at 
a velocity, V, in a medium of permeability u,, to the 
direction of propagation. This dependence upon 
direction of propagation causes some interesting and 
predictable phenomena in the low and very low part 
of the radio spectrum. 

The effect of the collision frequency, v, on the 
magnitude of abnormal components (i.e., electric- 
magnetic coupling coefficients, Tem, Tme) is also in- 
vestigated. A large angle of incidence, ¢;=80°, 82°, 
is considered and the resultant intensities of both 
"1 Contribution from Central Radio Propagation Laboratory, National Bureau 


of Standards, Boulder, Colo. “ 
2 Figures in brackets indicate the literature references at the end of this paper. 
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normal and abnormal components are noted at 
various directions of propagation. 


2. Theory of the Reflection Coefficient 


The essential nature of low- and very-low-radio- 
frequency waves propagated via the ionosphere can 
be described in terms of a reflection coefficient. 
Reflection coefficients or propagation constants 
(complex index of refraction of the ionosphere) for 
complex angles of incidence, ¢;, figure 1, such as 
would be required for ‘‘mode-type’’ calculations, 
Wait [3] (1960), can be readily evaluated from the 


Fiaure 1. Coordinate systems. 
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results of the previous paper, Jobler et al. [1] (1960) 
without analytic modification of the formulas. How- 
ever, for naive discussion purposes in this paper, a 
real angle of incidence, ¢;, will be employed. It 
should be noted, however, that the wave normal 
direction, ¥, relative to the direction of the earth’s 
magnetic field vector, Hm, in the ionosphere is never- 
theless complex [1]. The direction, 6, relative to 
the vertical, 2 direction, figure 1, is also a complex 
number as a consequence of Snell’s law: sin ¢;=7 
sin 6. The propagation downward from the iono- 
sphere (i.e., wave traveling from the ionosphere into 
the negative z, figure 1, region) is defined by four 
reflection coefficients, T.., Tem, Tme, Timm, Which relate 
the reflected radiation to the primary or incident 
radiation. The reflection coefficient, T,., refers to 
the vertical electric polarization of the incident plane 
wave and a similar vertical electric polarization of 
the reflected wave. The coefficient, T.,, describes 
the generation of the abnormal component by the 
incident vertical polarization (vertical electric- mag- 
netic coupling). Similarly, Tm, refers to the inci- 
dent horizontal electric polarization and the corre- 
sponding reflected horizontal electric polarization. 
Also, the abnormal component generated by hori- 
zontal electric polarization (vertical magnetic-elec- 
tric coupling) is described by the coefficient, T'ne. 
The reflection coefficients can therefore be defined 
(see fig. 1), 


Ey, _ Ey, 


where subscripts i and r refer to incident and re- 
flected waves, respectively. Expressions for these 
coefficients were developed in the previous paper [1]. 
In addition to the geometric parameters, figure 1, 
the reflection coefficient is completely defined by 
the magnetic intensity H,,=|H,,| gauss, the electron 
density N(el/em*), and the collision frequency 


3. Discussion 


A comparison of the reflection coefficients of the 
model ionosphere, |Zem| and phase, arg 
arg Tem, figure 2, and amplitude, |7mm|, |Tne| and 
phase, arg Tnm, arg T., figure 3, for north-east 
propagation, magnetic azimuth, ¢,=45°, with the 
corresponding reflection coefficients for south-west 
propagation, magnetic azimuth, ¢,—225°, is illus- 
trated. It is quite evident that a nonreciprocity 
exists in the propagation (i.e., the expected value 
of field intensity would not be the same if transmit- 
ter and receiver were interchanged), since neither 
normal nor abnormal components are precisely the 
same in both directions except at zero frequency. 
Indeed, the discrepancy is appreciable not only at 
LF(<300 ke/s) but also at VLF(<30 ke/s). 

The effect of the magnetic azimuth, ¢,, (reckoned 
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Figure 2. Comparison of north-east propagation with south- 
west propagation (vertical polarization) illustrating non- 
reciprocity with the aid of the model ionosphere reflection 
coefficients, amplitude, |T|, and phase, arg T. 
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FicurE 3. Comparison of north-east propagation with south- 
west propagation (horizontal polarization) illustrating non- 
reciprocity with the aid of the model ionosphere reflection 
coefficients, amplitude, |T|, and phase, arg T. 


N=870, v=4 (10°), 6j=75.08°, Hm=0.5, 1=60°, 2=45°, 225°. 


as clockwise angular measurement from magnetic 
north, fig. 1) is illustrated at a ‘‘temperate” mag- 
netic latitude, 7=45°, figures 4,5. Both amplitude, 
|T|, and phase, arg 7, are presented, ¢,=0—360°, 
at frequencies 10 and 100 ke/s. Note that the com- 
ponent, 7 .., is symmetrical about the eastward, 
¢.=90°, and the westward, ¢,=270°, propagation. 
Note also that the direction of maximum intensity, 
¢o,=90° (eastward propagation) is considerably 
greater than the corresponding reverse direction, 
$,=270° (westward propagation). 

Propagation via the ionosphere at the magnetic 
equator, J=0, tables * 1, 5, 9, 18, 17, exhibit small 


3 The integer to the right of each table entry, if present, indicates the power of 
at aad ten (10) by which the entry is to be multiplied. Thus, 6.4307—1= 
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Ficure 4. Model ionosphere reflection coefficients, amplitude, 
|T|, and phase, arg T, illustrating a dependence of the re- 
flected field upon the direction of propagation relative to the 
direction of the earth’s magnetic field vector. 


N=1.2 (108), v=106, 6;=80.397°, Hm=0.5, I=45°, f=10 ke/s. 
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Figure 5. Model ionosphere reflection coefficients, amplitude, 
|T|, and phase, arg T, illustrating a dependence of the reflected 
field upon the direction of propagation relative to the direction 
of the earth’s magnetic field vector. 


N=1.2 (108), v=108, ¢:=80.397°, Hm=0.5, [=45°, f=100 ke/s. 
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and substantially constant abnormal components 
(|T'm| and |Tme| are between 0.01 and 0.04) in all 
directions, ¢,, for a large limiting value of the angle 
of incidence, ¢;=82°, and a large collision frequency, 
v=2 (10’). It is interesting to note that the com- 
ponents differ in phase by z radians; i.e., arg T.m—arg 
Tmc= +p, where p is an odd integer, p=1, 3, 
5... .. The abnormal components, Tem and Tne, 
appear to interchange values abruptly at ¢,=90° 
and ¢,=270° at the magnetic equator. The abrupt- 
ness of this interchange of values decreases as the 
magnetic inclination, J, increases. It is of interest 
to note, tables 9, 13, for example, that north-south 
propagation, ¢,=180°, and south-north propagation 
¢a=0, do not exhibit a precise reciprocity at the 
magnetic equator, J=0. Note (tables 9, 13): 


ga=0 180° 


| Teel 

arn T., 
=|T eal, 

arg Tem =arg 
| =| 

ore Tins 

| Tne 

arg =arg 


Thus, at the “magnetic equator’, 7=0, the only 
nonreciprocity in the propagation along the mag- 
netic meridian is as might be expected, a phase shift 
of x in the abnormal components as each abnormal 
component vanishes at precisely eastward or west- 
ward propagation. 

Propagation from the model ionosphere at the 
magnetic north pole, J=90°, tables 4, 8, 12, 16, as 
might be expected, is independent of the magnetic 
azimuth, ¢,. A value, ¢,.=0, was computed at all 
frequencies, f=10-22 ke/s and a value, ¢,=240° at 


f=10 ke/s served as a check on the computation. 


The behavior of the propagation in the vicinity of 
the magnetic north pole, /=84.270° is illustrated, 
tables 3, 7, 11, 15. 

Propagation from the model ionosphere in “tem- 
perate’’ magnetic latitudes, J=45°, tables 2, 6, 10, 
14, figures 4, 5, exhibit abnormal components, 
Tem, Tme, Which in general differ not only along the 
magnetic meridian, but also at the various values of 
magnetic azimuth, ¢, (0 to 360°). Although the 
normal components, 7, Timm, in the north-south, 
¢,=180°, and the south-north, ¢,=0, direction 
(i.e., along a magnetic meridian) are identical, the 
propagation does not exhibit a precise reciprocity 
since the abnormal components, Tem, Tme, are not 
identical in these directions. However, reciprocity 
for propagation along a magnetic meridian can be 
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considered to be an approximately valid concept if 
the collision frequency, v, and angle of incidence, 
¢:, are sufficiently great such that the abnormal 
components, T's, Z'ne, become quite small. Indeed, 
the large values of the angle of incidence, $;=82°, 
employed in the calculations for the tables tended, 
in conjunction with the high collision frequency, 
v=2 (10’), to reduce greatly the abnormal compo- 
nents. 

The intensity of the abnormal component propa- 
gated from the ionosphere excited by a vertically 
polarized transmitter, |7'm|, has special interest 
at LF, f=100 ke/s, figure 5. The eastward propaga- 
tion, $,~90°, from a vertically polarized transmitter 

rovides a ‘dominant vertically polarized field, 
Te|, at the receiver; whereas, westward propaga- 
tion, ¢,~270°, from a vertically polarized trans- 
mitter can under certain conditions (see fig. 5 
specifications) provide a dominant horizontally 
polarized field, | 7'| (abnormal component) at the 
receiver. In practice, it is necessary to consider 
the influence of the displacement and conduction 
currents in the ground which can indeed discriminate 
between horizontal and vertical polarization in such 
a manner as to obscure the ionosphere phenomena 
described. 


4. Conclusions 


The propagation of low- and _ very-low-radio- 


frequency waves reflected by the ionosphere is 
indeed dependent upon the ‘Lorentz force’’ on the 
electron, and hence is dependent upon the direction 
of propagation relative to the direction of the 
earth’s magnetic field vector in the region of the 


ionosphere in which the reflection occurs. The 
intensity of the abnormal components is depend- 
ent upon the values of the various parameters 
(such as N, v, Hn, o2, J, 6) and in particular, the 
collision frequency, and the angle of incidence 
(distance) ¢; are especially important in the determi- 
nation of the precise value of the abnormal compo- 
nents, Precise reciprocity in the LF and VLF 
propagation (interchange of transmitter and re- 
ceiver) can occur only at the magnetic north or 
south pole, J=+90°. Otherwise, the propagation 
is dependent upon the direction, ¢,., and in general 
the intensities become greatest in an eastward 
direction, ¢,~90°. Furthermore, it is conceivable 
at low frequencies, if the influence on the ground 
can be neglected or separated from this consideration, 
that the horizontally polarized component excited 
by a vertically polarized transmitter can approach 
in magnitude and even exceed the vertically polarized 
component in the case of westward propagation, 
whereas it would appear that the corresponding east- 
ward propagation would exhibit a dominant vertical 
polarization. 
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Summary of tables 


Table | Polar- N y I da f 
ization 
el/em’ c/s | Degrees| Gauss | Degrees| Degrees; kec/s 
: e 1038 2(107) 82 0.5 a 0-300 | 10-22 
3 84. 27 
4 90 (2) 
5 m 0 0-300 
6 45 
84. 27 
8 90 (2) 
9 e 3(108) 0 0-300 
10 45 
ll 84. 27 
12 90 (2) 
13 ™ 0 0-300 
14 45 
15 84. 27 
16 90 (2) 
17 em 108 0-108 0 0 
® Independent of ¢a. 
565004—60-——6 65 
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Correlation of Monthly Median Transmission Loss and 
Refractive Index Profile Characteristics ' 


B. R. Bean and B. A. Cahoon 
(July 14, 1960; revised August 5, 1960) 


The difference in the monthly mean values of the refractive index at ground level and 
at one kilometer above the ground level is often used for the purpose of predicting the annual 
cycle of radio transmission loss. The present study investigates the possibility of utilizing 
differences to heights other than one kilometer. A comparison of 100 megacycles per second 
transmission ioss recorded over twenty-one paths with various refractivity differences from 
the surface to three kilometers reveals that the surface value of the refractive index yields as 
good a correlation as any of the refractive index differences due to the high correlation be- 
tween the surface values and these differences; therefore the more accessible surface value 
can be effectively substituted for the differences. Specifying the refractive index profile at 
two or three additive points and using multiple correlation techniques does not significantly 
increase the correlation. 

The use of radio data over the same path does not significantly improve the correlation 
over that obtained from only meteorological data, indicating the very practical result that 
inexpensive meteorological data may be used to predict the seasonal trend of VHF radio field 
strengths with as much accuracy as expensive radio path measurements. 


1. Introduction Taste | 
Current engineering practice in estimating the level | Path Callletters' Transmitter | 
of beyond-the-horizon VHF-UHF radio field strengths | 
{1 to 5] * consists of determining an average level for 1|KXyYz | Houston, Tex....- Austin, Tex._.| “96.5 | 237.9 | 20.32 
a given distance and then adjusting this average KIXL | Dallas Tex Austin’ 1048 | 2881 | 32.38 
level for climatic and seasonal differences by reference 4| KWKH | Shreveport, La....| Austin, Tex...| 945 | 446.4 | 47.36 
view, Tex... Austin, Tex__. 5. 5.6 | 40. 
to common meteorological measurements of pressure, 
“atur iditv , adi 6| WFAA | Dallas, Tex.......| Austin, Tex...| 97.9 | 280.3 | 32. 54 
temperature, and humidity as expressed by the radio | | | 
refractive index, n. Although past investigations 8| WCsI_ | Columbus, Ind..--| Urbana, Tl----| 93.7 | 223.7 | 20.88 
h llv tili d tl diff tw 9 KXOK St. Louis, Urbana, IIl_.__| 93.7 235.8 | 19.34 
ave principally utilized the ailference in n between 10| WHKC | Columbus, Ohio--| Hudson, Ohio_| 98.7 | 201.0 | 19.42 
the earth’s surface and 1 km above the earth’s surface, 
better results might be obtained by using differences 12] WJAS | Pittsburgh, Pa....| State College, | 99.7 | 188.5 | 30. 44 
4 a. 
to heights other than 1 km or between different levels 13| WEST | Easton, Pa......- State College, | 107.9 | 224.2 | 36.67 
on the profile. It is the purpose of this study to 14| WTOP | Washington, D.C_| State College, | 96.3 | 215.8 | 63.47 
explore this possibility via the correlation of the | 
annual cycles of monthly median radio transmission ' en ne | 
Joss and n gradients obtained over a variety of height 16| KFOR | Lincoln, Nebr.....| Grand Island, | 102.9 | 150.0 | 14.50 
increments. The sections that follow deal with the sc. | 
available radio and meteorological data, deriving | 
-orr. i iscussi ir 18| KDKA Pittsburgh, Pa....| Hudson, Ohio} 92.9 | 155.6 | 12.15 
correlations between the two and discussing their | Hartford: Millis | 129-9 | 16.29 
significance on a regional and a statistically-combined 2 | NBS-_ | Cheyenne Mt. Haswell, "| 92.0 | 158.8 | 6.400 
4 010. 
basis. 21 | WIP-FM | Philadelphia, Pa-.| Laurel, Md...) 93.3 | 167.7 | 12.04 
2. Radio and Meteorological Data 
and 21) the Gulf States (paths 1 to 6), Great Lakes 
The radio data used in the present studies are the | (paths 4 
the 21 transmission paths listed in table 1. All es 0 
data correspond to the hours of the 0300 and 1500 y It p 
G.m.t. radiosonde observations (1000 and 2200 = given on. fr 
the central and southeastern seaboard (paths 17 | 108s DeIng tow in the summ S and hig 
the winter months. It appears to be quite significant 
that similar annual cycles are observed in such differ- 
Contribution from Central Radio Propagation Laboratory, National Bu- | ent climatic regions as the Appalachian Mountains 
of Pennsylvania and southern California coast. 


? Figures in brackets indicate the literature references at the end of this paper. 
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Figure 1. Annual cycles of monthly median basic transmission 
loss for four geographically diverse radio paths. 


These annual cycles plus the data indicated by 
table 1 were used as a test of various radio meteor- 
ological parameters that might be used in propaga- 
tion studies. The refractive index gradients were 
prepared from 6-yr means of temperature, pressure, 
and humidity as a function of height by the National 
Weather Records Center of the U.S. Weather 
Bureau. These data were then converted to 


N=(n—1)10° (1) 


by means of the Smith-Weintraub relationship [6]. 
Simple differences of the value of N at the earth’s 
surface, N,, and its value at 0.1, 0.2, 0.3, 0.5, 0.7, 
1.0, 2.0, and 3.0 km above the earth’s surface were 
obtained. This procedure is illustrated for an ex- 
ponential curve on figure 2. For example, the value 
of the difference to 500 m, AN) 5, is determined by 


ANo.s=N.— (2) 


The variations of transmission loss were also compared 
with those of N,; since it has received such wide 
engineering applications [1 to 4, 7 to 9]. 


3. Choice of Statistical Model 


When analyzing two variables for a possible inter- 
relation one has available many statistical techniques 
such as linear, multiple, and polynomial regression. 
Individual graphs of transmission loss and AN; were 
prepared for each of the gradient levels on each of 
the propagation paths to aid in the choice of an 
appropriate statistical model. Visual investigation 
of these data, examples of which are shown on figures 
3 and 4, reveals the presence of a large random 
component that masks any obvious nonlinear relation 
between the variables. On this basis it was decided 
that use of nonlinear regression techniques was an 
unnecessary complication at this stage of the analysis, 
and the simple linear regression model was adopted. 
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Fiaure 2. Graphical representation of the determination of 
various N differences in an average atmosphere. 


4. Statistical Considerations 


We are here concerned with the degree to which 
the annual cycle of radio transmission loss may be 
predicted from the several AN, values described 
above under the assumption of a linear regression 
model. This is done on a quantitative basis by 
determining the linear correlation coefficient, 7ry.:, 
between a dependent variable, y, (in this case the 
transmission loss) and an independent variable z 
(whichever AN, we are considering). The reason for 
determining the correlation coefficient can be under- 
stood more easily in terms of the standard error of 
estimate, (SE),, which is given by 


2)", (3) 


where s, is the standard deviation of the dependent 
variable. Equation (3) allows one to evaluate quan- 
titatively the degree to which one may predict the 
variations of y by knowing its correlation with z (and 
the concomitant regression equation). Thus it is 
seen that if y is completely determined from 2, and 
ry-2=1, then (SE),=0. Conversely, if x and y are 
totally independent, then 7,..=0, and (SE),=s,, and 
there is no advantage gained by attributing variations 
of ytoz. We may define the improvement factor 


sy— (SE), 
k=100| (4) 


Fy 


TaBLE 2. Coefficient of correlation between monthly median transmission loss and indicated 


meteorological variable. 
Path of N, ANo.2 ANo.3 ANos ANo.1 ANi.0 AN?2.0 AN3.0 
ata 
1 38 0. 586 0. 504 0. 517 0. 525 0. 553 0. 589 0. 612 0. 626 0. 624 
2 36 618 . 462 . 527 552 . 593 . 614 . 632 . 659 . 651 
3 68 787 . 670 . 690 . 707 . 724 737 . 768 787 786 
4 41 774 . 619 . 623 . 640 . 654 657 . 706 714 7 
5 23 67 541 536 . 539 589 617 627 
6 7 . 940 . 806 . 850 . 851 . 883 . 890 . 909 916 
7 22 . 714 . 539 559 . 567 645 695 718 
8 48 . 726 713 . 729 . 729 . 742 . 763 . 766 767 745 
9 24 . 588 626 . 643 . 639 . 627 . 638 . 626 627 608 
10 74 - 506 527 . 540 540 . 550 . 540 556 518 
ll 64 . 525 . 542 . 552 . 554 575 539 
12 41 . 868 778 . 858 871 . 874 . 866 873 874 869 
13 26 - 785 590 . 724 734 725 . 743 778 800 
14 12 935 . 733 . 887 886 . 893 901 . 928 940 938 
15 28 . 797 —. 128 —. 228 —. 246 —. 134 . 731 . 802 787 | 814 
16 17 626 . 542 576 619 . 654 660 666 655 | 
17 17 . 373 . 247 314 . 370 384 466 455 421 
18 30 . 803 . 809 830 . 836 . 832 797 791 . 787 799 
19 13 oe 214 686 .779 782 775 760 . 778 768 
20 13 115 . 364 373 . 299 270 214 118 117 
21 38 652 579 806 . 781 761 722 688 655 7 
Mean values (omitting path 20) 
r(z) 667 0. 694 . 578 0. 633 0. 640 0. 655 0. 676 0. 702 0. 704 0. 700 
k[r@)] 667 28.0% 8.4% | 22.6% | 23.2% | 24.4% | 26.4% | 28.8% | 29.0% | 28.6% 
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Figure 3. Monthly median basic transmission loss versus Ns or 
AN, for Path 12 (Pittsburgh to State College, Pa). 


as the percentage reduction of uncertainty of pre- 
dicting y from knowledge of its correlation with z. 

It is frequently desirable to estimate the average 
value of the correlation coefficient when dealing with 
several correlation coefficients between the same 
variables but derived in different climatic regions. 
Since 7 is not normally distributed it is desirable to 
convert 7 to Fisher’s 2, 


z=}[log, (1+r)—log, (l—r)] 


whien is more nearly normally distributed, obtain 
the average of z, and convert this average 2 to an 
estimate of 7 [10]. 


5. Comparison of Annual Cycles 


Individual correlation coefficients of the monthly 
median transmission loss and the mean value of the 
various refractive index gradients were determined 
for each path and are listed in table 2. Figure 5 
depicts the relative behavior of 7 versus N difference 
to various heights for the same four radio paths as 
those graphed in figure 1. It can be noted from 
figure 5 and table 2 that the magnitude of the cor- 
relation coefficients is generally lowest when deter- 
mined versus ANp.; and rises with increasing Ah until 
the maximum value is reached at 1.0 or 2.0 km. 

Table 2 also lists the average correlation coefficients 
[10], indicating that the correlation of radio trans- 
mission loss with N, is as large as those obtained with 
any value of AN,. 

There are some notable exceptions to the above 
observations. For example, path 18 (Pennsylvania 
to Ohio) and path 9 (Missouri to Illinois) displayed a 
fairly uniform correlation for all height increments 
with a very small range of coefficients. Path 19, 
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where the 0.5 km level evinced the highest correla- 
tion, and path 21, where the maximum correlation 
appeared at 0.2 km, also had a very low range with 
the exception of the 0.1-km level where the cor- 
relation coefficient was significantly lower. 

The low correlations for path 20 might arise from 
the fact that this path is a borderline knife-edge 
diffraction case, as indicated by the low value of the 
angular distance, 0. 

Perhaps the most interesting case is that of path 
15, San Diego to Santa Ana, California, where the 
largest correlations are found between the trans- 
mission loss and N, and gradients for height incre- 
ments in excess of 0.7 km. The correlation co- 
efficients for height increments below 0.7 km are not 
only relatively small but are negative. This physi- 
cally unrealistic behavior, shown in more detail in 
figures 4 and 5, may be explained by the presence of the 
well-known southern California coastal inversion 
which occurs at about 0.7 km and is known to exert a 
dominating influence upon radio field strengths in 
that region [11, 12]. The interesting conclusion here 
is that the variations of radio fields follow those of 
N, and any of the N differences to heights above the 
base of the inversion but do not follow the inter- 
mediate differences. 
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6. Other Measures of the Gradient 


The preceding analysis has been concerned with 
the correlation of transmission loss and AN, when h 
takes on several different values. It is also profitable 
to examine the correlation of transmission loss and 
differences in index taken between various points on 
the N profile such as has been found useful in Ger- 


COLUMBUS-URBANA WTIC HARTFORD-MILLIS 
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08 | WALUES FOR Ns 4 
04 | LIMITS 
02+ | 
0 | 
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| VALUES FOR Ng 
Z 


90% CONFIDENCE 
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90% CONFIDENCE 
LIMITS 


CORRELATION COEFFICIENT 


0 1000 2000 300-0 1000 2000 3000 
HEIGHT IN METERS HEIGHT IN METERS 
Figure 5. Correlation coefficients of monthly median basic 
transmission loss versus height in meters of AN}, for four geo- 
graphically diverse paths. 


many by Grosskopf [13]. Another profile measure is 
Misme’s “useful gradient” [14], G, given by 


_2A4N,+AN, 


where AN» is the initial gradient and AN, is the 
gradient at the crossover height of the radio horizon 
rays. The correlation coefficients resulting from 
such a comparison are tabulated in table 3 and, it is 
seen that the average correlation coefficients are 
generally lower than those obtained previously. 
The general conclusion one reaches here, then, is that 
none of these other profile measures offers a significant, 
improvement over the widely used N, and AN, .0. 


G 


7. Use of Previous or Concurrent Radio Field 
Strengths to Predict the Annual Cycle of 
Radio Fields 


A perspective on the magnitude and usefulness of 
the above correlation coefficients may be gained by 
comparing the correlation of transmission loss 
recorded over nearly identical paths. Seven of the 
radio paths on table 1 had sufficient data so that the 


3In the course of a private conversation, P. Misme has pointed out to the 
authors that recent studies have led him to adopt a different value for the altitude 
at which the parameter dN/dh (referred to as AN, in our equation for @) is meas- 
ured. Furthermore, calculation by layers is in every case preferable as a method 
of obtaining the value of the “equivalent gradient.’’ Unfortunately, this com- 
ment reached the authors too late to be considered in this work. However, a 
forthcoming article by the authors will compare the variation of the equivalent 
«sorting snag other measures of atmospheric refraction with variations of trans- 
mission loss. 
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TABLE 3. Coefficient of correlation between monthly median 
transmission loss and indicated meteorological variable 


Path | G | ANo.2-0.5* | AN?.5-1.0 AN?.0-3.0 N; 
1 0.5477) 0.04999 —0. 2892 0.4042 | 5082 0. 586 
2 | —.1219 —.2191 2297 .5519 | 
3 6959 4511 2437 | . 6747 | 7425 | 787 
4 | | 3679 1349 5567 5925 774 
5 4166 1988 | 5455 672 
6 | 8383 5087 . 6355 . 7893 | 8666 . 940 
7 | 5708 1986 6610 | 714 
8 7351 6248 5134 6105 8102 
9 | 6583 ‘5700 6980 | 588 
10 | 5312 4951 4531 5694 5681 | 
11 | 5451 4933 4207 5743 | | 525 
12 8178 8 5985 | 7920 868 
13 6746 3432 3610 5836 6850 785 
14 —. 2005 8162 8775 935 
15 —.1262| —.3484 —. 4062 7793 1 797 
16 . 5636 3844 . 5682 . 6797 6889 626 
17 | =—.0350 3981 | 6309 373 
18 8213 7762 7643 6950 | 6336 803 
19 | 5428 | 8087 7599 566 7712 744 
20 | 3636 | 13519 1786 1188 115 
21 | | —. 09026 8065 6566 652 
Mean values (omitting path 20) 
| 
r(2)_ 0. 6224 | 0. 3548 0.3811 0.5461 | 0. 6624 0. 694 
kir(2)] 21.7% | 6.5% 7.5% 16.2% 25.1% 28. 0% 
| | | 


*Note: No..-0.5 indicates that the transmission loss was correlated with the 
difference in N taken at h=0.2 km and h=0.5 km. 


first 12 months of record could be used to predict the 
monthly median radio field strengths in subsequent 
years. The results of this study, including the mean 
correlation coefficient and the improvement factor, 
are given in category 3 of table 4. A comparison of 
the improvement factor with those of tables 2 and 3 
leads one to the conclusion that for any practice] 
prediction problem, where only past radio or meteor- 
ological data are available, one obtains at least as 
ood a prediction of the annual cycle of transmission 
oss from long-term meteorological data as from 
relatively short-term radio data. 

This conclusion arises from the fact that the cor- 
relation between paths 1 and 2 is low (0.348) even 
though the paths are so nearly identical that one 
would not hesitate to make such a comparison for any 


TaBLE 4. Correlation of annual cycles of transmission loss 
over identical or nearly identical paths 


Category! Paths Data 90% confi- 
dence limits 


oe: . See Successive 15-month peri- 30 | 0.348 | 0.295 to 0.399. 
ods. 
2 | 10&11__._| Coincident months. 58 . 838 | 0.829 to 0.847. 
a First 12 months of record 44 . 767 | 0.750 to 0.783. 
4. on subsequent years of 17 . 762 | 0.708 to 0.807. 
8. record. 24 a 0.730 to 0.795. 
3 |(9. 12 . 624 | 0.500 to 0.723. 
12. 18 . 737 | 0.683 to 0.783. 
15. 16 . 656 | 0.578 to 0.723. 
21, 14 . 890 | 0.859 to 0.917. 


=0. 751 
k[r(z)]=34.0% 


*j=pieces of data. 


engineering applications. The use of radio data 
from the same path (category 3) does yield a higher 
correlation than the use of meteorological data alone; 
but this type of prediction, running the actual circuit 
for a year, results in only a slight improvement in 
prediction at a cost increase of many orders of magni- 
tude over that of procuring the necessary meteoro- 
logical data. However, any actual radio path 
strength testing would probably fall into category 1 
where data from nearly identical paths are used, and 
the resulting correlation appears to be less than that 
obtained from meteorological data alone. Category 
2, use of simultaneous radio data, yields the highest 
correlation but is generally an impractical procedure 
for engineering applications. 

These results have a very important implication 
for future radio meteorological studies involving 
monthly median field strengths. It is hard to imagine 
any predictor of field strength variations that would 
be more effective than a nearby radio station at 
nearly the same frequency, recorded over nearly 
identical paths to the same receiver site. In this 
case it would seem that the maximum reduction in 
the uncertainty of predicting fields from monthly 
mean refractive index data could not exceed about 
35 percent. If this is indeed true, then further work 
along these lines would seem quite futile, and new 
approaches should be tried. Perhaps a profitable 
new approach would be to combine the effects of 
day-to-day variations of refractive index profiles in 
some manner similar to that suggested in the earlier 
works of the senior author [15]. 


8. Multiple Regression 


One might wonder if a combination of the above 
meteorological variables would produce a greater 
reduction of the uncertainty of prediction. Several 
combinations of meteorological parameters that 
might occur in practice were evaluated. These were 
transmission loss as a function of both N, and AN, y 
(assuming that both N, and AN could be obtained 
from climatic maps for the United States [1, 16] or 
the worldwide climatic maps in preparation by the 
International Consultative Committee for Radio) 
and also as a function of N, and AN, ., such as might 
be obtained from meteorological towers. These 
data, reported in table 5, indicate that very little 
improvement may be accomplished by using two 
meteorological variables determined from radiosonde 
data. Indeed, even when transmission loss as a 
function of all three of these meteorological para- 
meters is considered, the augmentation of the correla- 
tion coefficient does not appear to be significant. 
In fact, the examination of the linear relationship 
involving all three of the meteorological factors, 


ota 
ota 
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TaBLE 5.—Correlation coefficients 
of transmission loss versus multi- 
ple meteorological factors 


Path Ty.2y25" Ty.2,23" Ty.2y2923" 
1 0. 586 0. 613 0. 622 
2 636 633 . 669 
3 . 797 . 787 . 800 
4 . 794 . 812 . 814 
5 . 696 . 790 . 790 
6 - 958 . 958 959 
7 . 758 . 716 782 
8 . 733 . 818 822 
9 . 627 . 688 7 
10 . 528 . 607 610 
ll 542 626 628 
12 . 869 . 874 875 
13 . 787 . 785 787 
14 - 936 - 936 936 
15 - 823 . 816 837 
16 690 . 730 798 
17 435 . 592 7 
18 . 812 803 813 
19 . 749 7 77 
20 - 500 315 | 
21 718 | 77 
| 
Mean values (omitting path 20) 
r°(Z) 0. 717 0. 742 0. 758 
k[r@)] 30. 3% 33.0% 34.8% 


*For this table y=L», 1=Ns, 12=AN0.1, 

discloses that the attempt to specify the refractive 
profile by three points apparently does not lead to an 
appreciable increase in the correlation coefficient. 
However, the error analysis carried out in the next 
section indicates that the probable errors associated 
with radiosonde measurements are sufficient to mask 
the true correlation value to the extent that this 
conclusion might well be modified by the use of more 
accurate meteorological data. 

The above correlations do reflect systematic prop- 
erties of the atmosphere which in themselves present 
an interesting study. It is not too surprising to note 
that the mean correlation coefficient between N, 
and AN, ois high (0.9840). This interdependence of 
N, and AN, ois fortunate since it allows the computa- 
tion of AN,» from surface data, which is readily 
available throughout the world. 


9. Error Analysis 


It is really quite surprising that the initial gradient, 
ANo.;, did not yield larger correlation coefficients 
since the refraction of radio waves is known to be 
extremely sensitive to the initial gradient. It is 
quite possible that this low correlation reflects the 
difficulty of determining an N difference over a small 
increment with relatively inaccurate radiosonde 
sensing elements. We may analyze this error bv 
reference to common meteorological standard at- 
mospheres and standards of measurement set by the 
International Consultative Committee for Radio 
(CCIR). By setting 


(Ni te1)—(N2+e2) 


AN,-:= (hi —ho) (6) 
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where N, and N, represent the true value of N at h? 
and and and represent the errors determining 
N at h, and h, respectively, then the standard devia- 
tion of the error in AN,_, is 


+07 (€2) 
a(AN;-2) (7) 


assuming that f4i—h, is known exactly and that « 
and ¢€ are uncorrelated. The errors « and ¢ may 
be determined by noting that [6] 


N=K, pt (8) 
and that 
oN oN oN 
de+sp dP. (9) 


The partial derivatives may be evaluated from the 
ICAO standard atmosphere [17] assuming 60 per- 
cent relative humidity to yield approximate rela- 
tionships 


AN=—1.27AT+4.504e+ 0.274P (10) 
for h=0 and 
(11) 


for h=1 km. 

By assuming the relatively stringent accuracies 
in measuring P, T, and e called for by the CCIR 
[18] one obtains standard errors of 3.5 percent of the 
mean in determining AN;o, and 22.7 percent of the 
mean in determining ANo,1. 

Reference to standard treatment of the effect of 
measurement errors upon correlation coefficients [19] 
shows that if y is the dependent variable and the 
indpendent variable s*=zx-+e where is the true 
value and e¢ is the measurement error then the 
correlation coefficient is 


* Py.z (12) 


Vi+ 


assuming that « has zero mean and variance s*(e) 
and is uncorrelated with either z or y, the other 
terms having their usual meaning. It is seen that 
the effect of measurement error is to reduce the 
correlation coefficient by the factor 


1 


v1+s?(e)/s*() 

If we now assume that y (in this case the transmis- 
sion loss) is correlated with two different measures 
of the same variable, x; and x, (in this case ANo. 
and ANj») then from the ratio 


1 


| 
( 
5 
| 
| 


y 


one may determine the relative size of the measure- 
ment errors necessary to yield the same correlation 
coefficient by setting 


(15) 
Py.zg 


and obtaining 
8(€:)/8(€) (a2). (16) 


Returning to our data we find that this ratio aver- 
ages about 1.5 for the observed values of s(4No.)/ 
s(ANio) while the ratio of the errors determined 
from the ICAO standard atmosphere and the CCIR 
standards of measurement is about 10.0. Thus it 
is seen that the errors of measurement may indeed 
be the cause of the differences in these correlation 
coefficients. Although one cannot estimate the ac- 
tual size of the true correlation coefficients, one can 
estimate their ratio by reexpressing (14) as 


When we again use the standard-atmosphere-esti- 
mated errors and the actual values of s?(x,;) and 
#(22) for the 21 paths in our analysis, we find that 
the ratio py.2,/py.z, averages about 1.5 which implies, 
since TLy.ANn, >~09.7, the actual correlation with the 
surface layer TLy.AN), 18 Near unity, an observation 


that should guide future experimental research in 
the quest of more efficient radio-meteorological 
parameters. 

Although it is evident that the errors of measure- 
ment do make the choice of a radio climatological 
parameter based upon this study somewhat indefi- 
nite, we may use the present study as an objective 
criterion of the relative merits of such parameters 
when determined from radiosonde observations. 
Perhaps the advent of more accurate measurements 
will shed more light on the true role of the surface 
layer in radio meteorology. 


10. Extension to Other Regions 


Although this study encompasses only a small por- 
tion of the United States because of limited radio 
data, several of the more important climates of the 
world are represented by the 21 paths. Paths 1 to 6, 
14, 17, and 21 lie in the same humid subtropical 
region which covers parts of Argentina, China, and 
Yugoslavia; the climate of paths 7 to 13, 16, and 18 
is similar to that of Hungary, Romania, and parts 
of Japan; path 15 represents the Mediterranean re- 
gions and includes southern Australia, Africa, and 
central Chile; path 20 has meteorological conditions 
which resemble the middle latitude steppes of Canada 
and southern Russia while path 19 is similar to 
southern Sweden and central Russia. It is reason- 
able to assume that the general pattern of the mean 
profile of a climatic area would repeat itself, regard- 
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less of the country in which it was located. Such 
an assumption would permit the extension of any 
conclusions reached in this paper regarding the cor- 
relative aspects of various meteorological parameters 
versus transmission loss to other regions. In fact, 
studies which have been completed of N, versus 
transmission loss show good results not only in the 
Florida-Caribbean [3] area but also with trans- 
mission paths from Argentina to Uruguay [3], from 
Minorca to Sardinia in the Mediterranean [4], and 
from Tokyo to Osaka in Japan [20]. 


11. Conclusions 


In this study the calculation of the correlation 
coefficients of recorded transmission loss over 21 
radio paths versus various properties of the refractive 
index profile up to 3 km provides a means of evaluat- 
ing which radio-meteorological parameter is most 
efficient to predict average radio transmission loss. 
The best correlations appeared to be associated with 
the time-honored values of AN,» and N,, indicating 
that tbev continue to yield as good a correlation as 
any other parameter derived from currently available 
radiosonde observations. 

While comparing the interdependence of meteoro- 
logical parameters, the correlation of surface refrac- 
tivity, N,, with a gradient interval of 1 km was found 
to be high, explaining why the more accessible N, 
can be used with significant results in place of the 
1-km gradient in the prediction of transmission loss. 

An examination of the correlations of transmission 
loss recorded for 1 year versus that recorded in 
another year for the same radio path revealed corre- 
lations of the same order as those obtained from 
meteorological data alone, indicating the very prac- 
tical result that inexpensive meteorological data may 
be used to predict the seasonal variation of radio 
fields as accurately as expensive radio path meas- 
urements. 


The authors acknowledge the assistance of K. A. 
Norton, G. Richmond, B. J. Weddle, and P. C. 
Whittaker in this study. 
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Characteristics of Waveguides for Long-Distance 
Transmission’ 


A. E. Karbowiak and L. Solymar 


(May 23, 1960) 


Discussion of a waveguide communication system is given and the system is compared 


with existing communication media. The principal properties of waveguide 


the medium 


of communication—are discussed in some detail. The particular significance of helical 
and coated waveguides is pointed out and the design formulas included. The phenomenon 
of mode conversion-reconversion, which is peculiar to a waveguide communication system, 
is discussed in general and the basic theory as applicable to design is also discussed. Design 
features of components such as bends, transducers, tapers, ete., are analyzed. The effect 
of waveguide discontinuities is analyzed in some detail and various aspects of signal distor- 


tion are also considered. 


1. Introduction 


Man, in his constant efforts towards progress, 
finds it necessary to establish better, quicker, and 
more numerous communication systems. It is 
even said that complexity of communication systems 
grows in proportion to the development of civiliza- 
tion. As the conventional means of communication 
become gradually used to capacity and the traffic 
continues to increase, one, naturally, turns towards 
more exotic means of communication. Waveguide 
is a possible winner [1, 2, 3] ? and in this connection 
special acknowledgment should be made of the very 
important pioneer and continuing work of the Bell 
Laboratories in this field. Although the final assess- 
ment can only be made against the background 
of economy, a detailed performance assessment 
must be the first step. For a number of rea- 
sons, as we shall see subsequently, such a wave- 
guide is essentially a circular tube about 1 to 3 in. 
in diameter and the waves involved are in the 
millimetric band. Such a pipe—the medium of 
communications—would be laid in the ground or 
at the bottom of oceans for hundreds and thousands 
of miles. 

In order to bring out any advantages that a wave- 
guide communication system may have over the 
existing systems, the principal features are shown 
tabulated in table 1. From this it is evident that 
the main advantage of a waveguide communication 
system is its large communication capacity. Pro- 
vided that very short millimetric waves are used 
then, in principle, several hundreds of television 
channels, or the equivalent in speech channels, can 
be accommodated in a single pipe. The waveguide 
is also characterized by small attenuation, and in 
contrast to microwave links it is a screened system. 
Whether such enormous communication capacities 
will ever be required is another matter; but, if 
present statistics are anything to go by, then these 
capacities will be required in less than 20 years time, 


‘ Contribution from Standard Telecommunication Laboratories, Ltd., London 
Road, Harlow, Essex, E 


ngland. 
? Figures in brackets indicate the literature references at the this paper 
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if not for telephone communication then for pur- 
poses such as data transmission, color and cinema 
television as well as for purposes not yet conceived. 

A waveguide communication system is essentially 
simple (fig. 1). The video channels are grouped 
together (say by frequency division) and passed to 
a coder (C) where, for reasons to be discussed, it is 
coded in PCM (pulse code modulation). The coded 
signal is then made to activate an on-off modulator 
(M) operating on one of the many carriers in the 
whole spectrum carried by the waveguide. The 
carriers are combined through a suitable channel 
insertion filter and passed through a transducer (7) 
on to the waveguide. At convenient intervals along 
the waveguide (say every 10 miles) regenerative 
repeaters are inserted to amplify and restore the 
quality of the signal. 

For effective operation the terminal equipment 
(modulation/demodulation and coding/decoding) is 
designed to handle millimicrosecond pulses. This in 
itself is a complication but it is a necessary one to 
achieve PCM coding of the wide video band and 


["h C — CODER 
M— MODULATOR 
| P — SOURCE OF MICROWAVE POWER 
T — TRANSDUCER 
R — REPEATER 
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CHANNELS 
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Figure 1, A possible communication system using waveguide, 
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TABLE 1, 


Comparison of communication media 


Norte: One super channel=4 Me/s, 


Repeater spacing 


Type Operating frequency Channel capacity attenuation 
One pair 36 to 84 kc/s (go) 12+ 3+ (1) = 15 + (1) 
92 to 140 kc/s (return) 
Open Maximum As above 240 + (16) (0.4 db/mile) 
wire 16 pairs 
Balanced 
pair One pair 12 to 252 kc/s 60 
12 to 30 miles 
Maximum As above 1440 (4 db/mile) 
24 pairs 
a Up to 4 Mc/s 960 6 miles 
Coaxial (One super channel) 
G.B. 
Oe Up to 12 Mc/s 3 super channels 3 miles 
Up to about 3 Mc/s 720 8 miles 
Coaxial 
U.S. 
1 TV channel , 
Single wire About Probably a 


transmission line 


100 to 1,000 Mc/s 


few TV channels 


About 2 to 10 miles 


Microwave links 


Below 500 Mc/s 


Maximum of 60 


40 to 50 miles 


500 to 1,000 Mc/s 


Up to 120 


40 to 50 miles 


2,000 Mc/s 


240xX6 or 1TV X6 


30 to 40 miles 


4,000 Mc/s 


600xXo or 1TVX6 


25 to 30 miles 


6,000 to 8,000 Mc/s 


600X6 or 1TVX6 | 
Maximum 2 TV X 6 | 


25 to 30 miles 


11,000 Mc/s 


Less than 600 


Less than 20 miles 


Long haul waveguide 


(H,,) 


30,000 Mc/s 
up to about | 
100, 000 Mc/s | 


1,000 super channels 
= several hundreds 
of TV channels 
several hundred | 
thousands of 
speecii channels 


20 to 40 miles 
(attenuation 2 to 


4 db/mile) 
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derive the advantage of regeneration. Regeneration 
is a necessary requirement because of the large signal 
distortion, even with waveguides of moderate length 
(say 10 miles). It has been established beyond any 
doubt that for long-distance communication pur- 
poses only modulation systems capable of regenera- 
tion can be used, and of these PCM has definite 
technological advantages. 

This paper is not an attempt at a balanced sum- 
mary of the subject matter, but is a presentation of 
the opinions of the authors. 


2. Summary of Properties of Waveguides 


There are four principal factors influencing the 
erformance of a waveguide: (1) Its geometry both 

in shape and size (in terms of the wavelength), (2) 
its mode of operation, (3) the value of the surface 
impedance of its walls, and (4) the nature and 
magnitude of the tolerances on the cross section. 

For long-distance communication waveguides of 
circular cross section operating in the Ho, mode have 
been chosen. The choice is chiefly a matter of tech- 
nological preferences based on a number of consid- 
erations. Such waveguides have relatively low 
dispersion and low attenuation and can be con- 
structed to discriminate against unwanted modes, 
and the requirements on the tolerances on the cross 
section are less exacting than for other waveguides 
and modes. 

The exact waveguide properties and construction 
is a matter of compromise between a number of 
conflicting factors such as cost of terminal equip- 
ment and repeaters, permissible signal distortion, 
service reliability, etc. And, since these factors can 
be traded for each other it is impossible to give a 
specification for a waveguide under any given con- 
ditions except in the light of accumulated experi- 
mental experience. But, to fix orders of magnitude 
a Waveguide operating in the 4- to 8-mm band would 
be about 2 to 3 in. in diameter. 

With waveguides there are two principal sources 
of signal distortion: (1) Delay distortion [6], and 
(2) mode conversion-reconversion phenomena [1, 4]. 
Delay distortion arises (like with most communica- 
tion media) because of the nonlinearity of the phase 
characteristic, B(w). This leads to the undesirable 
feature that the high frequency components of the 
signal travel faster than the low frequency com- 
ponents, a phenomenon known as dispersion. <A 
corresponding distortion due to the variation in 
the attenuation characteristic a(w) is, with wave- 
guides, negligible [6]. 

The actual bandwidth limitation due to delay 
distortion depends to a large extent on the modula- 
tion method and is particularly troublesome with 
FM. PCM, on the other hand, is rather immune, 
but even here 10-musec modulation pulses would 
be unrecognizable after a journey down a typical 
waveguide, say, 30 miles long; regeneration is, 
therefore, called for. In addition, certain amount 
of equalization of the B(w) characteristic may 
be a necessity. 
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Mode conversion-reconversion is a new phenom- 
enon, peculiar to waveguide transmission systems. 
It arises because a typical waveguide suitable for 
long distance transmission can support a large 
number of modes apart from the desired one (H;), 
and because of the irregularities in the waveguide 
(small dents, changes in diameter and cross sectional 
shape, small bends and kinks, etc.). In consequence, 
at each minute irregularity the preferred trans- 
mission mode (Hg) becomes scattered and partially 
transformed into the parasite modes. The process 
is known as mode conversion and leads to a gradual 
diffusion of energy into the parasite modes and, 
therefore, increased attenuation. With some excep- 
tions this process of mode conversion is actually 
not particularly harmful. It cannot, however, exist 
on its own but, by reciprocity, it is accompanied 
by the complementary process of “reconversion.” 
In this process, the energy now partially carried in 
the unwanted modes becomes reconverted (by 
being scattered by subsequent irregularities) back 
to the Hy, mode. But since the group velocities of 
the parasite modes differ from that of the Ho, mode, 
the reconversion is not in phase with the signal 
carried in the Hy, mode. Further, since the degree 
of the distortion is phase dependent and the phase 
of the reconverted signal is related to the distance 
between the conversion and reconversion points in 
terms of the wavelength, the whole process of mode 
conversion-reconversion is frequency sensitive. The 
overall frequency attenuation curve is, therefore, 
very irregular, Jeading to considerable signal dis- 
tortion. Because of the nature of mode conversion- 
reconversion phenomena it is inconceivable that 
suitable terminal equipment could be constructed 
to counteract its effects. Mode conversion-recon- 
version effects are conveniently minimized by 
suitable waveguide design. 

There is one particular aspect of mode conversion- 
reconversion that calls for special attention: this 
arises through bending of the waveguide. The 
circular waveguide carrying the Hy, mode can 
support simultaneously the /;, mode, having identi- 
cal velocity.? In a straight waveguide this is of 
no direct consequence, but if the waveguide is bent 
in an are of a circle then the two modes become 
coupled with subsequent exchange of energy between 
the modes. This leads to large attenuation and 
signal distortion. Once again the effect is counter- 
acted by judicious waveguide design. 

It is possible to design a whole series of com- 
ponents [1, 5] such as specially tailored bends and 
corners, tapers, mode transducers, atc. But the 
design is not as straightforward as with components 
for conventional single mode rectangular waveguide. 
The waveguide can support a large number of 
modes and, therefore, any changes of the waveguide 
geometry must be made with great care bearing in 
mind mode conversion-reconversion phenomena. 
This requires a detailed study of wave transmission 


3 Such modes are termed degenerate. 


in multimode waveguides and the coupling effects 
due to irregularities; a field relatively unexploited 
in its analytical and experimental aspects. 

Despite all the planning difficulties, a waveguide 
when properly designed will give satisfactory service 
and will handle successfully a bandwidth well in 
excess of 20 kMe/s, which is equivalent to an intel- 
ligence bandwidth of about 1,000 Mc/s, a capacity 
not offered by any other existing communication 
medium. 


3. Delay Distortion in Uniform Waveguides 


As is well known, with waveguides both a(w) and 
B(w) are functions of frequency. It can be shown 
[6], however, that the effect of a(w) on signal dis- 
tortion is negligible in comparison with 6(w). Here, 
therefore, only the effect of the B(w) characteristic 
will be examined. 

The phase propagation coefficient B(w) is given by 


B= (w—«?)/c, (1) 


where c is the velocity of light, and w, the cutoff 
frequency of the waveguide, which is a geometrical 
parameter. These relations are illustrated in figure 2. 

The phase (v,,) and group velocities (v,) are 
defined by 


(2) 


These quantities are frequency dependent. 

If the signal is composed of several frequencies 
then, naturally, the signal will be distorted because 
the high frequency components will travel faster 
than the low frequency ones. The differential phase 
delay between the extreme components of the 
modulating envelope is a measure of signal distor- 
tion: it is the delay distortion term. For a wave- 
guide operated at a frequency substantially above 
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Figure 2, 6-w characteristic of a perfect waveguide, 


cutoff the bandwidth (in Me/s) can be expressed 
as follows [6]: 


f=(19/F) (fo (3) 
where 


F=cutoff frequency/carrier frequency 
l=length of waveguide in kilometers 
f.=carrier frequency in kilomegacycles per secona 
6¢= permissible phase delay between the carrier and 
the sideband in radians. 


It has been assumed that f is a small proportion of 
carrier frequency. 

For example, for a typical rectangular waveguide 
F=0.6, and therefore at the carrier frequency of 
10 kMe/s (and length 30 km) we get 

f=9.4 (1 rad)!?Me/s. 

The permissible amount of delay distortion (6¢) 
is a function of modulation method. Clearly, with 
AM if (6¢)=7/2 then inversion of sidebands takes 
place: a condition of severe distortion. For AM, 
FM, and many other modulation methods the maxi- 
mum possible value of 6¢ is only a small fraction of 
1 rad. For PCM, however, (é¢) can take quite 
substantial values; once again, therefore, PCM is 
superior to the more common modulation methods. 

For wideband application the modulation methods 
must be analyzed individually. For example, with 
pulse transmission systems, the energy contained in 
a pulse of finite duration gradually becomes spread 
out on the time scale. For a given waveguide 
length the degree of pulse distortion is a function of 
pulse duration. But up to a certain minimum 
pulse duration the pulse is quite recognizable and 
most of the energy is contained in the time interval 
equal to the duration of the transmitted pulse. 
For pulses of shorter duration the pulse distortion 
becomes quite severe and increases as pulsed dura- 
tion decreases [6]; simultaneously the energy be- 
comes diffused over a longer and longer time interval. 
The actual pulse width for which this happens is 
related to the geometry of the waveguide cross 
section and waveguide lengths (fig. 3). 


\4 : — 
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Ficure 3. Pulse distortion due to dispersion in waveguides, 
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Typical waveguides for long-distance transmis- 
sion have F=0.15, then at a frequency of 35 kMe/s 
pulses as short as 10 mysec can be transmitted 
without significant distortion over a distance of 
some 30 km. Such performance is sufficient for 
many applications, but if it were required to trans- 
mit even shorter pulses some equalization of wave- 
guide characteristic would then be necessary. Such 
equalization is best carried out at the terminal 
equipment. 


4. Properties of Uniform Waveguides and 
Their Uses 


4.1. Simple Waveguides 


The reasons for choosing for long-distance trans- 
mission the circular waveguide operated in the circu- 
lar electric mode (Hp,) have been explained in section 
2. Such waveguides can support, apart from the 
desired (Hp,) mode a multitude of other modes all of 
which are undesirable. The total number of modes 
that can possibly propagate in a waveguide with 
diameter D is known as the modulus of overmoding, 
M. This is approximately given by 


M=2.55(D/)2. (4) 


The larger the modulus of overmoding, the larger 
the number of possible modes and, consequently, the 
greater the engineering difficulties in preventing 
mode conversion-reconversion phenomena. Yet, a 
large ratio D/ is necessary to have low attenuation 
and dispersion. Of necessity, therefore, the propor- 
tions of a waveguide are a compromise. Typically, 
D/x is of the order of 10 or more and the modulus 
of overmoding is, therefore, at least about 200 to 400. 


Each of the modes propagates with a characteristic 
group velocity given by 


(5) 


The commonest type of waveguide is a plain 
metallic pipe. With such waveguides the wall of the 
waveguide presents a surface resistance FR, given 


by [7] 
R= V (afulo) V m (6) 


where o is the conductivity of the metal, u its permit- 
tivity, and ~/(u,/e,)=377 ohms, is the free space 
impedance. This surface resistance gives rise to 
attenuation which for E,,z7 modes is given by [7] 


De @) 
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and for H,,, modes by 


out) — 3:86 
mm D2 F?1—(m/x2)? 


R,. (8) 


These attenuations are given in decibels per meter if 
D, the diameter of the waveguide, is méasured in 
meters. 

The attenuation of the Hp, mode (the wanted mode) 
is given by [7] 


_ 8.686 F? 


(9) 


Thus for a copper waveguide at the free space wave- 
length of 8.7 mm, it is 1.2 db/mile. 

Most of the unwanted modes have attenuation 
and group velocities differing substantially from the 
Hy, mode. Thus attenuations of 10 to 50 times that 
of the Hj, mode is common and group velocities of 
most modes differ from the Hy, mode by as much as 
10 percent or more. A few of the modes, however, 
have properties similar to those of the Hy, mode. 
For example, the £,; mode has group velocity almost 
equal to that of the Hy, mode and a few of the modes, 
notably the An, Ags, Ay», Ai, Ay, and have 
attenuations only 3 to 10 times as much as the Hy, 
mode. These modes, therefore, tend to be particu- 
larly troublesome and the only way to deal with 
them satisfactorily is through judicious use of special 
waveguides. 


4.2. Special Waveguides 


Waveguides of special construction are necessary, 
as explained in the earlier sections, principally to 
combat or at least to minimize the mode conversion- 
reconversion effects. For such applications, disk 
waveguides were first to be suggested and various 
elliptic, corrugated, and more complex structures 
were tried for negotiating bends. After many years 
of research it is now clear that only helical, and 
dielectric coated waveguides are likely to be of any 
extensive use. Such waveguides behave as if their 
surface impedance were anisotropic. 

The theory of propagation in anisotropic wave- 
guides can be quite involved, but since waveguides 
for long distance transmission have, in general, small 
surface impedance an approximate treatment is 
adequate, and this will be given here [8, 9]. 

The propagation coefficient in such waveguides is 
given by 


Y¥=Yot 6y=JBo+ (a+ =a+ JB, (10) 


where the quantity 6y(=a+ 48) is composed of the 
attenuation coefficient (a) and a quantity 68 modi- 
fying the phase change coefficient (8,) of the wave- 
guide with zero surface impedance. The coefficient 
ais proportional to surface resistance (R,) and 68 to 
surface reactance (X,), these being respectively the 
real and imaginary parts of surface impedance, Z,. 
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With anisotropic waveguides the quantity Z, may 
be described in terms of the matrix of its components; 
most simply in terms of its two principal components 
Z, and Z; which are surface impedance components 
along the two principal axes. The sum of the pro- 
jection of these components along the ¢ and 2z direc- 
tions define the circumferential (Zs) and axial (Z,) 
components of the surface impedance; and it is these 
components that are used in the subsequent form- 
ulas. Z, can be regarded as the surface impedance 
presented to the currents flowing in the waveguide 
wall in the circumferential direction and Z, a cor- 
responding quantity for currents flowing in the axial 
direction. 

In such waveguides the attenuation of E-waves is 
given by [9] 


—2R,/Dy(1—F?) nepers (11) 


and for H-waves by 


| nepers (12) 


where the meaning of symbols is as before. __ 
In the above formulas if X,, Xs are substituted 
for R, and R,, respectively, then 68 will be obtained. 


4.3. Dielectric Coated Waveguides 


A metallic surface exhibits a surface impedance 
given by [8] 
Zm=(1+J) Rn (13) 


where R,, is given by (6). The surface of a copper 
tube (c=5.8X10’ mho/m), therefore, presents at 
\,=8.7 mm a surface impedance of (1+ 7) X1.28X 
10-*, relative to free space (about 0.048 ohms). 

If, now, the metallic tube is coated with a thin 
layer of dielectric then its surface impedance will be 
enhanced as given in table 2, where ¢ is the thickness 
of the dielectric layer, k,=2z7/d., €, relative permit- 
tivity, and 6 the loss angle of the dielectric. Clearly, 
by choosing suitable materials and the thickness of 
the dielectric layer, it is possible to obtain within 
wide limits any required surface impedance compo- 


The surface impedance of a dielectric coated metal 
surface 


TABLE 2. 


These value must be augmented by Zm=RmtjXm (Rm=Xm= €o/u0,) 
the surface impedance of the metal. 


Surface 
impedance Re part Im part 
Z.=RAjX, | R.=tke(tan 3/e) X.=tko(1—1/e,) 
Rem tan 8) 
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nents and, consequently, the attenuation and phase 
velocity of the H and H modes can be varied rela- 
tive to the H,, modes. By way of illustration, at 
the free space wavelength of 8.7 mm, for «,=2, 
tan 6=0.15, and t=0.046 mm, the surface impedance 
components become 7 1.710-? and 
Z=3.7 3.3 X107* over and above the value 
of Z,. If we substitute these values in eqs (11) and 
(12) we find that for a waveguide of 7-cm diam, 
there is a 20-fold increase in attenuation of FH and H 
modes (other than H,,) relative to the Hp, at the 
expense of only a 3 percent increase in attenuation 
of the Hj, mode. The surface reactance X, is, at 
the same time, increased by a factor of 140 leading 
to a substantial change in phase velocity of most 
modes. 

The increase in relative attenuation is necessary 
to cope with mode conversion-reconversion, due to 
irregularities in the waveguide, while the increase in 
the surface reactance, X,, is necessary so that the 
waveguide can cope with bends. 


4.4. Ring or Disk Structures 


Any periodic array of coaxial disks or rings of 
different surface impedances falls into this class. If 
the pitch of such a structure is small in comparison 
with the wavelength then a simple approximate 
treatment is to ng as follows: 

Let Z, and Z, be the surface impedance of the 
individual elements, thickness ¢; and ¢,, respectively, 
then the anisotropic components of the effective 
surface impedance are given by 


t Zot, 
titte 


Z,Z;(ti + te) 
Zt; +Zits 


Z, 


(14) 


Since in most applications |Z,|>|Z,| and t-<<t, 
these expressions can be further simplified to 


ts) 


(15) 

Thus given Z;, Z., ti, and t, the surface impedance 
components Z, and Z, can be calculated. In turn, 
the attenuation and phase propagation coefficients 
can be obtained from (11) and (12). 

The quantities Z; and Z, are determined from the 
astealins of the actual physical structure. This can 
be a lengthy calculation, but since in practical 
applications only an estimate of the value of the 
surface impedance is needed, it is permissible to make 
a number of rather drastic approximations. For 
example, for a corrugated waveguide (fig. 4), Z; is 
the copper intrinsic impedance and Z, is the input 
impedance to a parallel plane transmission line, 
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Figure 4. Corrugated waveguide. 


length 7, and separation between plates ¢,. Thus 


(y+ te) | +320 


where Z,, is the surface impedance of the metal 
(e.g., eq (13)). 

It will be observed that Z; has been increased a 
little, in proportion of (t,+¢,)/t, but that Z, is highly 
reactive. Such waveguides are particularly suitable 
for negotiating bends. 

Clearly, other structures can be calculated in an 
analogous manner. 


(16) 


4.5. Helical Waveguides 


It will be observed from eqs (15) that Z, and Z, 
can be adjusted within very wide limits to any desir- 
able values provided that the proportions of the 
microstructure of the waveguide surface are suitably 
chosen. This feature makes the disk waveguides, 
and their variants, extremely attractive. But, un- 
fortunately, such waveguides are difficult if not 
impossible to manufacture in quantity. It is for 
this reason that helical waveguides are preferred. 

It can be shown that if the pitch of the helical 
waveguide is very small, as would be with waveguides 
wound with a fine wire, and the surface impedance 
not too large (say less than 107', i.e., some 40 ohms), 
then the helical waveguide behaves as a correspond- 
ing ring structure. 

A particularly important structure is shown in 
figure 5; this shows a helix of insulated (copper) wire 
set in a suitable resin and backed by a resistive layer 
of surface impedance R, the whole structure is en- 
closed in a suitable protective jacket. If the wire 
diameter is ¢; and the gap between the wires ¢., then 
provided that ¢, is not too small, the structure 
behaves approximately as a corresponding disk 
structure (fig. 6) with /=t. 

Using the procedure as outlined before, eqs (15) 
can be used where Z,=Z,,=copper intrinsic im- 
pedance and Z, is taken equal to the input impedance 
of a transmission line length /, terminated in a 
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5. 


A helical waveguide. 
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Figure 6. Ring structure with a resistive layer. 


resistive impedance R. Clearly, therefore, 


te ty 
(17) 
Li=Lu 
1 


As a numerical example consider a helical waveguide 
with following constants: t=4h, 
Rn»=10~ (approximately copper impedance), R= 107? 
(approximately 4 ohms) then, neglecting the lay angle 
of the helix, at a free space wavelength of 8.7 mm 
the surface impedance components become 


Z,=0.003+ 0.021 
Zo=1.5(1+j) X 10-4. 
The propagation coefficients follow from eqs (7) 


and (8). Thus for a waveguide 7 cm in diam there is 
a 50 percent increase in the attentuation of the Ho, 
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mode, and a 30-fold increase in the attenuation of 
most H-modes. The attenuation of the lower order 
H-modes has also been increased by a factor of 30, 
while the surface reactance is increased 200 times. 


5. Basic Theory of Nonuniform Waveguides 


Theory of nonuniform waveguides is fundamental 
for thorough understanding of transmission aspects 
as well as design of many waveguide components 
peculiar to long distance transmission by waveguide. 
Thus, commonly, waves are conveyed from an os- 
cillator via a single mode rectangular waveguide 
through a transducer to a circular waveguide then 
through a taper to the main waveguide run. 

We shall see that a mode transducer, a bent wave- 
guide, as well as a taper are particular embodiments 
of nonuniform transmission lines. The main wave- 
guide run is also a nonuniform transmission line in 
the sense that it can carry a large number of modes 
and all minute irregularities act as coupling elements 
between the lines; each line representing one par- 
ticular mode. 

In general, we represent a nonuniform waveguide 
by its equivalent circuit: » coupled transmission 
lines, where n is the number of possible modes of 
propagation (see fig. 7). 


Z 
UNIFORM NON -UNIFORM UNIFORM 
WAVEGUIDE WAVEGUIDE WAVEGUIDE 
z=0 Z=L 
Figure 7, Schematic of a nonuniform waveguide. 


For convenience in analysis, nonuniform wave- 
guides can be subdivided into (1) straight wave- 
guides, and (2) bent waveguides. 


5.1. Straight Waveguides 


The first general solution to transmission in non- 
uniform waveguides was obtained by Stevenson [10] 
in terms of a series of transverse wave functions. 
Through such an approach it is possible to obtain 
equivalent generalized telegraphist’s equations [11] 
of the waveguide. Adopting Reiter’s approach 
[12, 13], it can be shown [17] that for the forward 
and backward traveling waves the following equa- 


tions hold 
d Aj 


1 d(in K;) 


43 + 38,47 + S845) 
~ Pp 


(18) 
where z is the coordinate along the waveguide axis. 


A* are amplitude coefficients of the forward (+) and 
backward (—) traveling waves on the ith line 
(corresponding to the ith mode). 

*® are the corresponding coupling coefficients be- 
tween the waves on lines 7 and p (a parameter 
related to the geometry of the waveguide). 

K; is the wave impedance of the ith line. 

8; the propagation coefficient. 


The above are exact equations for a set of lossless 
nonuniform transmission lines where S, A, and 8 are, 
in general, functions of z. These equations are exact 
representation of transmission in a nonuniform wave- 
guide. In principle, therefore, given the functional 
dependence of the waveguide geometry on the coordi- 
nate z, the propagation problem can be solved, but 
the computational difficulties are in most cases 
formidable. With practical waveguides, however, it 
is possible to proceed with a number of simplifying 
assumptions leading in many cases to rather simple 
solutions [14, 15, 16, 17]. 

In the first place, we assume that the waveguide 
geometry changes slowly and consequently S, K, 
and 8 are slowly varying functions of 2. The 
coupling between the lines per unit length will there- 
fore be small. If now the input at one end of the 
waveguide is a pure mode, amplitude A,, then the 
boundary conditions are simply 


A;,(L)=0 
A; (L)=0 


A;,(0)=A, 
A; (0)=0 


(19) 
(t#m) 


With these assumptions only the coupling between 
the line m and the remaining lines need be considered; 
the cross coupling between the lines can be neglected 
and so can the reflected waves. We then obtain 
a simplified solution in the form [17] 


At(z)=A, exp Bu | (20) 


The amplitude of the remaining waves are given by 


L L Z 
0 


Aj (0) 0 0 0 
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(21) 
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(21) 


In the derivation of these expressions it has been 
assumed that none of the modes are at the cutoff. 
When dealing with components, such as tapers for 
example, some modes are of necessity at the cutoff at 
some cross section along the taper, and it would 
appear that (20) and (21) are not a sound basis for 
design. Nevertheless, it appears that the use of (20) 
and (21) is usually justified [18]. 

It is to be observed that usually it is only the 
lower order modes that need to be at all considered 
(typically Hew, Ap, Hx, cf., section 4.1.). 

The 6 coefficients of these modes differ but very 
little from that of the Hy, and consequently (eq (21)) 
only the coefficients of the forward traveling waves 
(At) need be considered; the coefficients A~ are very 
much smaller and can, therefore, be neglected. 

Using eqs (21) an interesting proposition can be 
proven. In general, a pure mode at the input to a 
waveguide will produce at the output a combination 
of modes. For a good transducer the output consists 
of an almost pure mode of the desired kind. The 
purity of this mode can now be made as high as 
desired merely by increasing the length of the trans- 
ducer: the longer and more gradual the transducer 
the purer is the mode at the output. 


5.2. Propagation in Bent Waveguides 


The general coupled transmission line equations 
for a bent waveguide of arbitrary cross section were 
given by Shimizu [19] from which the relations for 
a circular waveguide follow. The particular case 
of a bent circular waveguide was analyzed earlier 
by Jouguet [20], who analyzed the field in a curved 
waveguide by means of perturbation calculus and 
obtained the field in terms of an expansion in powers 
of 1/R, where RF is the radius of curvature of the wave- 
guide. 

It can be shown that if a curved section of wave- 
guide is interposed between two straight waveguides 
(fig. 8) and a pure H, mode is incident at one end of 
the curved waveguide, then the wave appearing at 
the other end is composed of the Hq, E\,, and a series 
of H,, modes [19, 20, 21], but the exact proportion 


PURE Ho, — MODE 


Figure 8, Curved circular waveguide. 
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of each of the modes depends on the geometry of the 
bend. The equivalent circuit of the bend is there- 
fore n+2 coupled transmission lines: one for the Hy, 
one for the £j;, and n lines one for each of the Hi, 
modes. The coefficients of the coupling matrix 
are constants related to the geometry of the bend. 
But it can be shown, further, that if the radius of 
curvature is large, the series of the Hi, modes may 
be neglected, since these are of 1/R? or higher order, 
and therefore only the Hy, and E,,; modes need be 
considered. The equivalent circuit is then, simply 
two coupled transmission lines. 

Jouguet’s analysis can be further extended to wave- 
guides of arbitrary surface impedance [2]. It 
transpires that provided the radius of bend is not 
too small, the field inside a curved waveguide can be 
represented by 


(22) 


that is the field is a linear combination of the Hy; 
and £,, modes in a definite proportion. The pro- 
portionality factor can take two possible values 
given by 


7?—2pr—2=0, (23) 
where 
jZs. (24) 


If the surface impedance and the radius of curvature 
are small then p is small and 7 is approximately 


given by v2; therefore, the expression for the field 
[2] can be put into the form 


V= M, COs sin (25) 


where 


ko 
6Br= 
and hokey 2 (26) 


‘ 
(Bot +5 (azg+ay) 


with @ the phase propagation coefficient of the straight 
waveguide, ag and ay, the attenuation coefficients 
of the Z,, and Hy, modes in the straight waveguide. 

Two important conclusions must be made: (1) 
The energy is carried alternately in the Ho, and F,, 
modes, and (2) the mean attenuation of the complete 
normal mode (25) is equal to the mean of those for 
the £,,; and Ho, modes. 

Evidently, the waveguide is unsuitable for the 
transmission of the Hy, mode around the bends. 
However, the examination of the above equations 
shows that if p is a large real number a different 
situation arises. Large real values of p are obtained 
for highly reactive surfaces (X, a large number), 
and in such waveguides the two modes Hy, and FE), 


4 ef stands for the field of the Ho mode, and &j stands for the field of the 
Ey, mode. 
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propagate almost independently of each other, but 
the attenuation of Hy, mode is slightly increased. 
In such a waveguide, the attenuation of the Hy 
mode is given by [2] 


(27) 


where a is the attenuation in the straight waveguide 
and is given by (eq (12)) 


=aH+aR 


“=D Rg. (28) 
The additional attenuation due to the bend is given 
by 
D R,/\Z,|? 


Clearly, by choosing a sufficiently large surface 
reactance it is possible to bring this increase in atten- 
uation below any predetermined level. Such values 
of surface reactance can be conveniently obtained 
either by coating the waveguide surface with a layer 
of dielectric or by corrugating the waveguide surface 
in a suitable manner. 

The ratio az/a>=ve(R,/Rs) defines a quality factor 
of the bent waveguide. With the help of figure 9 it 
is possible to determine the required value of surface 
reactance for a given value of radius of curvature and 
quality factor vz. Clearly, helical waveguides of 
appropriate proportions (sec. 4) can be designed for 
the purpose of negotiating bends. 


(29) 
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Ficure 9. Quality factor, v, of a bent waveguide of 7-cm diam- 
eter as a function of wavelength, with |Z,| R as parameter. 
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With suitably designed waveguides, bends even of 
a few hundred feet in radius can be negotiated quite 
successfully. A further improvement can _ be 
achieved by increasing gradually the radius of bend 
by means of a suitable transition curve [22] (a 
“tapered transition’”’). This is yet another embodi- 
— of a set of nonuniform coupled transmission 
ines. 


6. Transducers 


A transducer is a waveguide configuration which 
achieves a required change of mode pattern and/or 
pattern size. 


6.1. Tapers 


A taper is the simplest illustration: it achieves 
a change in waveguide diameter and, if properly 
designed, leaves the mode pattern unaltered yet 
changed in size. 

The rigorous theory [23, 24] of tapers can be very 
involved indeed but for our applications a number 
of approximations are permissible resulting in an 
adequate and simple design procedure. 

Consider two coaxial waveguides joined by a 
conical taper and suppose that the larger waveguide 
(radius a) is operated substantially above cutoff 
(overmoded waveguide) and the smaller waveguide 
(radius 5) carries a pure Hy, mode. The taper will, 
evidently, give rise to a series of unwanted modes 
besides the Ho, mode, but because of the symmetry 
of the structure only the H,, modes will be present. 
If, now, the diameter of the smaller waveguide is 
sufficiently small to permit only the Hy, mode to 
propagate, the amplitudes of the H,, modes in the 
arger waveguide will be a simple function of the 
parameter o given by 


(30) 


where Z is the length of the taper and X the free 
space wavelength. 

These relations are obtained from simple analytical 
considerations [25] by expanding the fields in the 
circular waveguides in cylindrical harmonics and the 
field in the conical region in spherical harmonics and 
equating the fields in the cross sections radius a and 
b, respectively. 

It can be shown that out of the modes produced by 
the taper the Ho, mode is largest in amplitude. 
Considering, therefore, the conversion into the Hy 
mode only we get for the taper length 


0.57 a 

L ie i (a—b) (31) 
where p is power conversion coefficient, which is the 
ratio of the permitted power in the Hy, mode to that 
in the Hy; mode. 

For some applications the length of such a taper 
may be excessive, but a further improvement is 
possible by using a two-section conical taper (fig. 10). 
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Figure 10. Schematic of a two-section conical taper. 


For such a taper [26] by choosing correctly the pro- 
portions of the two cones the performance of the 
taper can be considerably improved. Optimum 
performance is possible only at one frequency, but by 
correct design it is possible to keep the sum of the 
energy carried in the Hy, and Hh; modes below any 
required level, over a substantial bandwidth [26]. 

More elaborate design procedures have been de- 
vised by several investigators [18, 24, 27]. But it 
appears [28] from results on antenna theory that 
provided the taper contour is “smooth and gradual” 
and that the derivatives at the taper ends are zero, 
there will be very little mode conversion even over 
a wide bandwidth. 


6.2. Mode Transducers 


With the exception of a few isolated instances 
there is no general design procedure for mode 
transducers. In a few simple cases it is possible to 
predict from the changes in the waveguide geometry 
the type of mode produced and with a little engi- 
neering intuition to design sometimes a reasonably 
satisfactory transducer. In most cases, however, 
even the intuitive approach is of no avail. Recently 
a satisfactory design procedure has been proposed 
[29] and confirmed experimentally. This is ex- 
plained below. 

Consider a mode transducer connecting two uni- 
form waveguides A and B. According to an estab- 
lished theorem, if the change from waveguide A to 
waveguide B is effected very gradually then a pure 
mode in waveguide A will produce a substantially 
pure mode in waveguide B: the purity of the mode in 
waveguide B will be the higher the longer the 
transducer. 

To determine the modes present at the ends of the 
transducer, it is necessary to choose several cross 
sections along its length and to determine the eigen 
functions in each. If the number of cross sections 
chosen is sufficient, then from the similarity between 
the successive field patterns of field distribution it is 
possible to examine the passage of any given mode 
along the whole length of the transducer, in detail. 

To design a mode transducer, the desired modes in 
the uniform waveguides are given, and it is then 
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required to determine the surface of the transducer 
for reasonable mode purity. We then choose a cross- 
sectional wave function VY (as a continuous function 
of z) such that it becomes the desired mode at each 
end of the transducer (at z=0 and z=L, where L is 
the length of the transducer). This wave function 
determines the surface of the transducer; and, 
provided that the function WV is a slowly varying 
function, the modes involved will be substantially 
pure. The design procedure is essentially numerical 
or graphical. 

As an illustration consider the Hp, circular to Hy, 
rectangular mode transducer. For convenience and 
better illustration we propose to design the trans- 
ducer in two parts: (1) Hy circular to Hy rec- 
tangular, and (2) Hy rectangular to H, rectangular. 

For the first transducer a convenient choice for the 
wave function is 


Yn= (DIA EVE cost (82) 


where to satisfy the requirements at each end of the 
transducer, the conditions are 


gi(0)=1 gi(I1)=0, 
g2(0)=0 = 


with L,;=the length of the transducer and 2z,y= 
Cartesian coordinates in the transverse plane. 

The diameter, d, of the circular waveguide and the 
width 2w of the rectangular waveguide are adjusted 
so that their cutoff numbers are equal (d/w=2.44). 

Evidently the function V,, reduces at the ends of 
the waveguide to the eigenfunctions of the circular 
and rectangular waveguides respectively. The sur- 
face of the transducer for any intermediate values of 
z is most easily determined graphically as follows. 
For any required value of z, curves of V,,=a con- 
stant are plotted: these represent the lines of electric 
field. But a permissible boundary must be normal 
to such lines and, therefore, the cross-sectional 
contour can be mapped. Figure 11 shows four con- 
tours of such a transducer, and since we know that 
the transducer must be gradual, in most cases four 
contours are sufficient to determine the surface of 
the transducer. 

Since the wave functions of the Hy. and Ho, modes 
have the same mathematical formula the obvious 
choice for the wave function of the second trans- 
ducer is [30] 


(33) 


V,,=COos — 


The curves V,,=a constant are simply straight 
lines and, therefore, the transducer has the shape 
shown in figures 12 and 13; the dimension 6 changes 
slowly from 6=0 to 6=h. 

Since h is not equal to the diameter of the circular 
waveguide a further taper is necessary to connect the 
two transducers. Thus the complete Hp; circular to 
Hw rectangular transducer is made up from the 
above three parts. 
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Schematic showing a gradual mode transformation 


Ficure 11. 
from the Ho; circular to the Ho2 rectangular (only one quad- 


rant showing). 


2w 


Fiacure 12, Cross section of anJHo: to Ho2 rectangular 
transducer, 


Perspective view of the complete Ho; to Hoe 
rectangular transducer. 


Figure 13. 


7. Slightly Irregular Waveguides 


The general problem of propagation in slightly 
nonuniform waveguides was investigated by Katzen- 
elenbaum [31] and in a number of specific cases by 
Morgan [32] and Tiguchi [33]. 

In the case of a waveguide carrying an Ho, mode 
at a frequency substantially above cutoff, a quasi- 
optical approach [34] can be adopted leading to 
simple derivation of relevant formulas. The problem 
can be stated with reference to figure 13. The cross 
section of waveguide A is circular and that of B 
differs a little but in an arbitrary manner. By 
assuming that the field in the aperture is substan- 
tially the field as produced by the Ho, mode in the 
waveguide A according to geometrical optics, simple 
integral relations can be derived for the intensity of 
the modes produced in waveguide B. There are 
only H-modes produced and the amplitude of the 
H; mode is given by [34] 


So 


where S, is the aperture, V7, is the gradient operator 
in the transverse plane, V; is the wave function of 
the H; mode in waveguide B, and Wp; is the wave 
function of the 4, mode in waveguide A. 

In the case of a junction formed by two offset 
circular waveguides of identical diameters only H,, 
modes will be generated if the exciting mode is a 
pure HZ), mode. If there is a small change in wave- 
guide diameter and the waveguides are coaixal then 
only H,, modes are produced. Similar deductions 
can be made for other discontinuities without 
actually solving the expression (35). 

For the two coaxial waveguides of slightly different 
diameter (see fig. 14) eq (35) leads to particularly 
simple relations. In this case the amplitude (A,,) 
of the H,, mode produced by the incident Zo, mode 
is given by 


Aon _ 
An (36) 
where x; is the zth root of J,(z)=0. 
(a) 
WAVEGUIDE A WAVEGUIDE 8 
> N 


(b) 


The junction of two waveguides of slightly different 
cross section. 


Figure 14, 
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Consider now two such steps in tandem. Clearly 
the contributions from the conversions must be 
combined taking the correct phases into account. 
It can then be shown that if the inserted length of 
the waveguide is one-half beat wavelength (of Hp, 
and Hy.) long and its diameter the mean of the two 
extreme waveguides, then the mode conversion is 
minimized [35]. The same procedure can be applied 
to other H,, modes and a suitable transition between 
two waveguides of different diameters designed. 


8. Mode Conversion-Reconversion 


As mentioned in section 2 mode conversion- 
reconversion is a phenomenon existing in any over 
moded waveguide, such as would be used for long- 
distance transmission. The level of mode conversion 
or reconversion on their own is a function of wave- 
guide tolerances, but the combined effect of mode 
conversion-reconversion, which determines the signal 
distortion is also a function of waveguide design or 
more precisely it is related to the surface impedance 
of the waveguide. 

For technological reasons one cannot reduce the 
mode conversion process beyond a certain limit, 
related to the waveguide tolerances. For given 
tolerances, however, the mode conversion-reconver- 
sion process can be reduced, in principle, below any 
predetermined level by choosing suitable surface 
impedance. The remaining problem then is a 
technological one, of realizing such surface impedance. 
The actual design is a compromise between what is 
desirable from the mode conversion-reconversion 
point of view and what is practicable or at all pos- 
sible to achieve using existing methods. 

Briefly, for good performance on bends, wave- 
guides must have high surface reactance in the axial 
direction (X, components) while the circumferential 
component (Xz) remains small. For good discrimi- 
nation against unwanted modes produced by typical 
small irregularities, large attenuation for unwanted 
modes and small attenuation for the Hy; mode is 
called for: that is large R, and small Ry. 

Due to waveguide irregularities a large number of 
modes become coupled together. If the coupling 
coefficients representing the aggregate effect of wave- 
guide irregularities contain a periodic component, 
considerable signal distortion takes place [4]: in effect 
the B(w) characteristic becomes broken up by regu- 
larly spaced absorption bands. But if the irregu- 
larities are spaced in a random manner no such 
extreme signal distortion will take place; but on the 
other hand, the @ and £6 characteristics become 
irregular in a random manner with a root mean 
square value proportional to (¢/Ag)? (ec is the rms 
value of the coupling coefficient c, and Ag the differ- 
ence between the phase propagation coefficients of 
the Hy mode and the parasitic mode) leading to 
effective crosstalk between the components of the 
coded signal. In general, the signal distortion due 


5 The conversion to Ho, modes has been neglected, since this is not usually 
troublesome. 
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to mode conversion-reconversion is a function of the 
quantity (c/AB)?/Aa where Aa is the difference be- 
tween the attenuations of the parasitic mode and 
the Hy, mode. For this reason the more the group 
velocities of the parasitic modes differ from that of 
the Hy, mode and the greater the attenuation of the 
parasitic modes, the smaller the signal distortion due 
to mode conversion-reconversion. As discussed in 
section 4, special waveguides can bring about such 
changes and hence their significance. 


9. Conclusions 


Waveguides for long distance communications 
already form a wide specialized subject, too wide to 
be discussed adequately in an article of this size. 
But we have discussed the main constituents of a 
waveguide communication system pointing out the 
necessity of regenerative modulation method, such 
as PCM. 

Special waveguides have been discussed in some 
detail because, in the light of present knowledge, 
such waveguides—either of the helical type or di- 
electric coated—are the only types likely to succeed. 

Mode conversion phenomena are intimately con- 
nected with overmoded waveguides, of the type used 
for long distance transmission. Deep understanding 
of mode conversion processes is essential for success- 
ful designs and accordingly a theory has been out- 
lined. The design of such components as tapers, 
mode transducers, and bends rest on mode conversion 
theory. 

Various sources of signal distortion have been 
analyzed in some detail, but it has been stressed 
that if PCM modulation is used and regenerative 
repeaters, as well as special waveguides, a most 
advanced system of unequalled capacity will result. 


Most of the work described in this paper was 
carried out as part of a research program of Standard 
Telecommunication Laboratories Ltd. The authors 
have benefited from discussions with a number of 
their colleagues: such discussions helped to crystalize 
many of the more involved ideas. 

Acknowledgment is also made to Standard Tele- 
communication Laboratories Ltd., for permission to 
publish the paper. 
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Useful Radiation From an Underground Antenna 


Harold A. Wheeler 
(February 26, 1960) 


An underground antenna delivers power to the surrounding conductive medium, and a 


fraction of the power goes out as radiation above the surface. 


his fraction is denoted the 


“radiation efficiency.”’ It is expressed in simple terms for two types of underground anten- 
nas. The first and simplest is a vertical loop in a submerged spherical radome. The second 
is a submerged horizontal insulated wire with each end connected to a ground electrode. 
In each case, the efficiency is the product of three simple factors: The first depending on the 
index of refraction between air and ground; the second proportional to the size (radius of 
the radome or length of the wire) ; the third giving the attenuation with depth. An example 
for 1 megacycle per second gives an efficiency of 0.0014 for an underground wire of specified 
dimensions, The radiation efficiency is applicable to sender or receiver. 


For protection from hazards above the ground, an 
underground antenna may be considered for some 
purposes, such as communication over a moderate 
distance. In a sender, such an antenna radiates a 
small part of its power out of the ground into free 
space. In a receiver, it intercepts a small part of 
the wave power that penetrates into the ground. 
Usually the most favorable path between two anten- 
nas is in the space above ground, even if both of 
the antennas are underground. 

It is the purpose of this note to give some simple 
formulas for the radiation efficiency and some other 
properties of underground antennas of two simple 
forms. These are based on previous publications 
of R. K. Moore [1] *, J. R. Wait [1a], and the writer [2 
to 6], which were directed to underwater antennas 
presenting essentially the same problems. 

By formulating the efficiency of radiation out of 
the ground, separate from other losses in the antenna 
and in propagation, the formulas are simplified and 
the contributing factors are made apparent. These 
formulas are used in computing the performance of 
an underground antenna to be described by L. E. 
Rawls [7]. They are also applicable to the trans- 
mission formulas for computing system performance, 
as proposed by C. R. Burrows [8], H. T. Friis [9,10], 
and K. A. Norton [11]. This “separation of the 
variables” is usually helpful in understanding the 
respective factors in performance. 

The radiation pattern of the cases here to be 
considered is that of a small vertical loop just above 
ground. The reference for evaluating radiation 
efficiency is a perfect small loop just above a perfect 
ground plane, as shown in figure 1(a). The power 
supplied to the loop is all radiated; this useful 
radiated power is symbolized by P,. In a sender, 
the loop image doubles the far field. In a receiver, 


! Contribution from Wheeler Laboratories, Great Neck, N.Y., and Develop- 
ment Engineering Corp., Leesburg, Va. 

2 Paper presented at Conference on the Propagation of ELF Radio Waves, 
Boulder, Colo., Jan. 26, 1960. 

8 Figures in brackets indicate the literature references at the end of this paper. 
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Figure 1. A perfect small loop as a basis for expressing the 
radiation efficiency of an underground antenna. 


a. Loop over perfect ground plane. 
b. Loop in submerged radome. 


the total interception area (3/2 radiancircle) is only 
half associated with the actual loop above ground, 
so it has an interception area of only % radiancircle, 
as referred to the power density in a plane wave 
along the ground. : 

The simplest case of an underground antenna is 
a small loop in a small spherical radome located 
underground at a certain depth, as shown in figure 
1(b). The radome, of radius a, is filled with perfect 
insulating material (such as air), and its center is 
submerged to the depth d below the surface. 

The relative properties of the ground and the 
space above (air) are indicated by the radianlengths 
in the respective mediums, 6 and \/27. Assuming 
that the ground behaves mainly as a conductor 
rather than a dielectric, the former (5) is equal to 
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the well-known skin depth. All properties are to be 
expressed simply in terms of ratios of the several 
“length” dimensions. 

The “index of refraction’ may be defined as 


n=)/2n8>>1 (1) 


which is the ratio of phase velocities in the respective 
mediums. It is assumed to be very great, which is 
conducive to propagation of a wave along the ground 
with little “wave tilt.” 

Assuming a perfect loop as the radiator in figure 
1(b), the power out of the loop (P.) goes into the 
ground, and a small part of it (P,) is radiated above 
ground. The ratio of these two quantities is the 


radiation efficiency: 
P, \x 

; (2) 

size depth 


immersion 

(1/n) (loss) (atten.) 
The three factors have the indicated significance. 
The middle factor is peculiar to this type, reflecting 
the exeess dissipation in the near magnetic field if the 
size of the radome is much less than the radiansphere 
in the ground. 

The same value of efficiency is applicable in a 
receiver by reciprocity. In fact, the efficiency is 
more easily derived for a receiver, since plane waves 
may be assumed above and below the surface. It is 
then necessary also to compute the resistance appear- 
ing in the loop as a result of the surrounding ground, 
and this is easily done for the small alesiant radome. 
In a receiver, the “noise temperature” of atmospheric 
radio noise above the surface may be multiplied by 
this efficiency to give the apparent temperature of 
such noise reaching the radome interior. This 
apparent temperature may be added to the actual 
temperature in the ground around the radome. 

Figure 2 shows another type of underground 
antenna, for which the radiation efficiency and some 
other properties will be evaluated. The radiator is a 
long insulated wire buried in the ground in a hori- 
zontal position. Each end is connected to a ground 
electrode. In a sender, the current in the wire is 
provided by a generator connected in series. 

For simplification, it is assumed that the wire has 
a small radius (r<<6) and a length (J) much greater 
than the skin depth but much ed than the wave- 
length in air (375<1</z); also that the depth (d) is 


P, 
d>8 
d > 
28 
INSULATED WIRE \\ 
1 >> 8 


Figure 2. A horizontal wire as an underground antenna. 
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greater than the skin depth (as indicated). Each 
ground electrode is assumed to be “perfect,” ap- 
proximating a large plane perpendicular to the wire. 

On the same basis as the preceding loop antenna, 
the long wire has the radiation efficiency: 


306 


size 


(gain) 


(3) 


As before, the middle factor is peculiar to this type, 
but here there is a gain instead of a loss, in proportion 
to the length of the wire. 

The resistance and reactance (R and X) of the 
underground wire, as derived by R. K. Moore [1], 
are expressed as follows: 


1 1 ,, 2nl 2rl l 

(5) 


The wave resistance of free space is represented by 
R,=377 ohms. 

The resistance is that of the ground return circuit. 
It is equal to the resistance of a cross section 26X26 
of ground conductivity, on the basis of uniform cur- 
rent density, as indicated in the diagram. This 
resistance is involved in evaluating the efficiency. 
If the wire depth is less than the skin depth, the 
resistance changes very little. In fact, an older 
derivation by i R. Carson [12] shows that the 
resistance has the same value for the special case of 
a wire located on the surface. 

The reactance is that of a coaxial line comprising 
the wire in a circular shield of radius about 0.8 6. 
The reactance is given for a wire of perfect con- 
ductivity and complete skin effect. For an actual 
wire small enough to neglect the skin effect, add % 
in the brackets in (5). If the depth is less than the 
skin depth, the reactance increases slightly; if on 
the surface, add '4 in the brackets in (5) for this 
change. 

Implicit in the radiation efficiency is the effective 
height of the loop formed by the wire and the 
ground return. The radiation above ground is the 
same as if there were a loop of the following effective 
height, with its center at the depth of the wire: 


h=6/y8. (6) 


It is remarkable that this is valid for any depth of the 
wire below the surface. 


An example of the long wire underground is 
described by the following dimensions and computed 
properties: 


f=1 Mess, 
wavelength \=300 m, 
ground conductivity_. o=0.01 mho/m, 
skin 6=5 m, 
d=1 m, 
/=100 m, 
r=3 mm, 
index of refraction____ n=10, 
effective height____-__ h=1.8 m, 
R=99 ohms, 
¥=900 ohms, 
power p=0.11, 


radiation efficiency_P,/Pa=0. 0014=—28.5 db. 

In this example, for each 1 kw of power into the 
ground, only 1.4 w will end up in useful radiation 
above the surface. 

The efficiency can be increased by locating a 
number of wires (N) in parallel in a grid, and sup- 
plying all with equal currents to form a sort of 

“current sheet.’”’ Ifthe wires are spaced more than 
the skin depth (or more than twice the skin depth if 
near the surface) the efficiency is multiplied by the 
number (N) approximately. 

Two practical problems have been ignored in the 
above. One is the power loss in the end connections. 
The other is the current variation along the wire. 

The ground electrode for each end connection 
should have enough surface to reduce its resistance 
much below the irreducible value of the ground 
return (given above). 

Some current variation along the wire is caused 
by the voltage applied to the capacitance to ground 
through the insulation. The voltage developed by 
the resistance is inevitable. However, the voltage 
developed by the reactance can be reduced by 
inserting series capacitors along the length. In the 
present case, for example, a capacitor of 50-ohms 
reactance (about 4 R) may be inserted at intervals 
of 50-ohms wire reactance to tune out this reactance 
at one frequency. This would require 18 series 
capacitors in the wire. 

The straight long wire does not have the highest 
efficiency possible for its length. A narrow loop 
of the same length, enclosed in a cylindrical radome, 
has greater efficiency by virtue of its wire return 
instead of a ground return. However, this loop 
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has a much smaller power factor, so it must be 
—— tuned and has a relatively narrow band- 
width 

These formulas for radiation efficiency and other 
properties are intended to indicate the behavior 
of certain types of underground (or underwater) 
antennas. They are expressed in simple terms, and 
they fit into the usual computations of terminal 
circuits and point-to-point transmission. 

The earlier work of Carson [12] and Moore [1] 
has contributed much to this subject, as mentioned 
in the text. The recent studies have been made as 
an aid to practical designs being proposed by Devel- 
opmental Engineering Corp. 
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Observation of F-Layer and Sporadic-E Scatter 
at VHF in the Far East’ 
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This paper describes properties of sporadic-F scatter and F-layer scatter observed over 
the Okinawa-to-Tokyo path (1480 kilometers) and the Philippines-to-Tokyo path (2850 
kilometers) operating at frequencies of about 50 megacycles per second. 

Sporadic-E scatter is often observed on the Okinawa signal in the evening hours and 
has the closest correlation (0.94 in correlation ratio) with the occurrence of sporadic FE 
characterized by the descriptive symbol ‘‘M”’ of all ionospheric factors. Bearing of the Hs 
scatter shows a regular diurnal variation similar to that of the normal E-layer scatter. 

F-layer scatter generally appears on the Philippine signa! in autumn when the F-layer 
at the path midpoint displays an anomaly denoted by the symbol “‘R” or ‘‘S” having a top 
frequency of higher than 14 megacycles per second. A pulse test exhibited a pattern of 
multipath signals extending over more than 1 millisecond. Bearing of the F-layer scatter, 
an evening phenomenon, gradually deviates westwards from the great-circle path with the 
lapse of time. 


1. Introduction The orientations of both paths were almost from 
southwest to northeast. The transmitting antenna 
used at Okuma was beamed to 35°13’ E of N devi- 
ating northwards by 8.5° off the great-circle path to 
Tokyo. Surface path length between Okuma and 
Tokyo is 1480 km, and Poro and Tokyo is 2850 km. 
So one-hop Fs propagation may exist over the former 
circuit and two-hop Fs or one-hop /2 propagation 
over the latter. 

The Yamagawa and Okinawa ionosphere stations 
are located 330 km and 160 km from the path mid- 
points, respectively, as seen in figure 1. 


Great progress has been made in the study of 
transmission by way of inhomogenities in the iono- 
sphere, commonly termed scattering in the E-layer. 
In addition, F-layer scatter in the Far East has been 
pointed out by R. Bateman and others [1]? and a 
peculiar signal enhancement has been reported re- 
cently by the authors [2] as a kind of scatter associated 
closely with sporadic E. Continuous measurements 
of signal strength, bearing, and angle of elevation at 
about 50 Mc/s have been carried out in Tokyo on 
Okinawa and Philippine signals transmitted until the 
end of October 1959 as the continuation of a part of 
the United States program for the IGY. 

Aside from characteristics of Es propagation which 
will be reported as an interest in the field of com- 
munications engineering, the following will mainly 
describe the types of scatter mentioned above. 


2. Scope of Observation 


2.1. Signals Measured and Experimental Circuits 


An outline of the experimental circuits is shown in 
figure 1. Continuous waves of 2-kw output were 
emitted towards the middle part of Japan at 49.68 
Me/s from Okuma, Okinawa, and at 49.84 Me/s from 
Poro, the Philippines. Each of the transmitting 
antennas was a horizontal 5-element Yagi, having a 
gain of 7 db relative to a dipole, a height of 12 m 
above the ground, and a half-power beamwidth of 
54°. The measurements were conducted at Fukuoka 


receiving station, near Tokyo, from October 1958 mit-point | 
through October 1959. Pulse signals, from Poro, | |ocQiwa ro 
were observed in autumn 1959. 


: Contribution from Kokusai Denshin Denwa Co., Tokyo, Japan. Ficure 1. Map of the Far East, showing the forward scatter 
2 Figures in brackets indicate the literature references at the end of this paper. paths and the ionosphere stations. 
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2.2. Method of Observation 


Three 5-element Yagis A, B, and C of same con- 
struction as the transmitting antenna were used for 
the measurements. Yagis A and B were mounted 
30.6 m high on two iron poles separated by 105 m, 
directivities being shifted by 20° westwards and 
southwards from the true bearing of Okuma respec- 
tively. Yagi C was mounted at a height of 12.6 m, 
being beamed to the same direction as Yagi A. Bear- 
ing and angle of elevation were obtained from the 
ratios of voltages induced on Yagis A and B and on 
Yagis A and C, respectively. 

The ratio of voltages was continuously recorded 
using two sets of voltage ratio meters [2] which are 
capable of measuring angles of arrival accurately 
even on a low level of the signal strength. Tbe out- 
put voltage of the meter was logarithmically com- 
pressed and indicated by a pen recorded as open- 
circuit voltage measured at 75 Q impedance in 
decibels above 1 uv. 


3. Observational Results 
3.1. Okuma-to-Tokyo Transmission 


Figures 2 and 3 show typical results observed on 
the Okuma signal in winter and summer months, 
respectively. Figure 4 exhibits an example of ap- 
pearance of the pen records during a period including 
all types of scatter. 

In winter, the signal received is mainly the normal 
E-layer scatter Se et meteoric bursts throughout 
the day. Attenuation of the scatter in the ionosphere 
is estimated to be between 72 and 82 db in the median 
value, showing the maximum at about 1800 J.s.t. and 
the minimum around midday. 

As seen in figure 4, the signal strength sometimes 
rises to a maximum of 30 db above the level of the 
normal value, late in the evening. This phenomenon 
will be referred to as F's scatter and its characteristics 
will be discussed later. 

Hourly median angles of arrival are indicated by 
dots, and ranges between upper and lower quartiles 
are shown by lengthwise lines in figures 2 and 3. 
Bearing data explicitly show a diurnal variation 
which deviates westwards off the true bearing during 
the period of 1800 to 1000 J.s.t. and slightly south- 
wards early in the afternoon. Seasonal and diurnal 
distributions of hourly median bearings are exhibited 
in figure 5 representing contours of equal bearings. 
The figure shows a feature of regular diurnal vari- 
ation with a tendency of deviating westwards a few 
degrees off the true bearing on the average through- 
out a year. The fact that the main lobe of the trans- 
mitting antenna 1s beamed northwards off the great- 
circle path may be one of the causes of the average 
deviation of bearing. However, attention must be 
drawn to the fact that the diurnal variation including 
the average deviation is quite similar to W. R. 
Vincent and others’ results of best antenna angle for 
meteoric bursts [3]. The seasonal variation of bear- 
ing showed a complicated but periodic behavior. 
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Figure 2. Diurnal curves of signal strength, angles of arrival. 
Okuma-Tokyo, 49.68 Me/s, January 1959. 
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Diurnal curves of signal strength, angles of arrival. 
Okuma-Tokyo, 49.68 Me/s, June 1959. 


FIGURE 3. 


No noticeable variation of angle of incidence was 
observed in the median value. The minimum, how- 
ever, was clearly observed during a period of 1100 to 
1200 J.s.t. during all the months except summer, 
when the normal £-layer scatter was mainly received. 

In summer when Js reflections are most common, 
and produce an increase of many times the normal 
signal level, hourly values of angles of arrival are 
distributed in a wider range than those observed in 
winter as seen in figures 2 and 3. 

The range of quartiles of signal strength, bearing, 
and angle of elevation shows the minimum excursion 
at about 0600 J.s.t. early in the morning. 
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Bearings in degrees 


Angle of elevation in degrees 


Received voltage in db above 1uV 


FIGureE 4, 


Ficure 5. Contours of equal bearing. 
Okuma-Tokyo, 49.68 Me/s, October 1958-October 1959. 


3.2. Poro-to-Tokyo Transmission 


The signal from Poro was received only a small 
part of the time, except in summer and autumn. 
Judging from the characteristic appearance of spo- 
radic /, most of the signals received in summer are 
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A typical record of angles of arrival and signal strength observed on January 12, 1959. 


inferred to be two-hop propagation by way of the Es 
specular reflection. Most hourly median bearings 
observed in summer were nearly the true bearing. 
Hourly median angles of elevation were between 3.5° 
and 3.9° and they showed lower values in the day- 
time. These angles observed are slightly lower 
than the calculated by estimating height of sporadic 
E to be 90 km. 

On the other hand, it is obvious that the signal 
received in autumn was propagated via the F-region. 
Behavior of this signal will be described in the later 
chapter. 


4. Sporadic-E Scatter 
4.1. General 


As mentioned above, a remarkable signal enhance- 
ment is nightly observed on Okuma signal, sometimes 
attaining a level 30 db higher than that of the normal 
E-layer scatter. The authors called the event ‘“#s 
evening enhancement”’ [2], because it is closely asso- 
ciated with sporadic EF and it is likely to occur late in 
the evening. The enhancement appears to represent 
sporadic F scatter because of the following features. 
The Es scatter does not undergo as rapid a change in 
signal strength as sporadic / propagation and shows 
a gradual variation over an hour to several hours. 
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When this takes place, variations of bearing and 
angle of elevation are always observed. These varia- 
tions are quite different from that of Zs propagation 
as shown in figure 4. The Es scatter closely cor- 
relates with some type of sporadic E at the path 
midpoint rather than the magnitude of the foHs 
itself. The relationship is now explored. 


4.2. Rate of Occurrence 


Judging from a statistical study of foHs and 
attenuation in the ionosphere, I’, the signal of T <45 
db is assumed to be reflected by sporadic E and the 
signal of T>70 db as the normal E-layer scatter. 
An analysis was made of the signal of 45 db <I' $70 
db, the range in which most of the Es scatter is 
included. 

Seasonal and diurnal distributions of signals of 45 
db <Is70 db are shown in figure 6(a) indicating 
contours, of equal occurrence rates of reception in 
percentage, Nr. One of the most distinctive features 
in the figure is that the maximum rate appears late 


(a) Nr(45db <P <70db), 


(b) Nm. 
Ficure 6. Contours of equal rates in percentage of Nr 
and Nm. 


Nr: rate of occurrence of Es scatter, and Nm: rate of occurrence of Es affixed 
by symbol “‘M”’ at the path midpoint. 


in the evening. The authors tried to find some 
ionospheric factor which had a close correlation with 
the distribution in figure 7(a). This was found to 
be the occurrence of the descriptive symbol ‘“‘M”’ 
associated with the “critical frequency” of sporadic 
E (foEs). Figure 6(b) indicates contours of equal 
occurrence of such sporadic E, Ny, observed at the 
Yamagawa ionosphere station. Comparing the tem- 
poral distributions of figure 6(a) and figure 6(b), it 
is seen that there is a reasonably good correlation 
between the two sets of data. 


20 24 


Ficure 7. Contours of equal rate in percentage of foEs 


above 7.5 Mc/s at Yamagawa. 
Vr (foEs > 7.5 Me/s). 


For reference, the occurrence of sporadic E with 
foEs above 7 Me/s, Nr, is shown in figure 7. This 
figure is pertinent enough for the purpose of investi- 
gating correlation because Ny for any value of fos 
gives a similar distribution to that of figure 7. 
Figures 6(a) and 7 have a resemblance to each other 
in seasonal variation, but not always in diurnal 
variation. 

In order to obtain quantative relationship between 
Ny and Ny, hourly mean values of Nr and Ny, for 3 
months every season are plotted as a correlation 
diagram shown in figure 8. From the figure, the 
following relation is calculated: 


Nr=0.42 Ni, 


where correlation coefficient between them is 0.9 and 
correlation ratio is 0.94. 

An analysis using data observed at the Kokubunji 
ionosphere station, which is located 750 km from the 
path-midpoint, did not result in as good an associa- 
tion as the above. This may be because the station 
is too far away to give a good correlation. 


4.3. Angles of Arrival 


Needless to say, the measured angles of arrival 
are influenced by the signals from all the modes of 
propagation present. The data were, therefore, 


3M, Measurement questionable because the ordinary and extraordinary 
components are not distinguishable. 
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Ficure 8. Correlation between Nr and Ny. 


Nr: rate of occurrence of Es scatter, and Nw: rate of occurrence of Es affixed by 
symbol ‘“M”’, 


arranged into three groups according to the previous 
classification of ionospheric attenuation, namely: 
>70 db, 45 db <I' $70 db, and <45db. Arrange- 
ment was made, in this manner, of the data obtained 
during a period in July and August 1959, for which 
wide dynamic range recordings were available. Fig- 
ures 9 and 10 are the results of measurements. 

From these figures the following features are noted: 
The normal £-layer scatter of [>70 db showed a 
regular diurnal variation of angle of arrival with a 
narrow range of quartile distributions; on the con- 
trary, the Hs-reflected signal of db displayed 
a wide distribution of angles of arrival with no clear 
diurnal variation and a westerly deviation of bearing 
of several degrees on the average. On the other 
hand, the bearings of the Es scatter of 45 db <I' $70 
db exhibited a diurnal variation similar to that of 
the normal scatter, being accompanied with a wider 
distribution than that of the normal scatter. 

From the results mentioned above, it is inferred 
that the Hs evening enhancement or Es scatter is not 
caused by the normal reflection from sporadic-E 
clouds but by an irregularity of sporadic Z, and that 
the diurnal variation of bearing of the Es scatter 
may be related to the same mechanism as that of the 
normal #-layer scatter. 


5. F-Layer Scatter 


F-layer scatter was pointed out by R. Bateman 
and others [1] as an evening enhancement observed 
in the Far East in autumn. The event frequently 
occurred over the Poro-to-Tokyo path in the evening 
hours in September and October. When this took 
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Ficure 9. Bearing in each range of T <45 db, 45 db<T 
< 70 db andT > 70 db. 


Okuma-Tokyo, 49.68 Mc/s, July 22-August 21, 1959. 
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Fiaure 10. Angle of elevation in each range of T < 45 db, 
45 db <T <70 db andT >70 db. 


Okuma-Tokyo, 49.68 Mc/s, July 22-August 21, 1959. 


place, F-layer measurements at the Okinawa iono- 
sphere station quite near the path midpoint always 
indicated either of the following conditions: 

(a) F2 critical frequency higher than the 18 Me/s 
required for the normal reflection in oblique incidence 
propagation at the used frequency, or 

(b) F2-layer critical frequency qualified by the 
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symbol “R” or “S’’,* where the top frequency was 
higher than about 14 Me/s. 

Figure 11 is a typical example illustrating such a 
correspondence of the signal strength with the 
“¢-plot’’ of ionospheric data. There is no doubt that 
the signal received in autumn is closely associated 
with the F2-layer. It is also obvious that iono- 
speric layers other than the F2-layer did not take 
part in the propagation during the period of autumnal 
observation where especially sporadic / was inactive. 
Figure 12 shows diurnal curves of reception of signal, 
Nr, and N;, the occurrence of the two ionospheric 
conditions, ‘“R” and ‘‘S’’, observed in October 1959. 
Fair agreement can be seen between Nr and N,;. 
The spread-F condition was rarely observed at 
Okinawa which is quite near the path midpoint, 
during the observational period. 
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Fieure 11. Diurnal curves of signal strength and f-plot of 
ionospheric data observed at the Okinawa ionosphere 
station. 


Poro, the Philippines—Tokyo, 49.84 Mc/s, October 7, 1959. 


4R, Measurement influenced by, or impossible because of, absorption in the 
vicinity of a critical frequency. 
S, Measurement influenced by, or impossible because of, interference or 


atmospherics. 
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Comparison of Nr and Ny. 


Nr: rate of reception of F-layer scatter. 
Poro, the Philippines-Tokyo, 49.84 Me/s, October 1959. 
Nr: rate of occurrence of the ionospheric anomalies at Okinawa. 


FicureE 12. 


Pulse signal of 50 microseconds duration was 
observed at Tokyo on September 23, 1959. Figure 
13 exhibits the typical A-display of the received 
pulse. It is composed of a number of components 
corresponding to various propagation paths spread- 
ing over a range of more than 1 msec. The maxi- 


Ficure 13. A typical A-display of the received pulse. 


Poro, the Philippines-Tokyo, 49.84 Me/s, 2203 J.s.t., September 23, 1959, time 
mark of 1 msec, 
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mum delay time is of the same order as that ob- 
served by R. Bateman and others [1] for the Poro- 
to-Onna, Okinawa, path (1347 km). The iono- 
spheric condition when the picture was taken 
corresponded to the case (b) above. From these 
observations, it is inferred that the received signal 
is propagated by way of scattering from an anomaly 
in the F2-layver peculiar to a region in the low 
latitude in addition to the normal /2-propagation. 

Figure 14 shows the results observed in October 


1959. The signals generally arrive from the true 
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Fieure 14. Diurnal curves of signal strength, angles of 


arrival. 
Poro, the Philippines-Tokyo, 49.84 Mc/s, October 1959. 


bearing early in the evening, but their bearings 
gradually deviate westwards with the lapse of time. 
It is one of the future problems to clarify whether 
the westerly deviation is caused by movement of 
the anomaly or by such mechanism as scattering 
from field-alined blobs of ionization in the F region 
as suggested already [1]. 


The authors express their thanks to Dr. S. Namba, 
Chief of Research Laboratory, K.D.D., for his 
encouragement in the study, to CRPL and The 
Voice of America in the United States Government, 
and the Radio Research Laboratories, Japan, for 
their cooperation and assistance in the continuation 
of transmissions and in the data arrangement. 
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A High-Resolution Rapid-Scan Antenna’ 
H. V. Cottony and A. C. Wilson 


(August 11, 1960) 


An electronically scanned antenna array is described. 
of seven Yagi elements spaced 1.4 wavelengths apart. 
Each converter is connected to a local oscillator differing by 
20 c/s from that connected to the converter for the adjacent element. 
The converter outputs are connected to a common IF amplifier. 


to a preamplifier-converter. 


locked in phase. 


It consists of a broadside array 
Each element antenna is connected 


All oscillators are 
This 


arrangement produces a 5.8-degree-wide beam swept over a i gg azimuth sector at 


a rate of 20 scans per second. 


he system operates at 40.92 Me/s. 


he output of the sys- 


tem is displayed on an oscilloscope, and data on the direction of arrival and character of a 


distant VHF signal are presented visually. 


Sample records of signal components propagated by ionospheric scatter, meteor trail 


reflections, and sporadic-F layer are cited. 


1. Introduction 


An antenna having a very narrow width of main 
lobe is a useful tool in the field of propagation studies. 
It can enable an investigator to resolve the modes 
of propagation arriving from different directions. 
Many of the important propagation phenomena are, 
however, subject to rapid variation with time. It is 
important, therefore, that, in addition to high resolv- 
ing property, the antenna be capable of being scanned 
over the desired sector at a rate more rapid than the 
change of the propagation phenomena being studied. 
A project for the design and construction of an an- 
tenna system demonstrating this principle was au- 
thorized in July 1959. This paper describes the 
accomplishments on the project to date. 

The application of an antenna is essential in se- 
lecting its parameters. This particular antenna had 
as its object the measurement of the direction of 
arrival of a VHF signal (40.92 Mc/s) transmitted 
from a point 1295 km (805 mi) away, and propagated 
variously by ionospheric scatter, meteor trail reflec- 
tions, and sporadic E-layer reflections. It was an- 
ticipated that the direction of arrival might depart 
by 10 and more degrees off the great-circle route. 
It was not, however, anticipated to be further than 
about 20 deg off because the reflection points had to 
be visible optically from both the transmitting and 
receiving points. The signal scattered by the iono- 
sphere is known to be fluctuating rapidly. The du- 
ration of the meteor trail reflection varies, but is 
likely to be of the order of 1 or 2sec. For this reason 
the duration of the scan had to be of the order of 0.1 
secorshorter. In addition, the maximum acceptable 
rate of scan was, in this application, limited by con- 
siderations of pulse transmission and beamwidth. 


1 Contribution from Central Radio Propagation Laboratory, National Bureau 
of Standards, Boulder, Colo. 
2 Figures in brackets indicate the literature references at the end of this paper. 


In order to measure the path delay, it was planned 
to utilize pulse transmissions for some of the measure- 
ments. Earlier work by Carpenter and Ochs [1]? 
showed that relative path delays greater than 0.5 
msec were rare, but that reasonably numerous delays 
in the time of arrival up to that value could be ex- 
pected. This would restrict the pulse repetition rate 
to 2000 pps. A pulse repetition rate of 500 pps was 
decided upon. The pulse repetition rate and the 
beamwidth of the main lobe set a Jimit on the maxi- 
mum rate of scan. Thus, if the beamwidth were 1 
deg and the pulse repetition rate —500 pps, for com- 
plete sector coverage without gaps the maximum 
scan speed should be limited to 500 deg/sec. For a 
2-deg beamwidth it could be 1000 deg/sec. If the 
sector of scan were restricted to 40 deg, 20 deg to 
each side of the great-circle path, the maximum per- 
missibla scanning rate would, for the two cases, be 
12 and 25 scans per second. This is the limit on the 
maximum rate of scan. 

The conventional means of obtaining a radiation 
pattern narrow in one plane is to make an antenna 
wide in that plane. While there are available 
methods for reducing the beamwidth of an antenna 
of a given dimension, these methods have limitations. 
The array described here was designed and employed 
conventionally. It consists of 7 horizontally polar- 
ized 5-element Yagi antennas in a broadside array. 
The spacing was made as large as possible without 


introducing harmful directional ambiguities. The 
spacing settled upon was 1.4 wavelengths. This 


resulted in ambiguities at an angle of 20.9 degrees 
off the great-circle path. For any other point within 
the 20.9°-0°-20.9° sector, the other point of ambi- 
guity lies outside the sector and may be presumed to 
be very improbable. Five-element Yagis were 
selected because their half-power beamwidth in the 
E-plane, 56 deg, permitted operation within the 40- 
deg sector without introducing serious inequality of 
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Figure 1. A view of the antenna system. 


response. At the same time they were effective in 
discriminating against the signals outside the sector. 
The gain of the individual Yagi antenna was approxi- 
mately 11.3 db relative to an isotropic radiator. 
The beamwidth of the array, when directed straight 
ahead, was computed to be 5.8 deg for 20-db side- 
lobe suppression. 

Figure 1 is a view of the array. It is located on 
the Table Mesa near Boulder, Colo. It is directed 
at a bearing of 84.66° pointing along the great-circle 
path, towards the transmitter at Long Branch, IIl., 
field station 1295 km away. 


2. Principle of the Electronic Scan 


Because the method of electronic scan introduced 
with this antenna is the important contribution to 
the project, it is believed desirable to describe it in 
some detail. 

Figure 2 is a diagram of a 7-element stationary 
array in which the directivity is obtained by proper 
phasing of the element currents. If the elements 
are equispaced with a separation d and the beam is 
to be directed at an angle 6 to the right, it is necessary 
to advance the phase of the current in elements 3,, 
2,, and 1, by phase angles 36, 26, and 6 where 
6=27(d/d) sin 8. The phases of the currents in the 
elements on the right side have to be retarded by 
corresponding angles. Such phasing can be accom- 
plished either mechanically, electrically, or electron- 
ically. Many ingenious methods of doing this have 
been tried [2]. All of them have some drawbacks. 
With this antenna, a method of phasing in a mixer 
circuit is introduced. 

It is known that when two voltages, ¢,=A sin 
2 af, and ¢.=B sin (22f,+6), are connected to the 
inputs of a nonilnear network, the outputs will con- 


FORWARD 
DIRECTION —Dyp 
EC, 


7, 
| MAly tone 


/ 


£0y, 
/ 
/ 
| 


Figure 2, A diagram of a seven-element broadside array 
illustrating the steering of the main beam by phasing. 


tain the sum and difference of frequency voltages, 
es=C cos [24(fit+fe) +4] and eg=C cos [22(f;—f2) — 4], 
in which the phase relationship of the original volt- 
ages is not lost. If e is the local oscillator voltage, 
the phase of the signal contributed by an element 
antenna is controllable by varying the phase of the 
local oscillator. The’ practical means of doing it is 
by varying 6 continuously in one direction so that in 
a time interval, ¢,, it changes from 0 to 27. Such a 
continuous phase change is equivalent to an incre- 
mented frequency change, Af=(1/t,). To obtain 
proper scan of the antenna, it is necessary to connect 
each to its own converter with the local heterodyne 
voltages differing in frequency by frequency incre- 
ment Af. If the local heterodyne voltage for antenna 
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0 is =A sin 2zf, the local heterodyne voltages for 
antennas 3,, 2;, and 1, will be es, and e,, equal 
respectively to Asin A sin 2r(fo—2A/) and 
A sin 2x(f.—Af). The element antennas on the right 
side of the array should have local heterodyne frequen- 
cies increased by the same increments. The beam of 
such antenna will rotate unidirectionally from left 
to right with a scan frequency equal to af. It 
should be noted that the range of scan is from one 
main lobe to another. If the element spacing, d, 
were equal to \/2, for 6= +7 there would be two main 
lobes, at B=—90° and +909°, respectively. Thus, 
the scan would begin at B=—90° and sweep to 
+90°. If the spacing, d, were less than /2, the main 
beam would form in the imaginery region of sin 6 
<—1, sweep through the real sector —90° to 0° to 
+90°, and disappear into the imaginary region of sin 8 
>+1. If d were greater than /2, several main 
beams might be present at all times and the entire 
ensemble of main beams would rotate clockwise. 
For d=1.4X and 6=z, there is a main Jobe in the 
imaginary region, B=sin~'—1.07, two real lobes at 
—20.9° and +20.9°, respectively, and another lobe 
in the imaginary region, B=sin™'+1.07. In one 
cycle the beams sweep through one interval between 
lobes. In this discussion, the fact that the radiation 
pattern has mirror-image symmetry about the line 
of the array has been ignored. In physical con- 
struction, a reflecting screen or unidirectional ele- 
ments would normally be employed. 

The proper operation of this scanning technique 
is very dependent upon finding a successful solution 
to the problem of generating an ensemble of local 
heterodyne voltages. These have to differ in 
frequency by the increment Af, but must be locked 
in phase. It is possible to accomplish this in a 
number of ways. The method selected here was to 
pulse the oscillator having a fundamental frequency 
fo at a pulse rate corresponding to Af. If the length 
of the pulse is short, a spectrum of frequencies is 
generated with the center frequency at fy and other 
frequencies at increments of Af up and down the 
frequency scale, all of these voltages will be in phase 
or 180 deg out of phase. If the length of the pulse 
is short compared to 2/(naf), the amplitudes of 
different frequency voltages will be approximately 
the same and the phases the same. 1 is the number 
of elements in the array. The frequencies can then 
be separated by crystal filters. In practice, it was 
found desirable to generate and separate the fre- 
quency ensemble in the 20 ke/s band, and to step up 
the frequencies after separation, by heterodyning 
against a common omlanes in two steps. The 
20-c/s frequency separation produced a corresponding 
sweep rate. This was rapid enough to observe the 
propagational variations with time. It was also 
sufficiently slow to observe pulse transmission 
without gaps in the sector coverage. 


3. Details of Construction 


The design for the five-element Yagi antenna had 
been worked out in detail a number of years pre- 
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viously for another application. Its gain was 
measured to be 11.3 db relative to an isotropic 
radiator. The approximate half-power width of the 
main lobe was measured to be 50 deg. Within the 
sector of interest, 21°-0°-21°, the gain varied 
within +1 db of the mean. The secondary lobe 
level was below —24 db. Each Yagi was matched 
to have a coaxial terminal impedance of 50.0 ohms. 

The Yagi antennas were separated by 33.65 ft 
(1.4 wavelengths). From earlier measurements, it 
is known that, at this separation, the effect of one 
of these antennas on the radiation pattern of the 
adjacent one is negligible. 

For communication over a 1295 km path via 
ionospheric scatter mode, the optimum angle of 
departure is computed to be approximately 6.7 deg 
and the height of the antenna, 2.15 wavelengths 
(approximately 52 ft). All of the antennas were 
connected to the equipment by equal lengths of 
RG-9B/U coaxial cable. It was contemplated that 
the antenna would be adjusted to have a Dolph- 
Chebyshev current distribution for a —20-db side- 
lobe level. On that basis, the beamwidth of the 
array was computed to be 5.84 deg. 

Figure 3 presents the computed radiation pattern 
of the array. The curve in dashed line is the radia- 
tion pattern of the single Yagi antenna. In com- 
puting this pattern, the antennas were assumed to be 
in phase; this produced the main lobe in the straight- 
ahead position. Other main lobes are at the 45°, 
135°, and 180° positions. 

The frequency of 10.70 Me/s was selected for IF 
amplification; therefore, for 40.92-Mc/s reception, 
the local oscillator frequencies had to be centered on 
30.22 Me/s. With the present state of filtering 
techniques, it would be difficult to separate com- 
ponent voltages spaced 20 c/s apart at this frequency. 
To avoid lengthy experimentation, the frequency 
of 18 ke/s was selected as being the highest at which 
ready separation of frequencies could be carried out. 


Figure 4 presents a block diagram of the electronic 
circuitry involved in the production of beam scan. 
An 18-ke/s crystal-controlled oscillator is pulsed at 
a rate of Af=20 pps. The length of the pulse is 
4 msec; a short pulse is required in order that the 
amplitudes of the required sidebands remain in 
phase and reasonably constant. The output of the 
pulsed oscillator contains frequencies of 18 ke/s, 
18 ke/s+af, 18 ke/s+2af,...18 ke/s—af, 18 
ke/s--2af, . . . etc. The output is connected to 
seven filters in parallel, each adjusted to one of the 
desired frequencies. The passband of each crystal 
filter is essentially flat for 6 cycles, centered at the 
individual crystal filter design frequency. The 
rejection of the unwanted sidebands by each crystal 
filter is at least 30 db. Each selected frequency, 
18 ke/s+raf, is stepped up by beating against a 
2-Me/s oscillator to 1.982 Me/s—raf. This fre- 
quency, in turn, is filtered from 2 Me/s and 2.018 
Me/s+-raf frequencies. 

The 1.982 Mc/s—raf frequency is again stepped 
up to 30.22 Me/s—raf by heterodyning with a 28.238 
Mes oscillator, and separated from 28.238 Me/s and 


| 
| 


RESPONSE. DECIBELS 


WORMALITED 


FIGuRE 3. 


AZIMUTHAL ANGLE, DEGREES 


EXPERIMENTAL PATTERN 
V7 OF FIVE ELEMENT YAGI 


\ 
4 
\ 


\ 
| 
\ 
i] 
\ 
! 
| 


Computed radiation pattern in E-plane of an array of 7 Yagi antennas. 


Currents adjusted to Dolph-Chebyshev distribution for -20-db side-lobe level. 


INPUT SIGNAL FROM YAGI ANTENNA, 


40.92 Mc/s 
| CRYSTAL con. |49:92 Mc/s | 
| FILTER MIXER | FILTER MIXER | FILTER Bo | 
18 kc/s + rAf 18 ke/s 1.982 Mc/s 30.22 Mc/s AMPLIFIER 
| + rAf - - rAf | | 
18 kc/s +0, +Af, +2Af,***: 
-Af, -2Af,----- 2 Mc/s 28.238 Mc/s 
10.7 Mc/s 
OUTPUT TO IF 
AMPLIFIER 
PULSE GENERATOR 18 kc/s 2.0 Mc/s 28.238 Mc/s 
20 pps S. CRYSTAL- CRYSTAL- CRYSTAL- 
4 msec PULSE CONTROLLED CONTROLLED CONTROLLED 
WIDTH OSCILLATOR OSCILLATOR OSCILLATOR 
Figure 4. Block diagram of one filter-mizer and preamplifier-converter channel 


26.256 Me/s+raf frequencies by filtering. From 
there it is connected to its appropriate, rth, converter. 
The seven crystal filters for the 18-ke/s frequency 
ensemble are the only critically adjusted items of 
the circuitry. Other filters and components are 


readily interchangeable with components of other | 


channels because frequency separation is reasonable. 
In order not to introduce phase differences, care was 
exercised in adjusting antennas to have identical 
characteristics, to cut all transmission lines to be 
equal in length, and to construct all equipment per- 
taining to the seven channels to be as nearly identical 
electrically as possible. This avoided gross differ- 
ences in phases or amplitudes. Provisions were 
incorporated in the frequency-generating circuits for 
compensating for small unavoidable discrepancies in 
phase. Provisions were also incorporated for adjust- 
ments in amplitude over a range of several decibels. 
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This was to provide for tapering the current distri- 
bution of the array. Figure 5 shows a view of the 
equipment mounted in the rack. 


4. Adjustment of the Antenna 


Before the antenna was placed in operation, it was 
expected that difficulties would be encountered with 
unequal phase and amplitude variations in the seven 
channels. It was also felt that the adjustments of 
the electronic equipment could not be expected to 
remain sufficiently stable for very long periods of 
time. Practical experience with the equipment, 
however, dispelled these apprehensions. Phase and 


amplitude differences were small and easily compen- 
sated for; the stability was such that adjustments 
remained satisfactory from one day to another, even 
The process of adjust- 


after equipment shutdown. 
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Figure 5. 


A view of the electronic equipment for production of 


rapid scan, 


ment itself proved to be easier than with a non- 
scanning array because the effect of each adjustment 
was visible immediately on the oscilloscope display. 

After some experimentation, a procedure for ad- 
justment was evolved. The phases and amplitudes 
of each channel were adjusted, one at a time, with 
reference to the center p aie as follows: A preci- 
sion coaxial attenuator with insertion loss equal to 
the level of the rth element for Dolph-Chebyshev 
distribution is inserted in the coaxial line leading to 
the Oth (center) element. The rth element is con- 
nected through its circuitry with an insertion of 180- 
deg delay (half-wavelength coaxial line). The 
outputs of the Oth and rth converters are then 
combined in the IF amplifier and displayed on the 
oscilloscope. The amplitude of the rth element is 
adjusted to obtain a null (or a series of nulls, depend- 
ing on the value of r). The phase is adjusted to 
move the null into the center of the display. This 
is the complete adjustment routine for one, rth, 
element. After all elements are adjusted, they are 
connected together through their respective convert- 
ers to the IF amplifier. It is desirable to check the 
pattern by a distant target transmitter. Minor ad- 
cnngee can then be made to improve the secondary 
obe level. A mobile target transmitter was also 
found to provide a convincing demonstration that 
the equipment really operated as intended. 


In lieu of a distant target transmitter a signal 
generator can be, and has been, employed. The 
antennas are disconnected. Each pre-amplifier is 
then connected through an individual, 20-db, preci- 
sion, coaxial attenuator to the signal generator. The 
distant target transmitter is thus simulated in the 
zero-degree position. If the combined input imped- 
ance of the seven circuits in parallel is matched: 
back to the 50-ohm value, a means of quantitative 
calibration is provided. The equipment described 
is provided with solenoidally operated switches to 
provide means for rapid check and calibration. 


5. Results 


Since the output of the antenna is cyclic with a 
20 c/s repetition rate, it lends itself readily to oscillo- 
scope presentation. The horizontal sweep is trig- 
gered by the pulse, and the output of the receiver 
is connected to the vertical deflection plates. The 
display shows the signal amplitude as ordinate, with 
the direction of arrival as azimuth. Since a distant 
VHF signal arriving by ionospheric scatter mode is 
distributed in the direction of arrival, the display 
pattern is a convolution of the antenna radiation 
pattern and the true signal or signal pattern distri- 
bution. The phenomenon is familiar in radio 
astronomy and has been discussed at length by 
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Bracewell and Roberts [3]. When the signal arrives 
from a point source, as from a target transmitter, 
the display shows the true radiation pattern with 
the response maximum in the direction of the signal 
source. Figure 6 presents three displays recorded 
while a local target transmitter, mounted on a 
vehicle, moved in front of the antenna from left to 
right. It is useful for demonstrating the proper 
functioning of the system. 

Since this system was put in operation, 2 months 
ago, a variety of propagational phenomena have been 
observed on the frequencies of 40.88 and 40.92 Me/s. 
During the greater proportion of the time, the 
observed signal arrives by a regular ionospheric 
scatter mode. Occasionally, strong meteor-trail 
reflections have been observed and recorded. These 
exceed the scattered signal by several orders in 
magnitude, but last only a second or two at a time. 
Signals propagated by sporadic-F layer have also 
been observed. These are very steady, high in 
amplitude, and usually arrive along the great-circle 
path. Figures 7 through 10 illustrate the phenomena 
as observed using this antenna system. It should 
be noted that the recording camera was not syn- 
chronized with the sweep of the oscilloscope. In 
order to get at least one complete sweep, the shutter 
speed was set at 0.1 sec. Therefore, most records 
show two almost complete scans. These are not 
believed to be objectionable; in fact, they convey 
the transient nature of the phenomena. The interval 
between sweeps is 50 msec. 

Figure 7 presents photographic records illustrating 
propagation by ionospheric scatter. While the 
maximum signal is generally at, or close to, the 
great-circle path (0-deg azimuth), it fluctuates, 
both in amplitude and in direction. Occasionally, 
it deviates by more than 10 deg. from the great- 
circle path. The convolution of the radiation and 
the signal patterns is noticeably broader than the 
radiation pattern by itself. This is clear evidence 
that the signal pattern is dispersed. A_ higher 
resolution antenna jis highly desirable for more 
detailed examination of the fine structure of the 
signal. 

Figure 8 is a record of several meteors. Although 
some are observed at zero azimuth, most are to one 
side or the other of the great-circle path. Because 
the duration of the meteor is very brief, it is not 
possible to adjust the sensitivity of the receivers to 
obtain a good full-scale deflection. Experiments are 
underway to adopt magnetic tape recording tech- 
niques for this application. The display pattern is 
very close to the radiation pattern of the antenna, 
indicating that the signal comes from a source 
which is virtually a point. 

Figure 9 shows a record observed while sporadic 
E-layer propagation was taking place. The display 
pattern appears to be identical to the radiation pat- 
tern. This is to be expected from a_ specular 
reflection by a smooth surface. The signal is likely 
to remain steady for a considerable length of time and 
is usually along the great-circle path. 
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Ficure 6. Display of responses obtained with a target trans- 
mitter in far zone moving from left to right. 


The response to a point source of radiation is the radiation pattern of the antenna. 
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Ficure 7. Display of responses with signal arriving via 
normal ionospheric scatter mode. 


The response is a convolution of the antenna radiation pattern and the angular 
distribution of the signal components. 
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Figure 8. Display of responses with signal arriving via 
meteor trail reflection. 
The propagation mode identifiable by enhanced signal strength, duration of 


the order of a second, and off-great-circle direction of arrival. 
The similarity between the response pattern and the antenna radiation pattern 


indicates that the reflection is a point source. 
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Fiaure 9. Display of responses with signal arriving via 
sporadic-E layer reflection. 


The sensitivity of the system is greatly reduced relative to that used with 
scatter reception (of the order of 90 db). 

Besides high intensity, this mode is characterized by great-circle path direction 
of arrival and stability over a period of an hour or greater. 


Simultaneously with the 20-kw signal at 40.92 
Me/s, used for the above experiments, there was a 
2-kw transmission at 40.88 Mc/s. Both originated 
at Long Branch, Ill. The 40.92 Me/s signal was 
radiated by a corner-reflector antenna with a half- 
power beamwidth of approximately 50 deg. The 
40.88 Me/s signal was radiated by a rhombic antenna 
with a beamwidth of about 5 deg. Because the area 
of the ionosphere illuminated by one was so much 
greater than that illuminated by the other, it was con- 
sidered desirable to compare signals received simulta- 
neously on two frequencies. Two IF receivers tuned 
to10.70 and 10.66 Mc/s were connected to the output 
of the antenna system at the same time. The 
outputs of the two receivers were connected to 
sections A and B of a dual-trace oscilloscope. 
Figure 10 illustrates the comparison between the 
two signals. While the study of propagational 
phenomena is beyond the scope of this paper, the 
application illustrates the versatility of the scanning 
principle. The principal limit on bandwidth of 
operation was, in this case, presented by the Yagi 
antennas comprising the antenna system. 

In the earlier discussion of the principles of the 
scan, it was mentioned that measurements using 
pulse transmissions were contemplated. In 1958 
and 1959, Carpenter and Ochs [1], using pulsed 
40.92 Me/s transmissions, with 3-ysec pulse, at 250 
pps and 700 kw peak power, measured relative 
delay time for meteor reflections. Delays of over 
1 msec relative to the normal scatter signal were 
observed. By recording the distribution, they were 
able to deduce the direction of arrival of reflected 
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Figure 10. Display of simultaneous responses to signal at 
40.88 and 40.92 Mc /|s. 


In addition to frequency difference, the 40.88-Mc/s transmission was character- 
ized by narrow (5-deg) beamwidth. 40.92-Me/s transmission used 50-deg-wide 
antenna beam; 10-db difference in transmitted power maintained intensity of 
illumination at approximately equal level. 


signals. The equipment employed by Carpenter 
and Ochs has been, to a large extent, incorporated 
in this project. It is planned to modify the display 
so as to present the direction of arrival as azimuth 
(as it is now) but to show the relative delay in 
microseconds as ordinate and with intensity modula- 
tion amplitude. The display is similar to that 
familiar in television. Figure 11 shows the raster 
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Figure 11. Scan ‘“‘Raster’’ for reception of pulses. 


Two-millisecond scan corresponds to pulsing period (500 pps). 
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form of the display. It should be noted that the 
line separation is 2 deg; this illustrates the discussion 
of the interrelations of minimum beamwidth and 
the maximum scanning rate earlier in this paper. 
Figure 12 presents a record of a signal arriving by a 
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FicuRE 12, Display of response to pulse transmission using 
intensity modulation of oscilloscope and vertical sweep 
sychronized with pulsing rate. 


sporadic-E mode. Because of lack of power in the 
pulse transmitter, reception of regular ionospheric 
scatter has been unsatisfactory. 

Work is under way to overhaul the transmitter 
for high-power pulse transmission. Work is also 
under way to increase the size of the array to 25 
elements. This would permit a beamwidth of 1.5 
deg with the present 20-db side-lobe level. A more 
detailed study of propagation phenomena can then 
be possible. 


6. Applications 


The significant characteristic of this antenna is 
the ease of obtaining very rapid rates of scan. The 
20 c/s rate was the slowest which was deemed readily 
realizable, without experimentation. As the scan 
rate is increased, the problem of separating the 
spectral frequencies of the pulse becomes easier. 
Scanning rates in the order of kilocycles, and even 
megacycles, appear realizable. 

The availability of very rapid scanning rates 
opens up an opportunity for use of sampling tech- 
nique for both measurement and communication. 
If the period of scan is shorter than the Nyquist 
interval for the message received, the information 
is fully preserved, even though the antenna is 
directed at the source only a fraction of the time. 
By proper gating, the information arriving from 
different sources can be directed into appropriate 


channels. The technique is analagous to that 
familiar in other applications. 

The antenna described in this paper scans in a 
single dimension. With a minor modification, the 
principle is equally applicable to a two-dimensional 
scan. Figure 13 illustrates this application. A 


Figure 13. Application of electronic sweep to scanning in 
two dimensions, 


rectangular array of mn elements is shown. 
These are arranged in m rows and n columns, 
n{(m—1)/2] elements to each side of the center row, 
and m[(n—1)/2] elements to each side of the center 
column. The elements may be identified by their 
position in the array, by row and column. 

The center element (0,0) is connected to a con- 
verter with a heterodyne frequency fo. The array 
can then be scanned in the 6-plane at a scanning 
rate Af if each converter receives its proper hetero- 
dyne voltage frequencies. For element (7,q), it 
should be f,,~=fotgA4f. To obtain scan in the 
¢-plane, at a scanning rate df, the heterodyne voltage 
frequency for the same element should be f,,=fo+ 
réf. To scan in both planes simultaneously, f,~= 
fytréftqAf. The scanning operation itself would 
be analogous to the scan in the television picture. 
The ratio Af/5f would determine the number of “‘lines”’ 
in the scan. It should be equal to, or greater than, 
the ratio of 6/W¢ where ® is the width of the sector 
scanned in the ¢-plane and Ws is the half-power 
beamwidth in that plane. Since Ws varies roughly 
inversely as the aperture, (m—1)d, the ratio Af/éf 
should equal approximately the number of rows, m. 

A special case of interest above is a possible 
application of Mills’ cross array. In this applica- 
tion, the rectangular array reduces to the central 
row and the central column with a common center 
element (0,0). The multiplicative process essential 
to the operation of Mills’ cross can be added in 
addition to the scanning. A complete picture of a 
phenomenon varying with time over an area can 
then be obtained instantaneously. This might be of 
interest, for example, in radio astronomy in the 
study of solar eruptions. 
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All of the above applications are from the stand- 
point of receiving antennas. The principle is ap- 
plicable to transmission, but the techniques have 
to be modified considerably. 


7. Conclusions 


The 2 months of operation of this antenna system 
provide adequate evidence of the practical value of 
the scanning principle described here. In the par- 
ticular application detailed here, it provided a visual 
view of rapid variations in amplitude and direction 
of arrival with time. The information value of such 
a display is considerably greater than that of single- 
dimensional variation of amplitude with time. 

The stability of the system with time, temperature 
changes, etc., appears to be satisfactory. It should 
be noted that, in many respects, the equipment as 
it is at present may be termed “breadboard.” With 
some refinements, a very high order of dependability 
should be realizable. The adjustment of the array 
is actually much easier than that of a conventional 
array without scan. Since the pattern can be dis- 
played on an oscilloscope, the effect of changing 
individual adjustment controls is instantly visible. 
The process of adjustment, therefore, becomes greatly 
expedited. 

The 5.8 deg beamwidth of the seven-element, 8.4 
wavelength-wide array is barely adequate to demon- 


strate the possibilities of the system. Reduction in 
beamwidth is, however, a matter of adding to the 
array. Work is under way to increase the number 
of elements to 25, and the overall array width to 
33.6 wavelengths. The beamwidth of the main lobe 
a = be reduced to 1.5 deg with 20-db side-lobe 
evel. 

The generation of a phase-locked ensemble of fre- 
quencies for heterodyne use may be accomplished in 
a number of ways other than that described here. 
Phase-lock oscillators can be employed. A single 
oscillator at a frequency Af may be used to generate 
harmonics which can be amplified and stepped up in 
frequency by beating against a common oscillator. 

Application to two-dimensional scan appears to be 
straightforward but, of course, remains to be tried. 
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Amplitude and phase of the low- and very-low-radiofrequency 
ground wave, J. R. Johler, L. C. Walters, and C. M. Lilley, 
NBS Tech. Note 60 (PB161561) (1960) $0.75. 


Graphs and tables of the low- and very-low-radiofrequency 
ground wave are presented as a function of frequency, 
100 c/s to 1000 ke. 


On the excitation of electromagnetic surface waves on a 
curved surface, J. R. Wait, JRE Trans. Ant. Prop. AP-8, 
446 (1960). 


The excitation and propagation of surface waves on a spheri- 
eal inductive boundary 1s considered. The source is taken 
to be a vertical elecuric dipole. The circumferential attenu- 
ation rate of the various modes are discussed where it is 
indicated that the dominant mode is very similar to the 
trapped surface wave for a plane inductive boundary. The 
results appear to conflict with those of Barlow, but are in 
sympathy with some numerical data of Elliott for the cir- 
cumferential attenuation rate of the dominant mode. 


Carrier-frequency dependence of the basic transmission loss 
in tropospheric forward scatter propagation, K. A. Norton, 
J. Geophys. Research 65, 2029 (1960). 


A further interpretation is given of certain Lincoln Laboratory 
data obtained in} an experiment using scaled antennas as 
published by Bolgiano [1959]. This paper has three objec- 
tives: (1) to consider the significance of these data as regards 
the theory of radio propagation through a turbulent atmos- 
phere; (2) io describe a suitable method for the measurement 
of the meteorological parameters entering the theory; (3) to 
apply a further statistical analysis to these data. On the 
basis of this analysis of the Lincoln Laboratory data, it is 
concluded (1) that the 


(u) (p/2)*Ku(p) 


correlation model provides a description of tropospheric 
turbulence adequate for the description of these data with 
uw assigned the fixed value 1, and (2) that the variation of 
the carrier-frequency dependence of the basic transmission 
loss from hour to hour arises simply from a lack of correlation 
in the hourly median losses on widely separated carrier fre- 
quencies rather than from changes of uw in this correlation 
model. 


The use of geostationary satellites for the study of ionospheric 
electron content and ionospheric radio-wave propagation, 
O. K. Garriott and C. G,. Little, J. Geophys. Research 65, 
2025 (1960). 


This paper describes a proposed radio propagation experi- 
ment utilizing a geostationary satellite. The experimental 
methods for determining the total electron content of the 
ionosphere are discussed, and the advantages and capabilities 
of the methods are summarized. Other possible studies 
utilizing the satellite transmissions are also outlined. 


Photochemical rates in the equatorial F, region from the 
1958 eclipse, T. E. Van Zandt, R. B. Norton, and G. H. 
Stonehocker, J. Geophys. Research 65, 2003 (1960). 


The F,-region electron-density data from the eclipse on 
October 12, 1958, at Danger Islands has been analyzed with 
anew method. It is shown that (1) temperature changes and 
transport of electrons were probably negligible during the 
eclipse, and (2) the rate of photoionization qg(h) and the 
linear recombination coefficient (‘attachment coefficient’) 
B(h) can be approximated between 290 km and 400 km by 


q(h) =880 exp [— (h—300) /186] electrons 
(h)=6.8X 10-4 exp [— (h—300)/103] 


where the height h is measured in kilometers. The rate of 
formation of ion pairs between 290 and 400 km is (.76)10! 
ion pairs/cm? sec, or about (0.25) ergs/em? sec. This implies 
a total rate of energy absorption in the F region of at least 
2 ergs/cm? sec. 


On the theory of the slow-tail portion of atmospheric wave- 
forms, J. R. Wait, J. Geophys. Research 65, 1939 (1960). 


The propagation of the slow-tail portion of atmospherics is 
considered from the waveguide-mode viewpoint. The 
source, which is a lightning discharge, is represented by a 
vertical dipole. The transient response of the distant electric 
field is then computed for various forms of the source current 
waveform. The results are then employed to reinterpret 
the experimental data of Hepburn. As suggested by the 
present theory it is found that the observed separation 
t, between the oscillatory head of the atmospheric and the 
maximum of the slow-tail amplitude varies with distance 
p to the source according to a law of the form 


1/2=A-+ Bp. 


The constant A is related to the pulse width of the source, 
and the constant B depends on ionospheric parameters. 
Values of effective ionospheric conductivities deduced from 
the theory are consistent with earlier results for the VLF 
band. The influence of nonvertical currents in the discharge 
channel is also briefly discussed. 


VLF attenuation for east-west and west-east daytime propa- 
gation using atmospherics, W. L. Taylor, /. Geophys. Re- 
search 65, 1933 (1960). 


Daytime attenuation characteristics have been computed 
by comparing the amplitude spectra of atmospheric wave- 
forms recorded at four widely separated stations. The 
results of these attenuation measurements are presented 
for the band of frequencies from 3 to 30 ke/s and involving 
distances of 1000 to 10,000 km. Nonreciprocity is evident 
from this study. Attenuation rates over sea water for east- 
to-west propagation were about 3 db/1000 km greater for 
frequencies below 8 ke/s, and about 1 db/1000 km greater 
for frequencies above 10 ke/s, than for west-to-east propa- 
gation. East-to-west attenuation rates over land were 
about 1 db/1000 greater than for over sea water. 


Influence of earth curvature and the terrestrial magnetic 
field on VLF propagation, J. R. Wait and K. Spies, J. Geophys. 
Research 65, 2325 (1960). 


An account is given of some recent work on the mode theory 
of VLF ionospheric propagation. Attention is confined to 
the behavior of the attenuation coefficient of the dominant 
mode. The ionosphere is assumed to be a sharply bounded 
and homogeneous ionized medium. It is indicated that 
earth curvature increases the attenuation rate by as much 
as a factor of 2 as compared with the corresponding atten- 
uation for a flat earth. The influence of the earth’s magnetic 
field is also shown to be important. In fact, east-to-west 
propagation paths suffer much greater attenuation than 
west-to-east paths. The theoretical results in the present 

aper appear to agree well with the experimental data of 
W. L. Taylor. 


Rapid frequency analysis of fading radio signals, J. M. Watts 
and K. Davies, J. Geophys. Research 65, 2295 (1960). 


Examples of frequency analysis of fading radio signals for 
long periods of time are demonstrated, and the method of 
obtaining them is explained. They include both regular HF 
propagation and VHF ionospheric forward-scatter samples. 
The procedure is also useful for the analysis of other natural 
phenomena having long time scales and slow variations. 


111 


| 
| 

= 


A study of 2-Mc’/s ionospheric absorption measurements at 
high latitudes, K. Davies, J. Geophys. Research 65, 2285 
(1960). 


Ionospheric absorption (L) at high latitudes is studied using 
the published data on 2.0 Me/s at five Canadian stations 
during 1957 and 1958. The seasonal and diurnal variations 
are considered, and it is found that a pronounced winter 
anomaly in noon absorption occurs at Churchill but not 
at Resolute Bay. The diurnal variations indicate that the 
dependence of absorption on solar zenith angle decreases 
with increase of latitude. The distribution of midnight 
absorption with latitude shows that, although the maximum 
occurs in the auroral zone, high absorption is also encoun- 
tered over the polar cap. A study of the duration of long- 
lasting blackouts shows that in summer the duration is 
longer as the latitude increases. 


On the excitation of electromagnetic surface waves on a 
curved surface, J. R. Wait, TRE Trans. Ant. Prop. AP-8, 
445 (1960). 


The excitation and propagation of surface waves on a spherical 
inductive boundary are considered. The source is taken to 
be a vertical electric dipole. The circumferential attenuation 
rates of the various modes are discussed where it is indicated 
that the dominant mode is very similar to the trapped 
surface wave for a plane inductive boundary. The results 
appear to conflict with those of Barlow, but are in sympathy 
with some numerical data of Elliott for the circumferential 
attenuation rate of the dominant mode. 


Diffractive corrections to the geometrical optics of low 
frequency propagation, J. R. Wait, Electromagnetic Wave 
Propagation (International Conference Sponsored by _ the 
Postal and Telecommunications Group of the Brusseis Uni- 
versal Exhibition), edited by M. Desirant and J. L. Michiels, 
p. 87-101 (Academic Press, New York, N.Y., 1960). 


The radiation fields are computed for an axial slot on a 
evlinder whose surface is regularly corrugated. The mathe- 
matical idealization is a magnetic line source on a cylindrical 
surface whose boundary impedance is specified. Special 
attention is given to surfaces whose radii of curvature are 
large compared to the wave-length. 


Radiation from a slot on a large corrugated cylinder, J. R. 
Wait and A. M. Conda, Electromagnetic Wave Propagation 
(International Conference Sponsored by the Postal and Tele- 
communications Group of the Brussels Universal Exhibition), 
edited by M. Desirant and J. L. Michiels, p, 103-109 (Aca- 
demic Press, New York, N.Y., 1960). 


The radiation fields are computed for an axial slot on a 
cylinder whose surface is regularly corrugated. The mathe- 
matical idealization is a magnetic line source on a cylindrical 
surface whose boundary impedance is specified. Special 
attention is given to surfaces whose radii of curvature are 
large compared to the wave-length. 


Atmospheric bending of radio waves, B. R. Bean, Electromag- 
netic Wave Propagation (International Conference Sponsored 
by the Postal and Telecommunications Group of the Brussels 
Universal Exhibition), edited by M. Desirant and J. L. Michiels, 
p. 163-181 (Academic Press, New York, N.Y., 1960). 


Meteorological data from many diverse climates and seasons 
indicate that the refractive index gradient over the first kilo- 
meter above the earth’s surface is highly correlated with the 
surface value of the refractive index. Further, 60% of the 
bending of radio rays passing completely through the earth’s 
atmosphere is accomplished in the first kilometer above the 
earth’s surface. These observations have led to the develop- 
ment of several models of the height distribution of the refrac- 
tive index that yield more realistic values of radio ray bending 
than the effective earth’s radius approach and have the fur- 
ther advantage of being adjustable for any season or location 
for which standard surface weather observations are available. 


Low and medium frequency radio propagation, K. A. Norton, 
Electromagnetic Wave Propagation (International Conference 
Sponsored by the Postal and Telecommunications Group of the 
Brussels Universal Exhibition), edited by M. Desirant and J. L. 
Michiels, p. 375-444 (Academic Press, New York, N.Y., 1960). 


During recent years, extensive progress has been made in un- 
derstanding the nature of ionospheric propagation in the fre- 
quency range from 30 to 1000 ke. In particular, Wait and 
Conda have developed suitable methods for determining the 
effect of the curvature of the earth on the illumination of the 
ionosphere by antennas radiating at angles near and below 
grazing incidence, Belrose has developed empirical formulas 
for ionospheric absorption in the range 70 to 250 ke, Watt has 
analyzed data which makes possible the extension of these 
formulas to still lower frequencies, and Bean has developed a 
radio standard tropospheric atmosphere which makes pos- 
sible good estimates of the bending of the radio waves in 
propagation to and from the ionosphere. Making use of these 
results, together with earlier studies made by the author of 
the absorption of radio waves in the band 540 to 1600 ke and 
some new results on focusing by the curved surface of a rough 
ionosphere, predictions are made of the propagation loss ex- 
pected in ionospheric propagation between short vertical 
electric dipoles for waves propagated by m reflections at the 
ionosphere. The predictions appear to be in general agree- 
ment with the available experimental data, although the 
paper clearly points up the desirability of more definitive 
studies of the influence of absorption and polarization, par- 
ticularly as regards the larger distances which involve propa- 
gation by more than one ionospheric reflection. The meth- 
ods in this paper indicate that the propagation loss increases 
only very slowly as the angle of departure approaches and 
goes well below grazing incidence with the earth; this pro- 
vides an explanacion for recent experimental data obtained 
with a pulse system as reported by Doherty which indicates 
very large intensities for one-hop sky waves on 100 ke both 
day and night at ranges up to 1800 miles, 


On the computation of diffraction fields for grazing angles, 
J. R. Wait and A. M. Conda, Electromagnetic Wave Propaga- 
tion (International Conference Sponsored by the Postal and 
Telecommunications Group of the Brussels Universal Exhibi- 
tion), edited by M. Desirant and J. L. Michiels, p. 661-670 
(Academic Press, New York, N.Y., 1960). 


The diffraction of electromagnetic waves by a convex cylindri- 
cal surface is considered. Attention is confined primarily to 
the region near the light-shadow boundary. The complex- 
integral representation for the field is utilized to obtain a cor- 
rection to the Kirchoff theory. Numerical results are pre- 
sented which illustrate the influence of surface curvature and 
polarization on the diffraction pattern. Good agreement 
with the experimental results of Bachynski and Neugebauer 
is obtained. 


The resonance excitation of a corrugated-cylinder antenna, 
J. R. Wait and A. M, Conda. Inst. Elec. Engrs. Mono, 386E 


(June 1960). 


Radiation from an axial magnetic line or slot source on the 
surface of a corrugated cylinder is considered. It is indicated 
that the power radiated in a given mode for the structure de- 
pends critically on the surface reactance and the circumfer- 
ence of the cylinder. In fact, for certain values of these 
parameters, particular modes are strongly excited and con- 
tain most of the radiated power. Numerical results are pre- 
sented for several interesting cases. 

The analysis is extended to an elliptic cylinder whose surface 
also possesses an inductive reactance. In order to facilitate 
the solution it is necessary to assume a special azimuthal vari- 
ation of the surface reactance. For the model as chosen, 
strong resonance characteristics are again obtained. This 
model may be adapted to study the problem of a corrugated 
panel on a flat metallic ground plane which is excited by a 
parallel slot source. 
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